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EXECUTIVE  SUMMARY 


This  document  is  the  second  part  of  a  set  of  notes  intended  to  help  a  new  Ph.D.  in 
physics  or  chemistry  (or  person  of  equivalent  training)  become  introduced  to  the  specialized 
literature  in  the  applied  fields  of  atmospheric  nuclear  weapons  effects,  missile  reentry, 
upper  atmospheric  physics,  and  remote  sensing  through  the  atmosphere.  Thus,  the 
document  is  not  a  substitute  for  a  course  in  high-temperature  gas  dynamics  and  in 
atmospheric  physics,  but  it  presents  a  simple  overview  of  some  key  physical  factors  and 
gives  references  both  to  relevant  textbooks  and  to  the  "gray"  hterature  specific  to  the  field. 

Pan  I  of  these  notes  (IDA  Document  D-487,  May  1990)  consists  of  a  general 
introduction  (Cliapter  1),  a  discussion  of  equilibrium  air  at  temperatures  to  approximately 
10,000  K  (Chapter  2),  a  discussion  of  the  radiation  from  high-temperature  air  (Chapter  3), 
and-as  a  bridge  to  the  present  Pan  Il-a  brief  discussion  of  some  specific  mechanisms  for 
heating  the  earth's  atmosphere,  ranging  from  shock  heating  through  nuclear  explosions  to 
solar  ultraviolet  radiation  (C!hapter  4). 

The  present  Pan  11  of  the  notes  begins  with  Chapter  5,  which  gives  an  introduction 
to  reentry  physics.  This  begins  by  identifying  the  problem  of  reentry  heating,  including  a 
sketch  of  a  typical  ICBM  trajectory.  There  follows  a  discussion  of  a  normal  shock  in  air, 
of  various  commonly  used  approximations  to  the  equations  of  gas  dynamics  including 
boundary  layer  flow,  and  brief  treatments  of  the  laminar-turbulent  transition  including  the 
differences  between  aerodynamic  and  atmospheric  turbulence,  and  of  ablation.  The  chapter 
finishes  with  a  simple  physical  discussion  of  the  optical  signature  of  a  reentry  vehicle. 

Chapter  6  is  a  primer  on  nuclear  weapons  effects.  It  surveys  a  typical  weapon 
output,  and  both  low-  and  high-altitude  phenomenology  including  a  discussion  of  some 
high-altitude  multiburst  effects.  The  chapter  closes  with  a  discussion  of  some  systems 
implications  of  multiple  high-altitude  bursts  and  an  indication  of  the  quality  of  predictions 
in  this  field. 

Chapter  7  addresses  selected  topics  in  the  propagation  of  electromagnetic  radiation 
through  the  atmosphere.  After  a  discussion  of  molecular  absorption  in  the  atmosphere  and 
of  the  physics  of  the  LOWTRAN  and  FASCODE  computer  codes,  there  is  an  overview  of 
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atmospheric  visibility  and  aerosols  followed  by  a  discussion  of  cloud  physics  and  of  cloud- 
free  line-of-sight  (CFLOS)  analyses.  The  chapter  closes  with  a  brief  discussion  of 
atmospheric  turbulence  in  its  effects  on  electromagnetic  signal  propagation  and  a  sketch  of 
some  non-linear  effects  that  can  become  important  for  high-energy  laser  propagation. 


5.0  SOME  TOPICS  IN  REENTRY  PHYSICS 


5.1  THE  NATURE  OF  THE  PROBLEM 

An  ICBM  or  a  satellite  entering  the  earth’s  atmosphere  from  low  earth  orbit  has  an 
effective  entry  speed  of  approximately  7  km/sec,  or  Mach  Number  M  ~  20.  During  its 
descent  to  the  surface  of  the  earth  a  large  fraction  of  its  kinetic  energy  of  directed  motion 
will  be  transformed  into  heat  as  a  result  of  frictional  resistance  of  the  atmosphere.  The 
kinetic  energy  involved  is  ~  2  x  107  joule/kg,  which  is  an  order  of  magnitude  greater  than 
the  latent  heat  of  vaporization  of  water,  or  of  the  same  magnitude  as  the  vaporization  energy 
of  carbon,  which  is  one  of  the  highest  values  known  for  this  quantity.  Thus  the  problem  of 
designing  a  reentry  vehicle  which  would  dissipate  or  reject  this  amount  of  energy  is  a  very 
severe  one;  it  was  solved  in  the  1955-58  time  frame. 

The  initial  plan  was  to  let  the  vehicle  reenter  relatively  gradually  so  that  the  thermal 
energy  could  be  radiated  away.  Experimental  work  with  shock  tubes,  rocket  exhausts, 
etc.,  showed  that  the  radiant  heat  load  on  the  vehicle  was  not  nearly  as  severe  as  had  been 
thought  because  a  1.0-cm- thick  layer  of  hot  air  does  not  radiate  nearly  as  much  as  a 
blackbody  at  the  same  temperature  (see,  e.g..  Section  3.8  of  Part  I).  It  proved  possible  to 
use  an  ablative  heat  shield,  in  which  the  entry  vehicle  is  coated  with  materials  that 
decompose,  rejecting  large  quantities  of  energy  per  unit  mass  of  the  ablator. 

For  orientation,  an  ICBM  is  designed  for  high-speed  reentry.  Its  peak  heating 
occurs  near  20  km,  and  its  peak  deceleration  is  ~  50  g;  it  uses  an  ablative  heat  shield.  The 
Space  Shuttle  or  other  manned  entry  vehicle  is  limited  to  very  much  lower  deceleration;  the 
Shuttle  flies  a  lifting  reentry  with  peak  deceleration  ~  1  g.  Its  peak  heating  occurs  near  60 
km  and  it  dissipates  frictional  heat  by  radiation,  using  special  tiles  for  heat  protection. 

A  typical  reentry  vehicle  has  dimension  of  order  1.0  m,  so  that  it  has  a  significant 
encounter  with  the  atmosphere  only  at  densities  sufficiently  high  for  the  gas-kinetic  mean 
free  path  Igk  (defined  in  Section  2.5  of  Pan  I)  to  be  smaller  than  1.0  m,  which  requires 
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altitudes  below  ~  130  km.^  Thus  a  shock  gradually  develops  in  front  of  the  vehicle  as  it 
descends  below  this  altitude;  a  boundary  layer  (which  is  perhaps  1.0  cm  thick)  can  only 
develop  where  the  mean  free  path  is  much  less  than  the  boundary  layer  thickness,  or  at 
altitudes  below  ~  85  km.  Section  5.2  gives  a  brief  overview  of  some  aspects  of  reentiy 
phenomenology,  and  the  rest  of  this  Chapter  expands  on  that  discussion. 

Section  5.3  reviews  ICBM  trajectories  and  sketches  how  the  Space  Shuttle  differs 
from  this.  Section  5.4  gives  a  discussion  of  a  normal  shock  for  the  case  of  a  polytropic 
gas.2  Section  5.5  gives  a  broad  overview  of  hypersonic  aerodynamics,  introducing  the 
various  dimensionless  parameters  that  arc  used  in  this  field;  it  is  supplemented  by  Appendix 
E  which  derives  the  equations  of  gas  dynamics  (which  are  representations  of  the 
conservation  of  mass,  momentum,  and  energy)  from  the  Boltzmann  equation.  This  leads 
naturally  to  the  Euler  and  Navier-Stokes  equations,  including  the  empirical  characterization 
of  viscosity  and  heat  conductivity  coefficients  in  laminar  flow.  Section  5.6  sketches  an 
application  of  the  equations  of  fluid  mechanics  to  boundary  layer  flow,  with  a  qualitative 
indication  of  how  the  airflow  around  a  hypersonic  vehicle  can  be  represented.  As  the 
vehicle  enters  the  denser  atmosphere  the  air  flow  in  the  boundary  layer  becomes  turbulent, 
and  so  the  transition  from  laminar  to  turbulent  flow  is  taken  up  in  Section  5.7,  where  some 
elements  of  this  large  and  difficult  topic  arc  discussed.  The  key  point  is  that  frictional  heat 
transfer  is  much  larger  in  turbulent  than  in  laminar  flow.  Some  characteristics  of  ablation 
and  ablative  materials  used  for  thermal  protection  of  high-performance  reentry  vehicles  arc 
mentioned  in  Section  5.8. 

The  various  topics  addressed  here  in  cursory  fashion  arc  intended  to  provide  an 
indication  of  some  of  the  factors  that  have  to  be  considered  in  the  phenomenology  of 
reentry  vehicles.  The  final  section  of  this  chapter  gives  an  introduction  to  the  optical 
(visible  and  infrared)  signature  of  a  reentering  body.  The  discussion  is  qualitative  and 
closes  with  an  introduction  to  the  Optical  Signatures  Code  (OSC)  which  provides 
quantitative  answers  to  these  questions.  This  code  embodies  a  great  deal  of  experience 
from  reentry  physics,  in  effea  codifying  the  various  elements  that  have  been  discussed  here 
and  their  overall  application  to  missile  reentry. 


^  At  higher  altitudes  there  are  occasional  collisicms  which  slow  down  a  spacecraft  only  gradually  and 
cause  a  satellite  in  a  200-kni  orbit  to  de-orbit  in  a  time  of  order  one  year. 

^  A  polytropic  gas  is  one  with  y  =  Cp/Cv  =  constant.  This  docs  not  apply  to  air  under  reentry  conditions  i 

because  of  the  dissociation  of  02"See  Part  I,  Section  2.4. 
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5.2  OVERVIEW  OF  REENTRY  PHENOMENOLOGY 


When  an  ICBM  or  low-earth-orbiting  satellite  traveling  at  hypersonic  speed^  enters 
the  earth's  atmosphere,  it  is  slowed  down  by  collision  with  air  molecules.  This  produces  a 
hypersonic  shock  wave  that  stands  off  in  front  of  the  vehicle.  At  still  lower  altitudes  (or 
higher  densities)  a  boundary  layer  is  produced  on  the  vehicle. 

During  the  slowing  down  of  a  vehicle  in  atmospheric  entry,  its  kinetic  energy  is 
transformed  into  thermal  energy  of  the  surrounding  air,  and  some  of  this  thermal  energy 
reaches  the  surface  of  the  vehicle.  Figure  5.1a  shows  schematically  how  the  atmospheric 
density,  velocity,  deceleration  and  heating  rate  change  with  altitude,  indicating  how,  as  the 
vehicle  enters  the  denser  lower  atmosphere,  its  velocity  gradually  decreases,  with  peaks  in 
deceleration  and  in  heating  rate  at  intermediate  altitudes  (30-60  km).  Figures  5.  lb  and  5.  Ic 
demonstrate  in  semi-quantitative  form  how  energy  is  converted.  In  Figs.  5.1b  and  5.  Ic  we 
show  the  energy  conversion  fraction,  which  is  that  fi'action  of  the  kinetic  energy  lost  by  the 
vehicle  which  reaches  the  surface,  as  a  function  of  altitude. 

At  very  high  altitudes  in  the  free  molecular flow  regime  where  the  gas  kinetic  mean 
free  path  Igk  is  not  small  compared  to  critical  vehicle  dimensions,  as  much  as  half  of  the 
kinetic  energy  lost  by  the  vehicle  appears  as  heat  input  to  the  body.  At  lower  altitudes  in 
the  continuum  regime,'*  a  boundary  layer  of  hot  air  develops  around  the  reentry  vehicle  and 
transfers  thermal  energy  to  the  vehicle,  mainly  by  conduction  and  convection  at  speeds  of 
around  7  km/sec.  Note  that  while  the  energy  conversion  fraction  here  may  be  only  0.01 
because  most  of  the  energy  goes  into  the  air,  yet  here  the  heat  load  on  the  vehicle  will  be 
much  greater  than  in  the  free  molecular  flow  regime  because  the  total  energy  transferred  is 
so  much  larger. 

Figure  5.1b  shows  a  representative  variation  of  energy  conversion  fraction  with 
altitude  for  a  relatively  blunt  vehicle  traveling  at  ICBM  velocity  of  about  7  km/sec.  Note 
the  transition  from  laminar  to  turbulent  convection  in  the  SO-km  altitude  range. 

If  a  vehicle  travels  at  ICBM  velocity  (or  low  earth-orbiting  velocity,  which  is  only 
about  10  percent  higher-cf.,  e.g.,  Gazley,  1961),  when  it  transfers  all  its  kinetic  energy 


^  Supersonic  means  that  the  vehicle  velocity  v  is  greater  than  the  local  sound  speed  a;  hypersonic  means 
that  the  velocity  v  is  very  much  greater  than  the  local  sound  speed  a.  For  low  earth  orbit  and  ballistic 
missile  reentry  we  typically  have  a  Mach  number  M  *  v/a  -  20. 

^  This  is  said  to  occur  when  the  gas-kinetic  mean  free  path  Ig^  is  very  much  smaller  than  a  critical 
vehicle  dimension. 
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into  theraial  energy  of  air,  the  peak  air  temperature  ("stagnation  temperature"  defined  in 
Section  5.4)  is  in  the  range  6(X)0  -  80(X)  K,  (depending  on  pressure)  with  O2  largely  and 
N2  partly  dissociated.  For  super-orbital  or  near  escape  velocity  of  1 1  km/sec  the  enthalpy 
of  air  is  increased  by  a  factor  (1 1/7)2  «,  2.5,  and  the  stagnation  temperature  is  somewhat 
higher  (8000-10,000  K)  but  now  radiative  heat  transfer  (which  goes  as  rather  than  as 
grad  T)  becomes  important  at  the  lower  altitudes,  in  addition  to  conductive  and  convective 
heat  transfer.  Figure  5.1c  shows  schematically  how  gaseous  radiation  effects  become 
important  at  low  altitudes  for  these  very  high  velocities,  changing  the  energy  conversion 
fraction. 

Figure  5.2  gives  an  overview  of  the  air  flow  around  a  non-lifting  (i.e.,  ballistic) 
reentry  vehicle.  It  demonstrates: 

1 .  The  shock  wave  A  associated  with  the  encounter  of  the  vehicle  moving 
through  the  ambient  air,  which  is  thereby  displaced  non-adiabatically  by  the 
vehicle.  Except  at  the  front  stagnation  point  B  where  the  shock  is  normal,  this 
is  an  oblique  shock  which  extends  downstream,  typically  for  many  vehicle 
lengths,  and  then  gradually  fades  out. 

2 .  The  stagnation  point  B  is  that  region  in  the  flow  in  which  all  the  kinetic  energy 
of  directed  motion  is  transformed  into  thermal  energy,  so  that  the  maximum 
temperature  in  the  flow  or  stagnation  temperature  is  attained  here.^ 

3 .  The  boundary  layer  C,  which  begins  in  the  stagnation  region  at  the  front  of  the 
vehicle  and  then  gradually  grows  as  it  moves  downstream  along  the  body.  In 
hypersonic  flow  this  is  a  region  close  to  a  solid  body  in  which  the  ambient  air 
interacts  with  the  solid  boundary  through  collisions,  certainly  through  firictional 
(viscous)  interactions  and  sometimes  also  through  ablation  or  other  chemical 
reactions.  Within  the  boundary  layer  the  flow  is  generaUy  subsonic  and  can  be 
treated  as  incompressible  but  viscous— i.e.,  collisional  interactions  with  the 
body  are  imponant.  Outside  the  boundary  layer  the  flow  is  treated  as 
supersonic,  i.e.,  compressible,  so  that  the  flow  equations  are  non-linear  and 
thus  difficult  to  solve,  but  the  detailed  interaction  of  the  flow  with  the  body  is 
not  important. 

4.  Past  the  downstream  end  of  the  body  the  flow  has  to  revert  gradually  to  the 
original  undisturbed  flow.  This  is  a  complex  process: 


^  See  Fig.  5.3,  which  shows  that  the  peak  or  stagnation  temperature  is  reached  at  the  blunt  nose  of  the 
RV.  Table  5.2  (below)  gives  the  energy  balance  behind  a  normal  shock  in  air.  and  shows  that  for 
ICBM  or  even  IRBM  entry  (Uq  =  7  km/sec  and  5  km/sec  respectively)  more  energy  goes  into  O2 
dissociation  than  into  any  other  degrees  of  freedom. 


4.1  Behind  the  body  the  boundary  layer  continues  as  the  shear  layer  which 
thickens  gradually  but  shrinks  down  to  form  a  minimum  width  or  neck  D 
before  expanding  gradually  to  form  a  wake  E  which  extends  downstream 
for  many  vehicle  lengths  before  gradually  fading  out  like  the  bow  shock 
A. 

4.2  Between  the  end  of  the  body  and  the  neck,  and  inside  the  shear  layer  there 
is  the  recirculation  region  F  which  is  filled  by  air  flowing  out  of  the 
shear/boundary  layer  proper,  the  character  of  the  flow  is  shown  in 
Fig.  5.2. 

4.3  Note  that  near  the  rear  of  the  recirculation  region  there  is  a  rear  stagnation 
point  G  where  the  temperature  becomes  almost  as  large  as  at  the  front 
stagnation  point.^ 

4.4  The  wake  itself  produces  a  trailing  shock. 

5 .  Note  that  there  is  a  PrandtUMeyer  expansion  fan  H  associated  with  a  turning  of 
the  boundary  layer  past  the  end  of  the  body. 

6.  There  may  also  be  a  secondary  shock  near  the  front  of  the  vehicle  where  the 
shape  changes  from  hemispherical  to  conical,  and  for  vehicles  of  more 
complex  shape  (e.g.,  sphere-cylinder-conical  flare)  there  will  be  similar  shocks 
and  expansion  fans  associated  with  every  turning  of  the  flow. 

7.  If  one  asks  for  the  mean  temperature  of  the  flow  field  surrounding  the  body, 
there  is  a  maximum  at  the  stagnation  point  (typically  7500  K  for  ICBM 
reentry);  in  the  boundary  layer  along  the  body  the  temperature  is  very  much 
less-typically  2000-3000  K  for  an  ablating  vehicle.^  Past  the  rear  of  the  body 
the  temperature  varies  with  location,  with  a  maximum  at  the  rear  stagnation 
point,  but  in  the  wake  proper  beyond  the  neck  the  temperature  gradually 
decreases  to  its  ambient  value. 

Figure  5.3  shows  the  temperature  profile  at  various  stations  along  a  "high- 
performance"  ablatively  cooled  entry  vehicle  at  different  altitudes.  At  very  high  altitudes 
the  peak  temperature  develops  first  on  the  rounded  (i.e.,  blunt)  nose,  and  then  gradually,  as 
the  boundary  layer  builds  up,  so  does  the  frictional  heating  along  the  side  of  the  body. 
This  heating  (which  at  flow  velocities  below  approximately  10  km/sec  is  due  mainly  to 
conduction  and  convection  rather  than  to  radiation)  increases  down  to  the  altitude  of  peak 
heating  which  here  occurs  at  70,000  ft  (21  km).  The  heating  rate  then  declines.  Note  the 


^  The  tempetauire  here  is  not  quite  as  large  as  at  the  front  stagnation  point  because  of  various  energy 
losses  in  the  flow. 

^  Figure  S.3  gives  some  representative  temperature  profiles  along  an  ablative  reentry  vehicle  at  different 
altitudes. 
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sudden  increase  in  temperature  at  station  36  in.  at  90,000  ft.  altitude.  This  is  due  to  the 
transition  from  laminar  to  turbulent  flow.  The  precise  character  of  transition  varies— 
Fig.  5.3  comes  from  a  computer  calculation. 

The  discussion  centering  on  Figs.  5. 1-5.3  demonstrates  many  of  the  significant 
physical  factors  that  contribute  to  reentry  phenomenology,  but  it  is  evidently  not  an 
adequate  account  of  hypersonic  aerodynamics.  However,  it  does  introduce  and  set  in 
context  a  number  of  the  physical  concepts  that  are  discussed  later,  in  particular  hypersonic 
shocks  (Section  5.4),  the  characteristics  of  the  equations  of  fluid  mechanics  and  their 
approximations  (Section  5.5),  the  concept  of  a  boundary  layer  (Section  5.6),  the  laminar- 
turbulent  transition  (Section  5.7),  and  the  significance  of  ablation  (Section  5.8). 

Note  that,  e.g.,  Anderson  (1986)  indicates  that  there  are  five  distinct  effects 
characterizing  hypersonic  flows: 

1 .  There  is  a  shock  wave  on  the  body  close  to  its  surface,  i.e.,  a  thin  shock  layer. 

2.  Viscous  effects  are  strong  and  can  dominate  the  flow  field  by  thickening  the 
boundary  layer  and  merging  the  shock  into  the  viscous  layer. 

3.  At  low  densities  there  can  be  temperature  and  velocity  discontinuities  ("slip 
conditions")  at  the  surface  of  the  body. 

4.  Hypersonic  flows  are  high-energy  flows,  leading  to  non-ideal  behavior  of  the 
gas  flowing  around  the  body.  The  principal  effects  are  vibrational  excitation 
and  dissociation  which  lead  to  a  change  in  specific  heats,  in  specific  heat  ratio 
Y  =  Cp/Cv,  and  also  in  chemical  composition. 

5.  In  hypersonic  flow  over  a  blunt-nosed  body,  the  shock  wave  around  the  body 
is  highly  curved  and  there  is  a  vorticity  (angular  momentum)  interaction. 

Items  1-3  and  5  are  discussed  in  the  aerodynamics  literature,  see,  e.g.,  Anderson 
1982, 1984, 1986,  while  item  4  is  discussed  in  Chapter  2  of  Part  I. 
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Altitude.  1.000  ft 


Density  Velocity  Dcccler-  Heoting 
p  u  otion  rote 


Energy  conversion  fraction  Energy  conversion  fraction 

(b)  (c) 


Figure  5.1.  Energy  Transfer  on  Reentry,  (a)  Schematic  of  Slowing-Down  and 
Heating-Up  as  Function  of  Altitude,  (b)  Energy  Conversion  for  ICBM 
Reentry  (u  -  7  km/sec)  (c)  Energy  Conversion  for  Ultra-Hlgh-Speed 
Reentry  (u  >  11  km/sec):  Effects  of  Hot  Gas  Radiation. 

(Source:  Gazley,  1961) 
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Figure  5.2. 


Schematic  of  Air  Fiow  Around  a  Reentry  Vehicie 
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11*10-M-S 


AXIAL  LOCATION,  in. 


Figure  5.3.  Surface  Temperature  Profiles  on  a  Quartz  Phenolic 
Ablative  Reentry  Vehicle  at  Different  Altitudes 
(from  R.  Bartle,  AVCO  Systems  Division) 
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Figure  5.3. 


(continued) 
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5 . 3  ICBM  TRAJECTORIES 


Figure  5.4  shows  some  representative  trajectories  for  an  Intercontinental  Ballistic 
Missile  (ICBM)  of  range  9000  km.  It  shows  the  altitude,  with  time  markers  every  30  sec 
from  burnout,  for  four  different  trajectories.  These  are  customarily  characterized  by  the 
reentry  angle,  which  is  defined  as  the  path  angle  to  the  horizontal  at  400,000  ft  (122  km) 
altitude  during  reentry  or  at  booster  burnout.  These  both  typically  occur  near  400,000  ft 
altitude. 

A  normal  or  minimum  energy  trajectory  corresponds  roughly  to  24-deg  reentry 
angle.  A  depressed  trajectory  (15-deg  reentry  angle)  would  make  the  ICBM  harder  to 
detect  from  space  because  of  its  lower  maximum  height;  more  energy  is  needed  to 
overcome  frictional  heating  in  the  atmosphere.  A  lofted  trajectory  (35-40  deg  reentry  angle) 
requires  more  energy  to  send  the  ICBM  higher  above  the  earth's  surface.  The  range  of 
reentry  angles,  15  deg  to  40  deg,  roughly  covers  the  range  considered. 

Table  5.1  gives  some  details  of  the  trajectory  for  a  "representative"  10,000-km, 
minimum-energy  trajectory.  It  is  shown  here  to  let  the  reader  see  how  various  parameters 
change  along  the  trajectory.  As  a  function  of  time,  we  show  the  horizontal  range,  altitude, 
velocity,  path  angle  y  to  the  horizontal,  and  a  representative  stagnation  temperature*  Ti 
along  the  trajectory.  In  the  reentry  phase  we  also  show  the  large  (negative)  acceleration  and 
the  ballistic  coefficient  |3  ^ 

P  =  W/CdA  ,  (5.1) 

where  W  =  mass  of  vehicle  ,  A  =  projected  area,  and  Cd  =  drag  coefficient,  P  changes 
slightly  along  the  entry  trajectory,  principally  because  Cd  changes  with  velocity,  including 
effects  of  the  laminar-to-turbulent  transition. 

For  a  representative  ICBM  trajectory  with  a  23-deg  reentry  angle,  peak  heating 
(for  a  vehicle  with  P  ~  2000  Ib/ft^)  occurs  at  23  km  (75,000  ft)  with  a  peak  deceleration  of 


*  The  values  of  Ti  shown  here  correspond  to  a  particular  value  of  spherical  nose  radius.  They  assume 
equilibrium  radiation  and  no  thermal  inertia,  so  that  the  numerical  values  are  representative  rather  than 
accurate.  They  do  display  exit  heating  as  well  as  entry  heating.  An  ambient  value  T  =  300  K  is 
assumed,  so  th^  significant  heating  only  occurs  for  Ti  different  from  300  K. 

^  P  is  normally  measured  in  Ib/ft^;  2000  Ib/ft^  =  9774  kg/m^  is  a  representative  value  for  the  present 
generation  of  ablatively  cooled  reentry  vehicles.  (For  a  satellite,  values  of  50-200  kg/m^  are  typical, 
giving  rise  to  lifetimes  of  order  several  years  in  400-500  km  orbits.) 

Cd  =  2  for  a  flat  plate  normal  to  the  flow,  and  Cd  =  2  sin^  for  a  cone  of  half-angle  0. 
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50-60  g.  By  contrast,  the  first  generation  of  U.S.  reentry  vehicles  (1957-1960  time  frame) 
used  a  copper  heat  sink  with  3  ~  100  Ib/ft^,  so  that  peak  heating  occurred  near  50  km 
(150  kft)  and  extended  over  a  much  longer  time  so  that  radiative  heat  transfer  was  adequate 
to  dissipate  the  slower  rate  of  heat  generation;  the  peak  deceleration  was  of  order  10  g.  The 
first  U.S.  manned  space  vehicle  (Mercury)  used  a  somewhat  similar  reentry  technology, 
but  all  subsequent  U.S.  manned  space  vehicles  have  flown  a  lifting  reentry  which  can  be 
much  gentler  than  a  ballistic  trajectory.  Thus  the  Space  Shuttle  reentry  angle  is 
approximately  1  deg  and  its  peak  heating  occurs  near  61  km  (200,000  ft);  its  peak 
deceleration  is  ~  1  g.  The  Shuttle,  with  its  much  slower  reentry  and  thus  lower  heating  rate 
than  an  ICBM,  uses  radiative  cooling  and  is  covered  with  tiles  which  dissipate  excess 
thermal  energy  by  radiation,  unlike  ballistic  missiles  which  have  ablative  heat  shields. 

5.4  NORMAL  SHOCK 

5.4.1  Introduction 

A  normal  shock  is  the  most  violent  elementary  discontinuity  that  one  studies  in 
aerodynamics.  It  is  a  one-dimensional  motion  which  is  observed  when  a  supersonic 
gas  flow  strikes  a  flat  plate  located  perpendicular  to  the  flow  direction;  or  equally,  when  a 
flat  plate  (piston)  moves  supersonically  into  a  stationary  gas.  The  flow  behind  the 
discontinuity  is  always  subsonic.  A  shock  is  adiabatic  but  non-isentropic. 

We  use  the  following  notation: 

upstrcam/supersonic  =  subscript  1 
downstream/subsonic  =  subscript  2 
p  =  density 

u  =  flow  velocity 

a  =  sound  speed 

e  =  internal  energy  per  gm  =  (Cv/M)dT 

h  =  e  +  p/p  =  enthalpy  per  gm  =  (Cp/M)dT 


M 


molecular  weight. 


ICBM  Tralectorles  (Source:  APS,  1987) 


Table  5.1.  A  Representative  ICBM  Trajectory 


Time(sec) 

Range 

Altitude 

Velocity  Acceleration 

P 

Path  angle 

T, 

(sec) 

(km) 

(km) 

(km/sec) 

(9) 

(kg/m^) 

(deg) 

(K) 

0 

0 

0 

0 

90 

300 

43 

5.1 

10.2 

0.61 

54.3 

970 

60 

13.8 

20.6 

1.01 

46.7 

1180 

74 

26.0 

32.4 

1.45 

42.2 

1240  (M) 

90 

47.4 

50.7 

2.12 

38.5 

1200 

105 

77.8 

74 

2.96 

36.2 

1 060  Stage  1  Burnout 

121 

116 

101 

3.06 

34.2 

475 

150 

194 

151 

3.44 

30.8 

300 

210 

398 

263 

4.64 

25.5 

263.5 

656 

385 

6.70 

23.4 

300  Stage  2  Burnout 

463 

1740 

857 

6.13 

19.6 

Stage  3  Burnout 

1176 

4930 

1580 

5.77 

0.0 

Apogee 

1890 

8120 

855 

6.12 

•19.7 

2100 

9260 

355 

6.75 

-23.6 

300 

2180 

9750 

129 

7.05 

1610 

-24.7 

300 

2182 

9770 

120 

7.06 

523 

•24.8 

540 

2193 

9830 

90.2 

7.10 

4070 

-24.9 

1480 

2203 

9890 

60.1 

7.14 

6730 

-25.0 

2950 

2213 

9960 

29.9 

7.04 

6.0 

7330  (M)  -25.2 

4950 

2218 

9990 

15.7 

6.01 

44 

7240 

-25.3 

5790(M) 

2220 

10000 

11.0 

4.94 

63  (M) 

7200 

-25.4 

5440 

2230 

10020 

0 

0.96 

17 

3380 

-27.6 

1750 

*  Data  kindly  provided  by  R.E.Finke,  IDA/STD. 
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Wc  consider  a  poly  tropic  gas  for  which 

Cp/Cv  =  constant  =  y=  1  +  1/s  ,  (5.2) 

where  s  =  1/(Y-1).  For  an  ideal  monatomic  gas,  Y=  5/3,  for  an  ideal  diatomic  gas  y=  7/5, 
but  for  dissociating  air  (i.e.,  at  temperatures  in  the  range  2000-6000  K),  y  '  ~  115  (see. 


e.g.,  Gilmore,  1955,  p.  67).  For  a  polytropic  gas 

p/p  =  RT/M  ;  e  =  sRT/M ;  h  =  (1  -i-  s)  RT/M  .  (5.3) 

Indeed,  e  ~  T  is  the  definition  of  a  polytropic  gas. 

5.4.2  Conservation  Laws 

The  conservation  laws  for  mass,  momentum,  and  energy  are 

Pi  ui  =  p2  U2  =  m  (constant)  (5.4) 

pi  +  Pi  ui^  =  p2  +  P2  02^  =  P  (constant)  (5.5) 

hi-»-(i/2)ui2  =  h2-i-(l/2)u22  =  [(l/2)(Y+l)/(Y-l)]a*2  ,  (5.6) 

where 

a*2  =  ui  U2  .  (5.7) 

5.4.3  Jump  Conditions 

It  is  convenient  to  introduce  a  density  ratio 

e  =  Pl/p2  =  U2/ui  ,  (5.8) 


which  characterizes  the  strength  of  the  shock.  Given  conditions  upstream  of  the  shock 
(ui,Ti,  pi,  or  pi)  and  the  equation  of  state,  then  when  the  density  ratio  e  is  specified,  so 
are  conditions  downstream  of  the  shock.  In  particular,  there  is  a  Hugoniot  condition 

h2-hi  =  (l/2)ui2(l-£2)  (5.9) 

or 

P2-Pl=Plui2(l-E)  .  (5.10) 

5.4.4  Some  Results  for  Polytropic  Gases 

These  results  are  specific  for  polytropic  gases,  but  comparable  results  apply  to  air-- 
the  Mollier  chart-see,  e.g.,  Feldman  (1957),  Logan  and  Treanor  (1957).  We  have 
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a*2/aj2  r=  A  +  B  Mi^  where  Mi  =  ui/ai  (5.11a) 

and  A  =  2sB  ;  B  =  l/(2s+ 1)  (5.11b) 

e  =  A/Mi2  +  B  (5.11c) 

T2/T1  =  l+Mi2(l-e2)/2s  (5.1  Id) 

P2/P1  =  (l/e)(T2/ri)  .  (5.11e) 

Notice  that  as  Mi  ->  oo ,  e  ->  B,  and  we  have  the  following: 

Monatomic  ideal  gas:  y  =  1 .67  ,  s  =  3/2  ;  B  =  1/4  (5. 12a) 

Diatomic,  polytropic  gas  (e.g.,  undissociated  air):  y  =  1.4  ,  s  =  5/2  ;  B  =  1/6  (5.12b) 

Dissociated  air  [very  rough,  see  Gilmore  (1955),  p.  67]: 

y  ’=  1.15  ;  s  =  20/3  ;  B  =  0.07  .  (5.12c) 


Thus  the  density  ratio  e  across  the  shock  remains  finite  but  the  temperature  ratio  can 
become  very  large:  T2/T1  ~  Mi^  — » 00  as  Mi  — >  00 ,  so  that  large  flow  velocities  give  very 
high  "stagnation"  temperatures  and  pressures,  but  only  moderate  increases  in  density. 

5.4.5  Energy  Transfer  and  Non-Equilibrium  Conditions 

A  shock  is  in  fact  not  an  instantaneous  discontinuity,  but  has  some  finite  thickness 
and  structure.  This  is  determined  by  the  transfer  of  energy  from  initial  kinetic  energy  into 
various  degrees  of  freedom.  Figure  4.2(c)  in  Part  I  shows  the  effects  of  kinetics: 
translational  and  rotational  equilibration  occurs  very  rapidly,  giving  a  high  translational 
temperature  immediately  across  the  shock  as  all  the  directed  energy  of  motion  is 
transformed  into  random  translational  and  rotational  energy.  Transfer  into  vibrational 
excitation  takes  longer  and  leads  to  a  lower  temperature  as  the  energy  is  distributed  among 
more  degrees  of  freedom.  Dissociation  and  electronic  excitation  take  longer  still,  giving  yet 
a  further  reduction  in  kinetic  temperature  while  ionization  is  slower  still,  yielding  the  final 
equilibrium  temperature.  Table  5.2  shows  the  actual  energy  distribution  for  a  normal  shock 
in  air  at  5-7  km/sec  where  most  of  the  energy  goes  into  dissociation,  a  significant  amount 
into  vibrational  excitation,  and  very  little  into  ionization.  By  contrast,  for  a  normal  shock  in 
Argon  most  of  the  energy  goes  into  ionization-see,  e.g.,  Bauer,  1969. 


I 
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Table  5.2.  Energy  Balance  Behind  a  Normal  Shock  In  Air  (Source:  Bauer,  1969) 


(A  velocity  of  5  km/sec  corresponds  to  IRBM  entry  while  7  knVsec  corresponds  to  ICBM  entry. 
The  ambient  conditions  are  po  >  0.01  atm.,  Tq  >  294  K, 
and  energies  are  listed  in  eV/molecule  or  atom.) 


Uo 

5  km/sec 

7  km/sec 

Upstream  of  the  shock 

Directed  kinetic  energy 

Thermal  enthalpy  CpT 

3.76 

0.006 

7.3 

0.006 

Downstream  of  the  shock 

Stagnation  temperature  (K) 

5800  (~  0.5  eV) 

7500  (~  0.65  eV) 

Enthalpy  of  translation  (3/2)  kT 

1.25 

1.62 

Rotational  enthalpy,  kT 

0.5 

0.65 

Vibrational  enthalpy 

0.37 

0.51 

Dissociation  energy 

1.51 

2.72 

Electronic  energy 

0.011 

0.075 

Ionization  energy 

0.002 

0.006 

Chemical  energy  (N2  +  O2  <«>  2NO) 

-0.02 

-0.01 

Equilibrium  composition  downstream  of  the  shock 

IN2] 

0.579 

0.304 

[NO] 

0.0164 

0.00615 

[O2] 

0.000974 

0.000124 

(Nj 

0.0840 

0.418 

(0) 

0.319 

0.272 

[e-l-[NO-^] 

0.000137 

0.000630 
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5.4.6  Numerical  Results  for  the  Stagnation  Temperature 


Some  results  for  T2  (  for  Ti  =  294  K  so  that  ai  =  343.8  m/sec)  are  shown  in 
Table  5.3  where  we  compare  the  results  for  an  ideal  monatomic  gas  (see  Eq.  5.12a),  a 
diatomic  non-dissociating  gas  (see  Eq.  5.12b)  and  calculations  for  dissociating  air  with 
y  =  1.15  (see  Eq.  5.12c)  with  the  numerical  calculations  for  air  and  Argon. 


Table  5.3.  Stagnation  Temperature  Computed  Under  Various  Conditions 
(Ti  =  294  K,  l.e..  ai  =  343.8  m/sec) 


ui  (km/sec) 

Mi  =  ui/ai 

5 

14.5 

7 

20.4 

10.7 

31.1 

Stagnation  Temperature  (K): 

Ideal  monatomic  gas  (Eq.  5.12a) 

19,600 

38,500 

89,300 

Argon  (see  Bauer,  1969) 

13,000 

15,700 

- 

Undissociated  air  (Eq.  5.12b) 

12,300 

24,100 

55,700 

Air  (see  Table  5.2) 

5800 

7500 

- 

Dissociated  polytropic  air  (Eq.  5.12c) 

4900 

9400 

21.600 

Note  that  the  results  for  dissociated  air  depend  on  density  and  are  only  schematic, 
because  for  dissociating  air  enthalpy  e  depends  on  both  p  and  T  in  a  complicated  way.  The 
value  of  Y '  is  not  a  constant  One  sees  that  in  the  case  of  dissociating  air  a  large  amount  of 
energy  goes  into  breaking  chemical  bonds  while  in  Argon  a  large  amount  of  energy  goes 
into  ionization.  Table  5.2  shows  that  the  other  available  internal  degrees  of  freedom  also 
take  up  significant  amounts  of  energy— evidently,  the  more  the  internal  energy,  the  lower 
the  stagnation  temperature. 

5.5  EQUATIONS  OF  GAS  DYNAMICS 

To  solve  fluid  mechanical  flow  problems,  one  has  to  write  down  the  equations  of 
motion  which  correspond  to  the  conservation  of  mass,  momentum,  and  energy,  together 
with  appropriate  collisional  (dissipative)  terms,  and  terms  arising  from  chemical  reactions 
(mainly  dissociation)  and  from  ablation.  These  are  five  equations  in  the  five  variables 
(three  velocity  components,  density  and  pressure),  but  they  are  a  set  of  coupled  non-linear 
partial  differential  equations  which  are  too  difficult  to  solve  without  approximation  using 
present  computational  power:  certainly  in  the  1930s  to  1960s  time  frame  when  most 
basic  work  in  this  field  was  done  it  was  necessary  to  make  some  very  significant 
approximations. 
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Appendix  E  sketches  how  one  can  derive  the  various  approximations,  in  particular 
the  Euler  and  the  Navier-Stokes  equations,  from  the  Boltzmann  equation,  and  also 
gives  simple  expressions  for  the  gas  kinetic  transport  coefficients  (viscosity,  thermal 
conductivity,  and  diffusion  coefficient).  Here  we  give  an  overview  of  the  physical  content 
of  these  equations.  The  Euler  equations  correspond  to  a  perfect  fluid  which  cannot  sustain 
tangential  (shearing)  forces,  only  normal  (pressure)  forces.  They  are  useful  under  many 
conditions,  but  they  do  not  work  near  surfaces  because  they  cannot  explain  the  drag  force 
on  a  surface  in  a  fluid  (because  they  do  not  have  a  viscous  drag  term).  The  Navier-Stokes 
equations  do  have  viscous  and  thermal  conduction  terms,  and  thus  are  more  generally 
applicable  but  correspondingly  more  difficult  to  evaluate.  In  the  boundary  layer,  the  flow 
interacts  with  a  surface,  so  that  the  viscous  terms  are  important  and  one  has  to 
use  the  Navier-Stokes  equations— but  because  the  flow  is  subsonic,  it  can  be  treated  as 
incompressible,  which  is  a  great  simplification  that  will  now  be  demonstrated. 

Compressibility  is  a  measure  of  the  change  of  volume  AV  under  an  external 
force,  here  the  pressure  p.  Thus,  we  may  define  a  modulus  of  elasticity  E  by  the  equation 

Ap  =  -E(AVA^o)  •  (5.13) 

For  a  gas,  E  =  po,  the  initial  pressure,  and  the  relation  between  changes  in  density  p  and 
volume  V  is 

Ap/po  =  -AVA^o  .  (5.14) 

so  that 

Ap/po  =  Ap/po  ,  (5.15) 

i.e.,  the  flow  can  be  treated  as  incompressible  if  Ap/po  «  1.  Now  from  Bernoulli's 
equation 

p  +  (1/2)  p  u2  =  constant  ,  (5. 16) 

where  u  =  flow  velocity,  we  find 

Ap/po  =  (1/2)  Po  u2  /  E  =  (1/2)  Po  u2/po  =  (1/2)  (u/a)2  .  (5.17) 

so  that  the  effects  of  compressibility  can  be  ignored  if  M2/2  «1,  where  M  =  u/a  is  the 
Mach  number. 

Traditionally,  aerodynamicists  use  similarity  solutions  to  a  great  extent,  with  a 
number  of  dimensionless  parameters.  The  common  dimensionless  parameters  are  listed  in 
Table  5.4  together  with  an  indication  of  their  physical  meaning.  The  most  important  ones 


1 1  Analysis  after  Schlichting,  1979,  p.  9. 
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are  the  Mach  and  Reynolds  numbers.  Note  that  the  Reynolds  Number  and  also  the 
Knudsen  number  contain  a  characteristic  length.  The  choice  of  characteristic  length 
depends  on  the  problem  at  hand;  it  may  be  distance  along  or  perpendicular  to  the  flow. 
This  will  be  discussed  further  in  Section  5.5.1. 

Note  that  many  aerodynamics  calculations  are  done  assuming  Le  =  Pr  =  Sc  =  1, 
which  tends  to  simplify  the  equations  to  be  solved. 


Table  5.4.  Common  Aerodynamic  Dimensionless  Parameters 


Mach  Number 

M  -u/a 

The  physical  meaning  of  is  that  it  compares  flow 
energy  to  internal  energy,  so  that  when  M«1  the 
thermophysics  is  not  coupled  to  the  flow 

Reynolds  Number 

Rex  =  pux/p  «=  (Inertial  force,  pu  3u/dx)  +  (Frictional  force,  p. 

where  x  =  distance  along  the  flow , 
y  >  distance  perpendicular  to  flow  (applies  to 
boundary  layer) 

Lewis  Number 

Leij  «=  Djj  p  cp/k 

Compares  diffusive  transport  of  energy  with  heat 

conduction 

Prandtl  Number 

Pr«  pcp/k 

Compares  momentum  and  heat  conductbn 
transport 

Schmidt  Number 

Sc  -  p  /p  D 

Compares  diffusive  and  nfX)mentum  transport 

Damkoehler  No.(1) 

Dami  -  x/u  tchem  (tchem  *  characteristic  time  for  chemical  reaction) 

Damkoehler  No.(2) 

Dam2  -  q/Cp  To  (q 

-  heat  added/gm;  Tq  >  characteristb  temperature) 

Knudsen  Number 

Kn  -  Igk/d 

(d  «  characteristb  distance  for  fbw;  Kn  >1  at  very  bw 
densities/  high  altitudes-implies  failure  of  continuum 
concept ) 

Richardson  Number 

Ri  *  -  (g/r)(3r/dz)  /  (3u/9z)2  (g  «  gravitational  acceleratbn,  z  =  vertical 

distance).  Ri  is  the  ratb  of  an  atmospheric  buoyancy 
term  that  may  promote  either  stability  or  instability, 
and  a  wind  shear  term  that  always  tends  to  produce 
turbulence.  Ri  is  important  in  atmospheric  fbws  such 
as  wind,  cumulus  convectbn,  etc. 

i 
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5.6  BOUNDARY  LAYER  FLOW 


Consider  the  flow  of  a  fluid  parallel  to  a  flat  plate.  Normally  there  is  no  slip  along 
the  surface  of  the  plate,  i.e.,  the  parallel  velocity  next  the  surface  is  zero,  while  of  course 
far  from  the  surface  of  the  plate  it  has  the  ambient  value  u(oo).  Thus,  the  velocity  profile 
parallel  to  the  surface  has  the  general  character  shown  in  Fig.  5.5  . 

The  fact  that  u  differs  from  u(oo)  close  to  the  wall  cannot  be  explained  by  the  Euler 
equations  (which  do  not  take  into  account  molecular  collisions  of  any  kind,  either  gas-gas 
or  gas-wall).  However,  the  Navier-Stokes  equations  which  contain  a  viscosity  term  to 
account  for  molecular  collisions  can  explain  how  the  fluid  can  support  a  shear  stress  and 
thus  must  be  used  within  the  boundary  layer. 


->x 


UU) 
— V 


PreeSpMc 


LmiifMr  B<iund«7  Layer 


Figure  5.5.  General  Characteristics  of  Boundary  Layers 


5-21 


The  thickness  5  =  5(x)  of  the  boundary  layer  as  a  function  of  distance  x  along  the 
surface  can  be  explained  in  terms  of  the  rough  equality  of  frictional  and  inertia  forces— see 
the  discussion  of  Reynolds  number  Re  in  Table  5.4— from  which  we  write 


Inertia  force,  p  u  du/Bx  ~  p  u^/x 

(5.18a) 

Friction  force,  |i  d^u/dy^  ~  p  u/52  , 

(5.18b) 

and  equating  these  terms  gives 

6  =  a7  (px/pu)  , 

(5.19) 

where  A  is  a  numerical  coefficient  whose  value  for  a  laminar  boundary  layer  is 
approximately  equal  to  5  (cf.  Schlichting,  1979,  p.  76.) 

When  the  Reynolds  number  becomes  sufficiently  large  (see  Section  5.7  below)  the 
flow  becomes  turbulent,  and  then  the  flow  velocity  profile  changes.  A  turbulent  boundary 
layer  has  a  laminar  sublayer  close  to  the  wall  and  a  smaller  velocity  gradient  du/dy  outside 
because  the  effective  (eddy)  viscosity  is  greater  than  the  molecular  viscosity  so  that  the 
shear  stress  x  =  |i  du/3y  remains  roughly  constant. 

In  physical  terms  we  can  understand  the  boundary  layer  as  that  region  of  the  fluid  in 
which  collisions  of  molecules  with  the  wall  are  important  Thus,  the  flow  of  a  fluid  past  a 
solid  surface  is  divided  into  two  regions.  Within  the  boundary  layer  we  must  consider 
collisions,  but  the  mean  flow  speed  tends  to  be  relatively  low  so  that  the  comparatively 
simple  equations  of  incompressible  flow  may  be  adequate,  while  outside  the  boundary 
layer  the  flow  speed  is  high  so  that  compressibility  effects  are  important,  but  the  effects  of 
collisions  do  not  have  to  be  considered  explicitly.  Thus  the  effects  on  the  flow  of  ablation 
will  only  be  important  in  the  boundary  layer  and  not  in  the  high-speed  flow  outside.  Some 
characteristics  of  ablation  are  discussed  in  the  Section  5.8  below. 

5.7  TRANSITION  FROM  LAMINAR  TO  TURBULENT  FLOW 

5.7.1  Introduction 

The  idealization  of  streamline  or  laminar  flow  is  observed  under  rather  specialized 
conditions,  but  once  the  characteristics  of  the  flow  itself  become  significant  (in  terms  of  lift, 
drag,  heat  transfer,  etc.)  the  conditions  are  normally  relatively  severe-e.g.,  high  fluid 
density,  high  speed,  high  shear,  or  whatever-so  that  the  smooth,  laminar  flow  generally 
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breaks  up  into  turbulent,  irregular  flow  with  greater  effective  frictional  drag  and  greater  heat 
transfer  than  in  the  idealized  case  of  laminar  flow. 

A  turbulent  flow  is  inherently  irregular  and  unstable  on  a  small  scale,  although  it 
has  some  statistical  regularities.  Thus,  the  theory  of  turbulent  flows  is  both  more 
complicated  and  also  less  well  developed  than  that  of  laminar  flows.  In  a  turbulent  flow  the 
energy  dissipation  rate  is  much  larger  than  the  molecular  viscous  dissipation  rate.  It  is 
often  convenient  to  introduce  an  effective  or  eddy  viscosity  coefficient  which  is  much  larger 
than  the  molecular  viscosity  |j.  which  is  introduced  in  Appendix  E,  Eq.  (E.18).  However, 
the  eddy  viscosity  is  not  just  a  property  of  the  medium  but  depends  on  the  history  of  the 
"air  parcel"  under  consideration.  One  cannot  simply  tabulate  and  look  up  an  eddy  viscosity 
in  general,  but  has  to  investigate  flows  of  the  same  medium  under  similar  conditions. 

To  ask  for  the  transition  from  laminar  to  turbulent  flow  one  has  to  specify  the 
conditions  of  the  flow,  and  here  I  give  three  distinct  and  different  examples.  The  numerical 
values  come  from  Schlichting  (1979),  pp.  39, 41, 26,  and  513.  In  particular,  conditions  in 
flows  involving  wails  (see  Sections  5.6.2,  5.6.3)  where  the  Reynolds  number  is 
significant  are  quite  different  from  free  atmospheric  flows  (see  5.6.4)  where  the  instability 
is  induced  by  the  interplay  of  buoyancy  and  wind  shear,  so  that  the  Richardson’s  number 

Ri  =  -  (g/p)  Op/az)/Ou/az)2wall  ,  (5.20) 

where  the  vertical  coordinate  z  increases  going  upward,  is  the  critical  parameter. 

5.7.2  "Poiseuille"  Flow  in  a  Tube 

Consider  the  flow  of  fluid  in  a  long  narrow  tube  of  diameter  d.  As  the  pressure 
head  and  hence  the  velocity  u  increases,  the  flow  becomes  turbulent  when 

Red  =  pud/n> Reel -2300-10,000  ,  (5.21) 

where  the  lower  limit  applies  to  a  tube  with  a  very  rough  walls  which  enhance  the  onset  of 
turbulence,  while  the  upper  limit  corresponds  to  exceedingly  uniform  conditions  in  the  flow 
(smooth  walls). 

5.7.3  Boundary  Layer  Flow 

Consider  flow  with  velocity  u  along  a  flat  plate.  A  boundary  layer  builds  up,  and  at 
a  distance  x  along  the  flow  for  which  the  Reynolds  Number  Rex  exceeds  a  critical  value, 
Re^'* 
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Rcx  =  p  u  x/p.  >  Rec2  ~  3.2  x  105  -  106  , 


(5.22a) 


the  flow  becomes  turbulent.  Note,  however,  that  this  can  be  expressed  differently.  In 
terms  of  the  boundary  layer  thickness  5  =  5(x)  of  Eq.(5.19),  the  flow  becomes  turbulent 
when 

Res  =  p  u  5  /p  >  Rec3  ~  Rcci  ~  2300  -  10,000.  (5.22b) 

5.7.4  Atmospheric  Turbulence 

In  free  atmospheric  flows  the  physics  is  different,  and  also  there  is  normally  no 
well-deflned  distance  such  as  d,  x,  or  5.  Thus  the  Reynolds  number  is  not  a  significant 
parameter.  Here  the  flow  becomes  turbulent  if  the  turbulent  energy  pr^uction  due  to  wind 
shears  is  large  compared  with  the  buoyancy  term.  In  this  case,  the  criteriuii  ior  turbulence 
is  expressed  in  terms  of  the  Richardson’s  Number  Ri  of  Eq.(5.20),  and  the  flow  becomes 
turbulent  if 

Ri<Ric  ,  (5.23) 

where  Ric  ~  0.042. 

5.8  SOME  CHARACTERISTICS  OF  ABLATORS  AND  ABLATION 

Ablation  is  a  self-regulating  mass  transfer  process  which  can  be  used  to  reduce  the 
frictional  heating  of  a  vehicle  during  entry  into  the  earth’s  atmosphere.  The  vehicle  surface 
is  coated  with  some  material  which  decomposes  by  a  variety  of  different  physical  and 
chemical  processes,  all  of  which  involve  taking  up  energy  from  the  gas  flow  and  thus 
reducing  heat  transfer  into  the  reentry  vehicle: 

(a)  Melting  and  vaporization  (two  phase  changes)-e.g.,  refractory  oxides  such  as 

quartz  or  silica 

(b)  Sublimation  (phase  change  directly  from  solid  to  vapor)"e.g..  Teflon 

(c)  Chemical  decomposition-here  a  variety  of  things  can  occur 

•  Depolymerization:  complex  (hydrocarbon)  polymer  breaks  down  into  a 
monomer  at  high  temperatures,  and  the  monomers  may  react  with  oxygen 
in  the  boundary  layer 

•  Pyrolysis:  thermal  decomposition  of  molecules  at  high  temperature.  For 
organic  compounds  this  leads  to  gases  (hydrogen,  methane,  acetylene,etc.) 
and  char  (carbon)  residue 

•  Combustion:  oxidation  of  material-e.g.,  graphite 
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•  Charting  ablators— such  as  carbon  phenolic,  quartz  phenolic,  silica  phenolic 
and  asbestos  phenolic-are  very  commonly  used. 

The  analysis  of  ablation  is  very  complex  and  highly  specific  to  the  panicular 
ablative  material  used.  It  is  necessary  to  analyze  the  changes  in  both  the  entry  vehicle  and 
also  in  the  boundary  layer,  considering  chemical  reactions  to  the  extent  necessary.  For  a 
historical  perspective  for  the  development  of  ablative  reentry  vehicles,  see  Sutton,  1982. 

5.9  THE  OPTICAL  SIGNATURE  OF  A  REENTRY  VEHICLE 

5.9.1  Introduction 

For  a  variety  of  applications  one  needs  to  know  the  range  of  values  of  the  optical 
(UV-LWIR)  signatures  of  generic  reentry  vehicles  (RVs).  There  exist  many  specific 
observations  and  calculations,  which  are  coordinated  through  the  USA/SDC-Teledyne 
Brown  Optical  Signatures  Code  (OSC),  but  before  one  goes  to  such  an  elaborate  model  it 
helps  to  have  a  simple  understanding  of  the  physical  factors  involved.  The  object  of  this 
discussion  is  to  provide  such  an  understanding. 

Explicitly,  as  representative  reentry  vehicle  we  consider  a  round-nosed  cone  of 
10-degree  half-angle  and  length  1.7  m,  which  is  described*^  in  Appendix  F.  Viewed 
broadside-on,  its  area  is 

Ab  ~  1.9  m2  .  (5.24) 

The  variability  in  actual  RV  signatures  is  due  to  the  following  factors:. 


(i) 

Size: 

the  RV  used  in  our  calculations  is  "average" 

(ii) 

Shape  and  weight: 

determines  peak  heating  altitude 

(iii) 

Ablative  coating: 

different  coatings  give  different  peak 
temperatures 

(iv) 

Reentry  angle  and  velocity: 

determines  peak  heating  altitude 

(V) 

Solar  illumination  geometry: 

this  is  very  important  between  UV  and  MWIR, 
where  solar  illumination  is  critical 

(Vi) 

Cloud  cover 

determines  the  reflection  of  sunlight  and  also 
has  an  influence  on  earthshine  emission. 

^  2  This  RV  has  the  same  ballistic  coefficient  P  as  the  RV  whose  temperature  profile  is  shown  in  Fig.  5.3. 
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As  a  result  of  factors  (i)-(iv),  the  LWIR  RV  signature  as  computed  here  is  likely  to 
vary  by  a  factor  2-3  at  high  and  low  altitudes,  but  by  a  factor  10  near  peak  heating  where 
the  factors  (ii)-(iv)  above  will  be  more  critical.  However,  in  the  shorter  wavelength  region, 
where  solar  illumination  predominates  (in  the  daytime),  variations  by  factors  of  100  or 
greater  are  possible.  There  may  be  effects  of  non-equilibrium  chemistry,  and  significant 
contributions  may  come  from  the  flow  field  of  the  reentry  wake. 

For  definiteness,  we  present  computations  for  three  different  wavelengths: 

•  X  =  0.5  |im,  visible 

•  X  =  3.5  ^m,  Mid-IR  (MWIR) 

•  X  =  12  |im,  Long-Wave  IR  (LWIR). 

These  are  all  atmospheric  transmission  regions,  so  that  atmospheric  transmittance  is 
generally  close  to  one.  Table  5.5  gives  the  numerical  values  used  here  for  a  variety  of 
parameters  at  these  different  wavelengths  The  parameters  are  defined  later  on  in  the  text, 
as  indicated. 

In  Section  5.9.2  we  discuss  the  shape  and  heating  of  the  RV;  in  Section  5.9.3  we 
compute  the  thermal  emission  signature  which  dominates  in  the  longer-wavelength  IR; 
Section  5.9.4  demonstrates  how  to  calculate  the  solar  scattering  (which  is  important  at  the 
shorter  wavelengths).  Section  5.9.5  gives  the  earthshine-scattering  contribution,  and 
Section  5.9.6  calculates  the  signature  due  to  sunlight  scattered  from  clouds  or  the  earth's 
surface.  Section  5.9.7  combines  all  these  results  and  Section  5.9.8  refers  to  the  Optical 
Signatures  Code  (OSC)  which  is  used  for  quantitative  reentry  vehicle  signatures. 

5.9.2  Typical  Shape  and  Heating  of  a  Generic  Reentry  Vehicle 

Figure  5.3  provides  the  computed  temperature  profile  for  a  particular  reentry  vehicle 
as  a  function  of  altitude  at  a  specific  range  and  reentry  angle.  In  free  space  [altitudes 
>  400,000  ft  (122  km)]  the  RV  is  at  a  uniform  temperature  of  approximately  300  K.  As  it 
descends  into  the  sensible  atmosphere,  a  bow  shock  begins  to  develop  so  that  first  the  nose 
heats  up:  thus  at  300,000  ft  (92  km)  there  is  considerably  less  radiation  than  is  observed  at 
lower  altitudes  (but  significantly  more  than  in  free  space).  For  the  particular  case  shown 
here  the  nose  temperature  has  reached  its  peak  value  (due  to  stagnation  point  shock 
heating-cf.  Section  5.4)  by  250,000  ft  (76  km).  As  the  RV  descends  further,  its  skirt 
temperature  increases  as  the  boundary  layer  builds  up.  Note  that  the  temperature  profile  at 
90,000  ft  (27  km)  altitude  shows  an  abrupt  rise  in  temperature  along  the  body.  This 
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indicates  the  transition  to  turbulence  (Section  5.7)  in  the  particular  computation  scheme 
used  by  AVCO  in  deriving  the  results  of  Fig.  5.3. 

For  computational  purposes,  we  break  up  the  target  (viewed  broadside-on)  into 
eight  segments  j  =  1,  2,  8,  each  of  fractional  area  aj  (i.e.,  Lj  aj  =  1).  The  details  of  the 
computation  are  given  in  Appendix  F. 


Table  5.5.  Parameters  Used  Here 


•  Transmittance: 

-  Because  these  are  all  transmission  bands  and  most  of  the  atmosphere  is  between  the 

earth  and  the  target,  we  shall  take  xs  = 

1  and  TE  =  0.8  at  all  wavelengths  considered. 

*  Emittance  and  Reflectance: 

X(pm) 

0.5  3.5 

12 

Target  emittance  e 

0.5  0.7 

0.9 

Target  reflectance  p(=  1  -e) 

0.5  0.3 

0.1 

Earth  emittance  eE* 

0.3  0.9 

0.9 

•  Earth  Surface  Temperature: 

From  LOWTRAN-5  Midlatitude  Summer  Atmosphere  we  take  Te  =  287  K 

*  From  McCtatchey  at  al.,  1 972-see,  e.g..  Part  I  of  these  Notes,  Fig.  3.3. 


5.9.3  Thermal  Emission  Signature 

If  we  divide  the  target  into  segments  aj,  then  for  a  reentry  vehicle  having  an 
emissivity-area  product  eA,  the  temperature  history  of  Fig.  5.3,  the  thermal  emission  in  a 
spectral  band  k  at  height  h  is 

Lr.e  (W/sr-pm)  =  eA  2^  aj  B[Xk;Tj(h)]  ,  (5.25) 

where  Tj(h)  =  temperature  of  segment)  at  height  h  (see  Fig.  5.3) 

Xlc  =  center  wavelength  of  band  k  ()im) 

B(X,T)  =  blackbody  spectral  radiance  (W/m2-sr-p.m)  for  wavelength  X  and 

temperature  T. 
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The  results  for  Lj.e  for  all  three  wavelengths  and  four  altitudes  are  shown  in 
Appendix  F;  the  results  will  be  set  in  context  and  discussed  in  Section  5.9.7. 

5.9.4  Diffuse  Scattering  of  Sunlight 

At  the  shoner  wavelengths,  a  target  is  seen  largely  by  reflected  sunlight,  and  here 
the  magnitude  of  the  signature  depends  on  the  solar  spectral  irradiance  Ms(X)  and  also-to  a 
critical  degree— on  the  character  of  the  scattering.  The  reflectance  of  both  sunlight  and 
earthshine  may  generally  be  taken  as  diffuse  Lambertian  characterized  by  a  bi-directional 
reflectance  so  that  Lsca>  the  radiance  of  the  surface  (W/m2-sr-|im)  is  related  to  Ejnc,  the 
irradiance  on  the  surface  (W/m^-jim)  by  the  relation: 

Lsca  (W/m2-sr-|im)  =  (p/Jt)  Fine  (W/m^-pm)  .  (5.26) 

The  computation  of  the  signal  due  to  sunlight  scattered  diffusely  by  the  target  is 
sketched  in  Fig.  5.6.  Explicitly,  the  solar  (spectral)  irradiance  Mg  (W/m^-pm)  is  attenuated 
by  the  transmission  factor  scattered  diffusely  (see  Eq.  5.26)  by  the  target  of  effective 
area  A  to  give  the  signal 

Lt,s  (W/sr-pm)  =  A  Mg  xg  cos())i  (p/n)  cos(>2  ,  (5.27) 

where  the  target  is  modeled  as  a  flat  plate  with  A  =  Ab  of  Eq.  (5.24)  and  the  angles  are 
defined  in  Fig.  5.6.'^  This  is  then  attenuated  by  the  transmittance  XpaUi  from  target  to 
sensor  to  give  a  target  radiance  Lx ^  Xpaih- 

Table  5.6  gives  the  solar  scattering  at  X  =  0.5,  3.5,  and  12  pm  for  a  specific  set  of 
numerical  assumptions  that  are  called  out  in  the  text  and  foomotes.  Note  that-unlike  the 
thermal  emission  signal  discussed  in  Section  5.9.3-the  results  are  independent  of  altitude  h 
or  RV  temperature. 


^  ^  Note  that  there  is  a  large  contribution  due  to  sunlight  scattered  from  the  earth  or  clouds-see  Section 
5.9.6. 

^  p  =  1  -  e,  where  e  =  emissivity  of  the  surface.  Numerical  values  used  are  listed  in  Table  5.4. 

^  ^  All  transmittance  factors  are  taken  equal  to  one  in  the  numerical  examples  because  we  are  dealing  with 
transmission  windows  in  the  atmosphere. 

1 6  >ve  use  =  ^  =  45°,  i.e.,  cos  4>i  x  cos  <>2  =  0.5. 

*  ^  Except  for  possible  differences  in  scattering  geometry  and  atmospheric  transmittance  with  altitude. 
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EXITANCf 


Figure  5.6.  Computing  the  Target  Irradlance  Due  to  Solar  Scattering 
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Table  5.6.  Representative  Numbers  for  Solar  Scattering 


A.  (pm) 

wm 

3.5 

■■ 

Solar  irradiance  Hs(X,)  (W/m?-pm) 

940 

14.6 

0.12 

Bidirectional  reflectance  p 

0.5 

0.3 

0.1 

Scattered  Signal,  Lj.s  (W/sr-pm) 

280 

2.6 

7.1  E-3 

5.9.5  Earthshine  Scattering 

The  computation  of  the  signal  due  to  earthshine  scattered  diffusely  by  the  target, 
Lt^  ,  is  sketched  in  Fig.  5.7.  This  computation  is  somewhat  more  complex  than  that  of 
Lt^S  because  one  has  to  integrate  contributions  from  different  portions  of  the  earth's 
surface,  in  contrast  to  the  sun  which  can  be  treated  as  a  (distant)  point  source. 

The  signal  from  the  target  due  to  earthshine  scattering,  Ljje,  is  obtained  by  starting 
with  earth  radiance  Le  =  Ee  B(Te),  tracking  the  irradiance  of  the  target,  and  then  going 
through  the  diffuse  scattering  process  of  Eq.  (5.26).  Figure  5.7  sketches  the  procedure. 
Given  Le,  one  attenuates  by  the  atmospheric  transmittance  factor  to  the  target 

‘CE(<l>a)  =  XE^ota  ,  (5.28) 

multiplies  by  the  solid  angle  factor  dAE/[RE(0a)]^  reduced  by  the  projection  cos<t>i,  and  then 
integrates  over  the  whole  accessible  surface  of  the  earth.  This  somewhat  intricate 
procedure  is  discussed  in  detail  in  Appendix  G.  Next  we  have  diffuse  (Lambenian) 
scattering  as  in  Eqs.  (5.26)  and  (5,27),  which  gives  the  following  expression  for  the 
earthshine-scatteted  signal  of  the  target  as  seen  by  the  sensor,  Ljje, 

=  Ab  JdAE  Le  COS<|>i  [RE(<l>a)]"^  (P/’l)  cos<})2 

=  Ab  2  eE(X)  B(X,Te)  P  cos(|>2  M'^e.P) 


(5.29) 

(5.30) 


MflUNCf  AT  TAMH  ME  TO  EAlTNtlMI 


Figure  5.7.  Computing  the  larget  Irradiance  due  to  Earthshine  Scattering 


where  the  integral  IjC'tE.P)  is  discussed  in  Appendix  G,  and  the  angle  of  attack,  P,  as  well 
as  the  angles  (|)i  and  ^2  defined  in  Fig.  5.7.*®  Numbers  analogous  to  Table  5.6  are 
listed  in  Table  5.7.  Numerical  values  for  the  various  parameters  used  are  listed  in 
Table  5.5. 


Table  5.7.  Earthshine  Scattered  Target  Signal 


X(pm) 

m 

3.5  12 

0.3 

0.9 

B 

Earth  Radiance  Term,  B(A.,Te)  (W/m^-sr-tim) 

1.1  E-34 

0.14 

Earthshine  radiance,  Lj.e  (W/sr-pm) 

3.0  E-36 

7.0  E-3 

IB 

*  Note  that  tE  =  0.8,  i.e.,  b  =  In  (1/te)  =  0.22.  and  thus  frlb)  * 
Appendix  G. 

7  X  10“2  from  Fig.  G.2  of 

5.9.6  Signal  due  to  Sunlight  Scattered  from  the  Earth  or  Clouds*^ 

At  the  shorter  wavelengths,  an  important  contribution  to  the  target  signature  comes 
from  sunlight  that  is  scattered  from  the  earth  or  from  cloud  tops  and  is  then  reflected 
diffusely  from  the  target.  Reference  to  Fig.  3-3  of  Part  I  shows  that  the  reflectance  of 
various  natural  surfaces  in  the  visible  and  near-IR,  Pnat>  is  variable  but  quite  significant. 
Thus  the  solar  irradiance  Ms  (W/m^-pm)  is  attenuated  by  a  transmission  factor  t’s  and  the 
signal  signal  reflected  from  the  earth  or  from  cloud  tops  is  Ms  x's  (pnai/^).  so  that  from  an 
analysis  similar  to  that  for  earthshine  scattering  in  Section  5.9.5,  [see  Eq.(5.30)],  one  finds 
that 

Lt.SS  =  Ab  Ms  T’s  (Pnat/Jt)  P  COS())2  lT(b)  .  (5.31) 

Table  5.8  gives  values  of  Ms,  Pnai.  P,  and  Lx,ss  for  X  =  0.5 , 3.5  and  12  pm;  in 
the  calculation  we  have  used  the  same  numerical  values  for  the  various  parameters  as  in 
Sections  5.9.4  and  5.9.5,  namely  cos  <j)2  =  0.707,  t’s  =  Xe  =  0.8,  lyCb)  =  7  x  10“2. 


*®  We  use  P  =  01  =  <t>2  =  45  deg. 

*  9  I  wish  to  thank  Dr.  Rich  Loda  for  pointing  out  the  impotance  of  this  term. 
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Table  5.8  Signal  due  to  Sunlight  Scattered  from  the  Earth  or  Clouds 


A.  (pm) 

■a 

3.5 

12 

Solar  Irradiance  Hs(A.)  (W/m^-pm) 

940 

14.6 

0.12 

Earth  (or  cloud)  reflectance,  pnat 

0.5 

0.2 

0.1 

Target  bidirectional  reflectance,  p 

0.5 

0.3 

0.1 

Radiance  due  to  scattered  sunlight.  Lj.sslW/sr-pm) 

240 

2.1E-2 

2.8E-5 

5.9.7  Total  Signature  of  a  Reentry  Vehicle 

With  the  exception  of  specular  reflection  or  glints, the  signature  of  a  reentry 
vehicle  is  the  sum  of  contributions  due  to  solar  scattering,  Lt,s.  thermal  emission  Lj,©, 
and  earthshine  scattering  Lt£  (which  turns  out  to  be  small).  Combining  the  results  of 
Tables  5.6  through  5.8,  we  show  the  results  in  Table  5.9,  which  are  quoted  in  units  of 
W/sr-pm: 

(a)  In  Table  5.9  the  only  figures  which  vary  with  altitude  (i.e.,  temperature-- see 
Fig.  5.3)  are  those  for  the  thermal  emission,  Lj,©,  which  are  largest  at  peak 
heating  near  75,000  ft  (23  km)  altitude. 

(b)  Lt.©  varies  strongly  with  temperature:  note  that  the  visible  and  MWIR  radiance 
is  very  large  because  at  a  peak  temperature  near  2600  K  the  RV  glows  brightly. 

(c)  In  fact,  Lt.©  is  much  larger  than  the  solar  scattering  contribution,  Lt.s.  in  the 
visible  near  peak  heating. 

(d)  Recall  that  non-equilibrium  considerations  and  wake  contributions  are  not 
included  here. 

5.9.8  Introduction  to  the  Optical  Signatures  Code  (OSC) 

The  OSC  is  a  computer  program-now  (1989)  in  its  15th  version-which  has  been 
developed  during  the  last  15  years.  Its  objective  is  to  develop  IR  optical  signatures  of 
reentry  vehicles,  PenAids,  chaff,  etc.,  in  Midcourse  and  Reentry  phases,  in  a  computerized 


The  signal  due  to  glint  is  Lj.s*  =  Ab  r  Ms  Ts  cos^i  x  cos4>2  instead  of  Eq.  (5.27),  where  r  =  specular 
reflectivity,  which  may  be  close  to  one.  Reference  to  Table  S.6  shows  that  this  can  be  significantly 
larger  than  Lj^s  .  However,  the  actual  reflected  power  has  to  be  multiplied  by  the  fraction  of  the  toud 
area  that  gives  such  specular  reflection,  which  is  normally  very  much  smaller  than  one.  Accordingly, 
this  contribution  to  the  signal  is  not  considered  here. 
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Table  5.9.  Total  Signature  of  Reentry  Vehicle 


X(4m) 

0.5 

3.5 

12 

ti=.aQ0Kft(.9l  hm). 

Solar  scattering, 

l-T,S 

280 

2.6 

7E-3 

Thermal  emission, 

l-T,e 

2  E-8 

38 

45 

Earthshine  scattering. 

Lt,e 

3  E-36 

7E-3 

0.1 

Solar-Earth/Cloud  Scatt. 

Lt,SS 

240 

1E-2 

2.8E-5 

Total  signal 

Ltot 

520 

41 

45 

h  =  150l<ftf46  km) 

•-T,S 

280 

2.6 

7E-3 

•-T,e 

840 

6500 

340 

Lt.e 

<  — 

—  smaD  — 

—  > 

Lt,SS 

240 

small 

small 

Llot 

1300 

6500 

340 

h  =  75  kft  (23  kml  Peak  Heating. 

Lt,s 

280 

2.6 

7E-3 

LT,e 

53,000 

76,000 

1400 

Lt,e 

<  — 

—  smal  — 

—  > 

Lt,SS 

240 

small 

small 

Ltot 

29,000 

41,000 

730 

h  =  25Kn(7.6km) 

LT,S 

280 

2.6 

7E-3 

Lt.© 

5E-9 

840 

180 

l-T,E 

<  — 

—  smal  — 

—  > 

Lt,SS 

240 

small 

small 

Ltot 

520 

840 

180 

form  which  is  accessible  to  the  user.  Because  of  long-term  continuity  by  the  developer 
(Teledyne  Brown  Engineering,  HuntsviUe,  AL  35807,  Tel.  205-726-xxxx),  there  is  a  lot  of 
collective  memory  (by  John  Davis  x  1125,  S.  Murty  x  1109,  et  al.);  but  because  of  the 
computer  format  it  is  frequently  easier  to  get  numbers  than  human  judgment.  In  the 
production  of  the  code  much  laboratory  and  field  data  have  been  absorbed  and  analyzed. 
There  is  in  addition  a  great  deal  of  information  on  material  properties  (optical,  thermal, 
etc.),  on  reentry  flow  fields,  heat  transfer,  ablation,  and  on  optical  signatures.  Murty 
(1987)  has  presented  some  very  nice  lecture  notes  on  endoatmospheric^^  heating  theory. 

Information  about  the  code  may  be  obtained^  from  Ellen  Montenegro,  TBE-1413, 
or  from  Tom  Prestwood,  USASDC,  Passive  Sensor  Division,  DASD-H-YP,  Huntsville, 
AL  35807,  Tel.  205-895-4572.  While  the  code  is  largely  Unclassified,  the  validation  data 
make  the  whole  thing  Secret.  The  code  as  a  whole  consists  of  a  number  of  different 
codes.^  To  get  results  that  can  be  trusted,  significant  effort  on  the  user's  part  is  required. 


"Endoatmospheric”  means  within  the  atmosphere,  i.e.,  at  altitudes  below  60-100  km,  as  contrasted 
with  "exoatmospheric,"  or  above  the  sensible  atmosphere. 

These  names  and  phone  numbers  are  correct  as  of  October  1989. 

23  Many  of  the  individual  codes  were  developed  by  various  reentry  contractors  during  the  1960s,  so  that 
there  are  {xobably  deficiencies  in  corporate  memory. 
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6.0  PRIMER  ON  NUCLEAR  WEAPONS  EFFECTS 


6.1  INTRODUCTION 

When  a  nuclear  weapon  is  detonated,  a  very  large  amount  of  energy  is  produced, 
either  by  the  fission  of  uranium  or  plutonium  or  by  fusion  reactions  involving  D,  T,  or  ^Li. 
All  these  reactions  produce  neutrons  and  sustain  chain  reactions;  fission  produces  on  the 
order  of  200  MeV  per  atom  while  a  fusion  reaction  produces  3-17  MeV,  and  the  yield-to- 
mass  ratio  obtained  by  either  mechanism  is  the  same  within  a  factor  3  or  so.  One  kiloton 
(kt)  of  energy24  is  10^2  calories  or  4.2  x  lO^^  Joules,  and  corresponds  to  the  complete 
fission  of  0.057  kg  of  uranium.  The  detonation  is  by  no  means  100  percent  efficient;  a 
representative  mass  of  a  nuclear  weapon  is  of  the  order  of  lO^-lO^  kg.  Thus  we  see  that 
most  of  this  material  remains  as  Bomb  Debris.  This  material  is  heated  to  a  very  high 
temperature  during  the  detonation,  being  largely  vaporized  and  ionized.  A  significant 
fraction  of  the  energy  of  the  detonation  appears  as  debris  kinetic  energy.  The  debris 
expansion  velocity  u  is  (see,  e.g.,  Holland,  1977,  p.  162.) 

u -(Y/M)  1/2 10^-10*  cm/sec  ,  (6.1) 

where  Y  (Megatons,  Mt)  is  the  total  yield  of  the  weapon  and  M  is  the  total  mass 
participating  in  the  expansion. 

Most  of  this  energy  is  emitted  in  a  time  of  the  order  of  0.1-1  p.sec.  In  this  time  the 
bomb  radiates  essentially  all  of  its  energy  as  a  blackbody  of  1-10  keV  blackbody 
temperature,  which  corresponds  to  EM  radiation  of  wavelength  12-1.2  A,  i.e.,  soft  X-rays. 
The  total  energy  output  of  a  nuclear  weapon  is  distributed  roughly  as  is  shown  in  Table 
6.1.25  Note  that  the  largest  part  of  the  energy  output  is  in  the  form  of  (3-30  KeV)  X-rays, 
while  debris  kinetic  energy  provides  most  of  the  resL^^ 


Which  corresponds  roughly  to  the  heat  of  combustion  of  1000  tons  of  TNT. 

25  See,  e.g.,  Holland,  1977,  p.  162  (but  more  X-ray  yield  than  listed  there),  also  ENW,  1964  or  1977, 
Ch.1-2. 

25  For  orientation  on  mass-yield  relations,  note  that  (a)  If  Y  =  20  kT  =  8  x  10*^  J,  a  velocity  of 
102-10®  cm/sec  means  that  a  mass  of  order  100-1000  kg  is  involved;  (b)  The  U.S.  8-in.  (203-mm) 
Howitzer  Ml  15  fires  projectiles  in  the  90-95  kg  range.  One  of  these  is  the  nuclear  M422  round  whose 
yield  is  in  the  5-10  kt  range  [Jane's  Armour  and  Artillery  (1986-87)  p.  673.] 
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Table  6.1  Typical  Energy  Distribution  from  a  Nuclear  Explosion  In  Vacuum 


X-ray  yield 

70% 

Debris  kinetic  energy 

25% 

Delayed  nuclear  radiation 

5% 

Prompt  neutrons 

0.1-1% 

Prompt  gamma  rays 

0.01-0.1% 

It  is  clear  that  only  a  small  fraction  of  the  mass  of  the  device  undergoes  fission,  and 
only  a  small  fraction  of  the  yield  comes  out  as  nuclear  radiation.^^  The  nuclear  radiation  is 
frequently  not  considered  in  the  analysis  of  vulnerability,  and  in  most  discussions  of 
phenomenology  it  does  not  matter  whether  the  weapon  is  a  fission  or  fusion  device,  since 
most  of  the  energy  comes  out  as  X-rays  and  debris  kinetic  energy. 

If  the  energy  is  assumed  to  be  distributed  isotropically,  then  the  X-ray  flux  W(R) 
due  to  a  warhead  of  yield  Y(Mt)  at  range  R  (km)  is^ 

W(R)  =  6400  (Y/R2)  caVcm2  ,  (6.2) 

and  this  is  the  principal  kill  mechanism  in  the  high-altimde  range.  Thus,  for  a  vehicle 
hardened  to  1  cal/cm^,  the  kill  radius  for  a  1  Mt  bomb  is  270  km. 

A  typical  total  exposure  due  to  a  nuclear  weapon  detonated  in  vacuum  is  listed  in 
Table  6.2.  The  shielding  listed  for  the  y-ray  dose  is  not  large,  but  is  representative  for  a 
space-based  sensor. 

Within  the  atmosphere,  the  bomb-generated  X-rays  are  absorbed  rather  rapidly. 
Table  6.3  indicates  the  atmospheric  extinction  of  different  bomb  radiations  in  the 
atmosphere.  We  see  that  the  (soft)  X-rays  are  absorbed  below  ~  70  km,  while  neutrons 
and  gamma  rays  are  transmitted  down  to  ~  20  km.  Note: 


2^  Electrons  and  protons  and  other  positive  ions  interact  so  strongly  with  ambient  atoms  and  molecules 
that  they  do  not  travel  far  enough  to  be  significam.  but  the  uncharged  neutrons  and  gamma  rays  interact 
sufficiently  weakly  with  the  ambient  environment  that  they  travel  a  long  way  and  thus  can  be 
significant  from  the  standpoint  of  their  effect  on  biological  organisms  and  on  electronic  solid  state 
devices,  which  are  very  sensitive  to  ionizing  radiation. 

28  This  is  the  standard  equation,  which  assumes  that  80%  of  the  bomb's  yield  goes  into  X-rays  distributed 
isotropically.  (The  figure  of  80%  is  conventional  here,  as  compared  with  70%  or  other  values  in 
Table  6.1.) 
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(a)  The  absorption  and  extinction  (=  absorption  +  scattering)  coefficient  of  X-rays 
in  air  as  a  function  of  energy  is  given  in  Fig.  4-3  of  Part  I  of  these  notes 
(Bauer,  1990).  Note  that  the  extinction  coefficient  falls  off  very  rapidly  with 
increasing  energy,  and  scattering  is  only  important  at  the  highest  energies . 

(b)  Typically,  the  X-ray  spectrum  of  a  nuclear  weapon  corresponds  to  a  blackbody 
(Planck)  distribution  of  approximately  1  keV  temperature.  The  Planck 
spectrum  is  of  the  form  29 

F(u)  =  15  7C-4u3/(eu- 1)  ,  (6.3) 

where  u  =  E/kT;  note  that  F(u)  has  a  maximum  at  u  =  2.8  . 


Table  6.2.  Total  Exposure  In  Vacuum  from  a  Nuclear  Weapon  of 
Yield  Y(Mt)  at  a  Range  R(km)* 


X-ray  flux 

6400  (Y/r2) 

cal/cm2 

y-ray  dose 

2.1  {Y/R2) 

Mrad(Si)  behind  23  mil  Ta  and  2  g/cm2  Ar* 

y-ray  dose  rate 

1.9  x  10^3  (Y/R2) 

rad(Si)/sec 

Neutron  flux 

1.6  X  10^5  (Y/R2) 

cm~2 

*  Numbers  kindly  provided  by  William  Chadsey  and  Rich  Sutton  of  SAIC. 

*  *  1  rad(Si)  corresponds  to  the  absorption  of  1 00  erg/gm  in  Si.  If  we  assume  that  all  the 
■y-rays  of  effective  energy  1  MeV  -  1.6  x  10^  ergs  are  absorbed  in  1  g/cm^,  then  1 
Mrad  corresponds  to  ay-ray  flux  of  10*/1 .6x10^ cm“2.  6.25x1 0''3cm~2  23milTa 
and  2  g/cm^  Al  is  a  minimum  level  of  shielding  normally  used  in  spacecraft. 

We  see  from  Table  6.3  that  at  altitudes  below  50  km  the  X-ray  energy  (which 
represents  most  of  the  bomb  output)  is  absorbed  near  the  burst,  giving  rise  to  a  well- 
confined  firebalpo  whose  mass  is  actually  of  order  of  10^  tons/Mt  yield.  Table  6.4 
indicates  schematically  how  the  fireball  is  created  and  enumerates  the  damage  mechanisms 
within  the  atmosphere,  which  are  sketched  out  in  Section  6.2.4,  following  a  discussion  in 
Section  6.2.3  of  near-surface  bursts. 
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See,  e.g.,  Eq.  (3.1),  also  Figs.  3.1  and  3.2,  all  in  Section  3  of  Part  I  of  these  Notes  (Bauw,  1989). 

30  In  fact,  the  debris  is  also  stopped  near  the  burst  point  in  the  lower  aunosphere;  at  high  altitudes  the 
(largely  ionized)  debris  follows  the  geomagnetic  field  lines. 
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Table  6.3.  Atmospheric  Transmission  for  a  Horizontal  Path  Thickness 
of  10-km  Length  for  3  keV  X*rays,  1  MeV  Neutrons, 
and  1  MeV  Gamma  Rays  at  Different  Altitudes. 


Altitude 

(km) 

Dens.(kg/m3) 

for; 

Transmission  through  10  km  horizontal  path 

3  keV  X-rays  1  MeV  neutrons  1  MeV  gamma  rays 

(Ext.  coefft.): 

(170  cm2/g) 

(1.7  X 10-2^  cm2/atom) 

(0.1  cm2/g) 

0 

1.23 

0 

0 

0 

10 

0.41 

0 

2  X  10-‘'8 

4  X  lO-"'"' 

20 

0.089 

0 

1  X  1 0-"^ 

5.5  X  10-3 

30 

0 

0.17 

0.72 

40 

0.0040 

0 

0.67 

0.79 

50 

0 

0.90 

0.94 

60 

1  X 10-23 

0.97 

0.98 

70 

8.8  X  10"5 

3  X  10-7 

0.99 

1.00 

80 

0.03 

1.00 

1.00 

90 

0.58 

1.00 

1.00 

100 

5.0  X  10-7 

0.92 

1.00 

1.00 

Table  6.4.  Effects  of  a  Nuclear  Weapon  Exploded  in  the  Lower  Atmosphere 


X-rays  absorbed  in  Atmosphere 

Shock/Btast  wave  ->  Damage 

Rreball 

Debris  slowed  down  in  Atmosphere 

Thermal  Radiation  -♦  Damage 

Nuclear  Radiation  (neutrons  and  y-rays*) 

- )  Damage  (Secondary) 

*  Only  weakly  absorbed  in  atmosphere  -  see  Table  6.3. 


At  altitudes  above  75-100  km  the  ambient  density  is  so  low  that  thermal  X-rays  are 
no  longer  absorbed  within  the  atmosphere,  and  so  the  effects  of  a  nuclear  explosion  are 
qualitatively  different  than  at  lower  altitudes.  Table  6.5  sketches  the  phenomenology  and 
damage  mechanisms  in  this  altitude  regime. 
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Table  6.5.  Effects  of  a  Nuclear  Weapon  Exploded  In  the  Upper  Atmosphere 


X-rays* 

*  Emitted  upwards,  propagate  freely 

-» 

Damage  (Primary  Kill) 

•  Emitted  downwards,  absorbed  to  heat  air 

Heated  Air**  •  Ionized 

-» 

Radio  Blackout 

•  Radiates 

-> 

Redout 

Debris*** 

Transfers  energy  to  heat  air 

Nuclear  Radiation  (propagates  freely) 

Damage  (Secondary) 

*  Distributed  roughly  isotropicaliy. 

**  May  be  moved  upwards  by  Heave  (see  Section  6.3.3) 

***  Travels  along  (possibly  perturbed)  geomagnetic  field  lines. 


Overall  high-altitude  phenomenology  is  discussed  in  Section  6.3.  Section  6.3.1 
describes  single-burst  phenomenology  at  high  altitudes;  unlike  the  limited  extent  of  nuclear 
fireballs  within  the  atmosphere,  the  very  tenuous  upper  atmosphere  is  disturbed  out  to 
distances  of  hundreds  of  kilometers,  and  thus  not  only  the  long-range  effects  of  individual 
bursts  but  also  the  interaction  between  multiple  bursts  become  very  important.  Multiburst 
phenomenology  is  discussed  in  Section  6.3.2  and  a  particular  phenomenon  called  Heave, 
which  involves  the  rise  of  relatively  large  masses  of  air  heated  by  the  absorption  of  X-rays 
and  UV  radiation,  is  treated  in  Section  6.3.3. 

At  high  altitudes  there  is  a  highly  ionized  region  of  large  extent  which  interferes 
with  radio  wave  propagation;  this  "blackout"  is  discussed  in  Section  6.4.2;  there  is  also  a 
large  region  of  heated  air  which  radiates  (mainly  in  the  in&ared)  to  produce  an  enhanced  IR 
background,  the  so-called  "redout"  which  is  discussed  in  Section  6.4.3. 

It  must  be  stressed  that  the  high-altitude  multi-burst  phenomenology  is  based  on 
models  which  are  not  based  on  actual  nuclear  test  observations.  Section  6.5  discusses 
overall  uncertainties  in  nuclear  effeas  predictions. 
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NUCLEAR  EXPLOSIONS  IN  THE  LOWER  ATMOSPHERE 


6.2.1  Introduction 

The  previous  discussion—see  in  particular  Tables  6.3  and  6.4— points  out  that  a 
fireball  is  produced  at  sufficiently  high  ambient  densities  or  sufficiently  low  altitudes.  The 
phenomenology  of  fireball  formation  and  evolution  in  the  atmosphere  is  described 
qualitatively  in  Section  6.2.2,  which  sketches  the  evolution  and  cooling  of  the  fireball  as  it 
goes  from  the  initial  X-ray  temperature  of  1-10  keV  (10^-10*  K)  to  ~  1  eV  (12,000  K) 
within  some  10  sec  for  a  Mt-range  weapon  (for  more  detail,  see  Erode,  1968).  Section 
6.2.3  gives  a  simple  numerical  model  of  the  evolution  of  the  fireball  and  its  blast  and 
thermal  damage  mechanisms.  Two  special  issues  for  a  near-surface  burst— crater  formation 
and  atmospheric  dust  loading— are  discussed  briefly  in  Section  6.2.4.  Some  numbers  for 
damage  mechanisms  due  to  a  low  air  burst  are  sketched  in  Section  6.2.5  and  Electro¬ 
magnetic  Pulse  (EMP)  effects  are  outlined  in  Section  6.2.6. 

Let  us  begin  with  some  numbers  for  orientation.  As  a  measure  of  the  size  of  a 
nuclear  fireball,  its  effective  radius  at  thermal  minimum^^  for  a  near-surface  burst  of  yield 
Y  (Mt)  is 

Reff-RaY04  ,  (6.5) 

where  Ra  ~  430  m,  i.e.,  Reff  -  90  m  for  a  20  kt  weapon  or  430  m  for  a  1  Mt  weapon  (see 
ENW,  1977,  p.  70). 

Figure  6. 1  shows  how  a  low-altitude  fireball  expands  with  increasing  time,  and 
Fig.  6.2  shows  the  stabilized  cloud  rise  height  of  the  nuclear  cloud  as  a  function  of  weapon 
yield  at  mid-latitudes.  To  get  a  feel  for  the  order  of  magnitude  of  the  energy  of  a  nuclear 
detonation,  note  that  a  20-100  kt  weapon  rises  to  the  mid-latitude  tropopause  (-11  km). 
This  can  be  understood,  since  reasonably  large  thunderstorms— which  have  a  comparable 
energy  content— rise  to  this  same  height.^2 


■m 


^  ^  This  will  be  defined  in  Section  6.2.2. 

The  energy  that  drives  a  thunderstorm  is  the  latent  heat  of  vaporization  of  the  atmospheric  moisture, 

which  is  in  the  range  600-700  cal/g  which  is  (2.5-3)  x  10^  J/kg  or  0.6-0.7  kt(energy)/kt(mass).  The 

water  vapor  density  in  a  cloud  is  1-10  g  H20/m3  or  10^-10^  tons  H2O  per  km^,  so  that  the  latent  heat  • 

of  a  thunderstorm  of  volume  1  km^  is  -  0.6-7  kt  of  energy. 
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TIME  AFTER  EXPLOSION  (SECONDS) 


Figure  6.1.  Fireball  Expansion  (ENW,  p.  74)  "Breakaway"  corresponds  to  the 
situation  when  the  hot,  incandescent  Interior  of  the  fireball  can  be  seen 
through  the  air  shock,  which  cools  and  becomes  less  opaque 
as  It  expands  •  see  discussion  in  Section  6.2.2. 


Figure  6.2.  Nuclear  Cloud  Rise  Height  as  Function  of  Yield. 
[Synthesized  by  H.  Mitchell  (Falcon  Research)  from 
various  sources,  see,  e.g.,  Bauer  (1979)] 
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6.2.2  Qualitative  Description  of  a  Low  Air  Burst  [quoted  from  J.  Zinn 
(LANL--see  Barasch,  1979)]. 

Within  the  first  microsecond  after  detonation,  the  entire  nuclear  yield  is 
deposited  in  the  bomb  materials,  heating  them  to  extremely  high 
temperatures  in  the  range  107  k.  Some  of  this  energy  is  immediately 
radiated  away  in  the  form  of  X-rays  and  extreme  ultraviolet  radiation,  but 
this  radiation  is  immediately  reabsorbed  in  the  air  within  the  first  few  meters 
surrounding  the  device.  The  air  is  thereby  heated  to  a  very  high  tempera¬ 
ture,  in  the  range  of  106  K,  with  a  corresponding  increase  in  pressure  to 
several  thousand  atmospheres.  This  heated  air  and  the  hot  bomb  vapors 
(debris)  in  the  interior  constitute  the  initial  fireball.  The  fireball 
subsequently  expands  through  a  combination  of  radiative  and  hydrodynamic 
processes.  An  intense  shock  wave  (thin  region  of  highly  compressed  air)  is 
formed  at  the  surface,  which  expands  outward  at  a  high  velocity.  Inside  the 
shock,  the  energy  is  rapidly  redistributed  through  emission  and  reabsorption 
of  ultraviolet  radiation. 

Fig.  6.3  shows  the  thermal  power  radiated  from  a  nuclear  burst  at  different 
altitudes  as  a  function  of  time.  At  the  low  altitudes  under  discussion  here, 
the  power  is  emitted  in  two  pulses. 

The  time  variation  of  light  emission  during  the  first  pulse  of  Fig.  6.3  is 
controlled  by  the  radial  growth  of  the  shock.  Although  the  shock  surface  is 
intensely  luminous,  the  shock  is  opaque  to  visible  light,  thus  concealing  the 
air  and  bomb  debris  in  its  interior.  As  the  shock  expands,  it  engulfs  cold 
air,  causing  its  temi^rature  to  decrease,  with  a  corresponding  decrease  in 
brightness.  The  optical  energy  emission  rate  (power)  is  pro|wrtional  to  the 
fireball  brightness  times  the  area  of  its  surface.  After  detonation,  the  optical 
power  first  increases,  because  the  increase  in  surface  area  outweighs  the 
decrease  in  brightness.  Later  it  decreases,  because  the  decrease  in 
brightness  is  dominant. 

During  the  early  expansion  and  cooling  of  the  fireball,  the  ultraviolet 
radiation  in  the  interior  decreases  rapidly  until  it  cannot  effectively 
redistribute  the  energy.  Thereafter  the  shock  temperature  falls  more  rapidly 
than  the  temperature  in  the  interior.  As  the  expansion  proceeds,  the  shock 
continues  to  cool,  and  as  it  does  so,  it  becomes  less  and  less  opaque  to 
visible  light,  while  also  becoming  less  luminous.  A  point  is  reached, 
corresponding  to  the  minimum  temperature  region  in  Fig.  6.3,  where  the 
shock  becomes  sufficiently  transparent  to  allow  light  from  the  hotter  interior 
to  begin  to  escape,  causing  the  optical  output  to  begin  increasing  again33. 
Further  growth  of  the  shock  results  in  increasing  transparency  and  a  further 
increase  of  the  luminosity  of  the  inner  region. 

As  the  expansion  continues  further,  the  temperature  of  the  inner  fireball 
decreases,  due  to  the  combined  effects  of  hydi^ynamic  expansion  and  loss 
of  ener^  by  radiation.  This  cooling  eventually  results  in  the  decrease  in 
luminosity  following  the  second  maximum  in  Fig.  6.3.  After  the  second 


3  3  This  is  denoted  by  "breakaway"  in  Fig.  6.1 . 
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maximum,  the  inner  fireball  gradually  becomes  transparent,  revealing  the 
bomb  debhs  inside,  whose  brightness  also  decreases  with  time. 

The  majority  of  the  total  energy  radiated  by  the  fireball  comes  from  the 
second  peak,  which  is  not  very  different  in  its  instantaneous  brightness 
relative  to  the  first  peak,  but  last  about  100  times  longer.  By  the  time  this 
"thermal  pulse"  is  over,  the  available  energy  has  either  been  radiated  away, 
trapped  in  the  debris  and  entrained  air  of  the  fireball,  or  it  is  propagating  in 
the  shock  wave,  which  is  now  well  outside  the  visible  fireball. 

Note  that  with  increasing  height  of  burst,  the  minimum  in  radiated  power 
shown  in  Fig.  6.3  declines,  disappearing  above  40  km. 

Figure  6.4  gives  an  alternative  representation  of  what  is  shown  in  Fig.  6.3  by  listing  the 
evolution  with  time  of  the  fireball  temperature  from  a  20  kt  burst. We  note  that  the 
fireball  remains  in  the  1000-10,000  K  temperature  range  for  times  of  3-30  sec,  i.e.,  long 
enough  so  that  its  phenomenology  has  considerable  implications  to,  for  example,  strategic 
defense  systems.  For  a  discussion,  see  ENW,  1964,  Ch.  1  and  2,  in  particular  pp.  87- 
100.35 


3^^  Comparable  results  are  obtained  for  other  yields.  Minimum  and  maximum  temperabires  are  essentially 
independent  of  yield,  but  the  time  at  which  these  temperatures  occur  scale  approximately  with  the 
0.4th  power  of  yield-see  ENW,  1977,  p.  70. 

35  Sec  also  Zeldovich  and  Raizer,  1%6,  Vol.  II,  p,  61  If,  for  a  description  taken  from  the  first  (1957) 
edition  of  ENW,  which  Zeldovich  and  Raizer  describe  as  better  than  that  in  the  1962/64  edition.  Note 
that  there  now  exists  a  3rd  (1977)  edition  of  ENW.  This  is  poorer  still  fi'om  the  standpoint  of 
describing  what  actually  was  observed  during  a  nuclear  explosion,  because  by  1977  the  observations  of 
aunospheric  nuclear  tests  of  1962  and  earlier  had  been  massaged  to  provide  a  homogenized  but  not 
necessarily  more  accurate  description  of  what  goes  on. 
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Figure  6.4.  Variation  of  Apparent  Fireball  Surface  Temperature  with  Time 
in  a  20*kt  Explosion.  (Source:  ENW,  1964,  p.  75) 
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6.2.3  A  Simple  Model  of  Low- Altitude  Fireballs  (Source:  Ewing  & 
Barrett,  1989) 

6.2.3. 1  Fireball  Size  and  Motion 

This  discussion  assumes  spherical  fireballs  (i.e.,  ignores  the  development  of  a 
toroid),  applies  for  air  bursts^^  below  50  km  burst  height  and  to  times  up  to  60  seconds 
after  burst. 

A  simple  estimate  for  the  radius  is 


R(t) 

= 

0.85  YO-33  (po/p)0.3  +  0.037  YO-3  (po/p)0-2  (t-t*)  (km) 

where 

Y  = 

total  weapon  yield  (Mt) 

Po  = 

air  density  at  sea  level 

P  = 

air  density  at  burst  altitude 

t  = 

time  after  burst  (sec) 

fpe  = 

(Y  po/p)^^  (sec) 

t*  = 

tpe  +  Y  (sec). 

Eq.  (6.6)  applies  for  times  greater  than  t*;  note  that  tpe  is  an  approximation  to  the  time  of 
pressure  equilibrium,  i.e.,  at  earlier  times  there  are  major  pressure  forces  giving  rise  to 
vigorous  motions.  The  first  term  in  Eq.  (6.6)  is  the  initial  fireball  radius,  the  second 
accounts  approximately  for  expansion.  An  approximate  expression  for  the  riye  velocity  of 
low  altitude  fireballs,  again  for  the  first  60  sec  after  burst,  is 

Vrise  =  0.1 1  (po/p)0.15  (km/sec)  (6.8) 


6. 2. 3. 2  Blast 

The  partitioning  of  energy  between  blast  and  thermal  output  begins  to  change  above 
10  km  altitude,  and  thus  a  simple  blast  efficiency  factor 

BefT  =  Max[  Mind;  1.15 -0.015  H);0]  (6.9) 

v.'here  H  (km)  =  height  of  burst.  Thus  Beff  =  1.0  at  H=  0,  0.85  at  20  km,  0.4  at  50  km, 
and  0  for  H  >  77  km. 


An  air  burst  is  one  detonated  so  high  that  the  fireball  does  not  interact  with  the  ground;  in  the 
following  Section  (6.2.4)  we  discuss  bursts  that  do  interact  with  the  ground. 
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Peak  overpressure  (traditionally  expressed  in  psi)  is  obtained  by  scaling  the 
effective  blast  yield  BeffY  (Mt)  to  a  standard  1  kt  explosion  at  sea  level.  Figure  6.5  is  a  plot 
of  this  standard  free-air  burst  taken  from  Glasstone  &  Dolan,  1977;  the  dashed  curve  is  a  fit 
of  the  functional  form  that  follows  [Eq.(6.10)]. 


With  horizontal  range  from  the  burst  expressed  in  km,  peak  overpressure  for  the 
standard  burst  can  be  approximated  by 


where 


[18.7/(log  Dst+ 3.72) -5.0] 


PST 

=  10 

(6.10) 

PST 

=  overpressure  for  1  kt  at  sea  level  (psi) 

Dst 

=  horizontal  range  at  sea  level  (km) 

Dst 

=  D/YscaLe1/3  (km) 

(6.11a) 

D 

=  actual  horizontal  range  (km) 

YsCALE 

=  scale  blast  yield  (kt) 

YsCALE 

=  1000  Bgff  Y/pamb  (^^0 

(6.11b) 

Pamb 

=  ambient  pressure  (atm.)  =  2.83  Tamb(K)  Pamb(g/cro3)  . 

This  formulation  assumes  that  the  ambient  pressure  is  approximately  constant  between  the 
burst  altitude  and  the  altimde  of  interest 


Overpressure  at  the  point  of  interest  is 

P  =  PSTPamb  (psi)  (6.12) 

Dynamic  pressure  (=  kinetic  energy  per  unit  volume  of  air  just  behind  the  shock 
front)  is^^ 

q  =  (5  p2/2)/(7  X  14.7  pamb  +  P)  (6. 13) 

and  the  shock  velocity  at  any  altitude  is 

V  =  a  {  1  +  [6p/(7  X  14.7  Pamb)]  (nt/sec)  (6.14) 

where  a  =  sound  speed  at  this  altitude  =  340  (To/T)^/2  (m/sec)  if  Tq  =  temperature  at  sea 
level  and  T  =  temperature  at  relevant  altitude. 


This  formula  assumes  that  y  =  Cp/Cv  =  1.4  and  that  sea  level  pressure  is  14.7  psi. 
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DISTANCE  FROM  BURST  (FEET) 


Figure  6.5.  Peak  Overpressure  from  a  1  Klloton  Free  Air  Burst  for  Sea  Level 
Ambient  Conditions.  [Solid  curve-Giasstone  &  Dolan,  1977,  p.lOS; 

dashed  curve-Eq.(6.10)] 
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6. 2. 3. 3  Thermal  Fluence 


At  altitudes  up  to  about  250  km,  this  can  be  expressed  as 
Ruence  =  10^  f  Y  t/(47cR2)  (cal/cm^) 

where  x  =  atmospheric  transmission 


(6.15) 


R  =  range  from  burst  point  to  field  point  (km) 

f  =  fraction  of  total  weapon  energy  emitted  as  thermal  radiation,  which 


can  be  approximated 

as: 

f  = 

0.5 

for  0  <  H  <70  km 

(6.16a) 

f  = 

1.2-H/lOO 

for  70  km  <H  <  115  km 

(6.16b) 

f  = 

0.05 

forll5km<H  <250  km  . 

(6.16c) 

6.2.4  Near>Surface  Bursts 

When  a  bomb  is  detonated  so  close  to  the  ground  that  the  explosion  interacts  with 
the  ground,  a  crater  will  be  formed  and  a  large  fraction  of  the  energy  goes  to  produce  this 
crater  and  to  send  a  shock  wave  through  the  ground.  Because  the  earth  density  is  larger 
(by  a  factor  ~  1000)  than  the  air  density,  the  crater  radius  will  be  significantly  smaller  than 
the  fireball  radius;  thus,  for  a  1  Mt  near  (but  below!)  surface  burst  the  radius  of  the  crater  is 
on  the  order  of  200  m,  while  the  radius  for  a  1-Mt  airburst  at  thermal  minimum  is  ~  400  m 
(see  ENW,  1977).  This  ground  shock  due  to  a  surface  burst  will  lead  to  structural  damage. 
As  the  height  of  burst  decreases  so  that  the  interaction  of  the  explosion  with  the  ground 
becomes  stronger,  more  of  the  energy  goes  to  produce  a  crater  and  to  send  a  shock  through 
the  ground.  For  a  discussion  of  ground  and  subsurface  bursts,  see  ENW,  1977,  Ch.  6. 

When  a  bomb  is  detonated  so  close  to  the  surface  that  the  fireball  touches  the 
ground,  a  crater  will  be  formed  and  a— possibly  large— dust  cloud  will  be  sent  into  the 
atmosphere.  Figure  6.6  shows  the  anticipated  mass  of  dust  that  is  lofted  for  a  given  yield 
as  a  function  of  the  "scaled  height  of  burst"- which  scales  with  yield  as  Figure  6.6a 
shows  the  data  base,  and  Fig.  6.6b  a  currently  accepted  model.  Typically,  the  dust 
particles  will  be  quite  small-mainly  less  than  1  jim— but  the  results  are,  of  course,  highly 
variable,  depending  on  the  local  geology. 

There  are  very  few  larger  panicles  (>  1-10  |im),  but  they  are  particularly  important 
for  public  health  because  a  large  fraction  of  the  radioactivity  produced  mainly  by  fission 


ANHYDROUS  NASS  LOADING  (ML  OUST/Nt  UCATON  YIELD) 


Figure  6.6.  Dust  Mass  Lofted  into  the  Stabilized  Cloud  by  a  Near-Surface 
Nuclear  Explosion.  (Source:  Rausch  et  al.,  1988) 

(a)  Data  Base,  (b)  Model. 
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and  by  neutron  activation  condenses  on  their  surface.  These  large  particles  drift  in  the  wind 
and  are  eventually  deposited  as  radioactive  fallout,  which  produces  a  serious  long-term 
health  hazard  for  people  (see  ENW,  1977,  Ch.  9  and  12). 


6.2.5  Damage  Mechanisms  in  the  Lower  Atmosphere 


In  the  lower  atmosphere  a  nuclear  explosion  creates  a  fireball  which  produces  both 
a  strong  blast  wave  that  carries  some  60-65  percent  of  the  total  energy  of  the  detonation 
and  is  a  major  source  of  damage  and  also  thermal  radiation,  which  typically  carries 
30-35  percent  of  the  energy  of  a  near-surface  burst.  Nuclear  radiation,  principally  neutrons 
and  gamma-rays,  both  "prompt"  and  "delayed,"  produce  some  damage,  while  some  long¬ 
term  injuries  to  people  are  caused  by  radioactive  fallout. 


Some  representative  damage  radii  on  the  ground  for  a  1  Mt  bomb  detonated  at  an 
altitude  of  1  km  above  the  surface  are  the  following: 

•  Fireball  radius  at  thermal  minimum^*  0.44  km  (ENW,  1977,  p.  70) 

•  Fireball  radius  at  breakaway^^  0.53  km  (ENW,  1977,  p.  70) 

•  Peak  dynamic  pressure  at  the  surface  (ENW,  1977,  p.  1 17): 

1001b/in2  1.7  km 

101b/in2  2.4  km 

1  lbAn2  4.3  km 

•  Radiant  energy  (~  35  percent  of  energy)  (ENW,  1977,  p.  291): 

30  cal/cm2  7  km 

3cal/cm2  16  km 


•  Neutrons  (prompt  and  delayed)  (ENW,  1977,  p.  347): 

300  rad  (tissue)  2. 1  km 

30  rad  (tissue)  2.6  km 

•  Gammas  (prompt  and  delayed)  (ENW,  1977,  p.334): 

300  rad  (tissue)  2.7  km 

30  rad  (tissue)  3.3  km 


For  more  detail  and  an  indication  of  damage  levels,'*®  see  ENW,  1977,  or  EM-1. 


3  ^  Defined  in  Section  6.2.2. 

Defined  in  Section  6.2.2. 

The  rad  is  a  unit  of  absorbed  radiation,  characterized  by  the  absorbing  material  (tissue  for  biological 
effects.  Si  for  electronics,  etc.).  1  rad  =  100  erg  nuclear  or  ionizing  radiation  per  gm  absorbing 
material.  See  also  footnote  **  on  Table  6.2. 


6.2.6  Electromagnetic  Pulse  (EMP) 


When  a  nuclear  weapon  is  detonated,  a  large  amount  of  ionization  is  produced  very 
rapidly  from  the  deposition  of  X-rays  and  y-rays,  in  a  time  of  order  lO"^  sec.  If  the 
ionized  region  is  not  spherically  symmetric,  typically  as  a  result  of  the  effects  of  the  ground 
or  of  the  variation  of  atmospheric  density  with  height  (~  7  km  scale  height),  the 
electromagnetic  field  produced  by  this  ionization  produces  a  large  RP  radiated  signal,  the 
Electromagnetic  Pulse. 

The  character  of  this  field  varies  significantly,  depending  on  the  presence  of  the 
ground  (conducting  ground  plane)  and  on  the  ambient  atmospheric  density.  Thus,  for 
surface  and  very  low  air  bursts  the  maximum  field  strength  is  on  the  order  of  keV/m, 
extending  out  to  some  10-20  km  from  the  burst  because  most  of  the  X-rays  and  y-rays  are 
absorbed  in  this  region.  For  intermediate  altitude  air  bursts  (20-50  km  height  of  burst)  for 
which  there  is  no  ground  effect,  the  field  strength  is  much  smaller,  in  the  range  of  tens  to 
hundreds  of  eV/m,  and  with  an  extent  which  may  be  slightly  larger  because  of  the  lower 
ambient  density.  However,  for  high-altitude  bursts,  where  the  ambient  air  density  is  so 
low  that  no  fireball  is  formed,  the  maximum  field  strength  can  again  be  in  the  range  keV/m, 
but  the  geometrical  extent  can  be  very  large.  Thus,  for  a  Mt- yield  burst  at  500  km  altitude, 
where  the  ionization  peaks  at  altitudes  of  100-200  km,  there  will  be  large  electromagnetic 
fields  on  the  surface  of  the  earth  out  to  the  geometrical  tangent  altitude  range,  which  here 
would  extend  some  2500  km  from  the  burst  point. 

These  electromagnetic  fields  can  have  a  very  deleterious  effect  on  unprotected 
electronics,'^^  but  these  effects  can  be  minimized  by  proper  design  (screening). 

For  a  discussion  of  EMP,  see,  e.g., 

•  Ch.ll  of  ENW,  1977. 

•  IEEE  Transactions  on  Electromagnetic  Compatibility,  Vol.  EMC-20,  No.  1. 
part  1,  Feb.  1978  (also  published  as  IEEE  Trans,  on  Antennas  and 
Propagation,  Vol.  AP-26,  No.  1,  Jan.  1978). 

•  Soper,  1985-for  a  Classified  description. 

A  related  phenomenon  is  System-Generated  Electromagnetic  Pulse  (SGEMP), 
which  refers  to  the  transient  electromagnetic  pulse  created  by  the  electrons  emitted  from  the 


Or  even  on  (long)  power  lines;  note  that  following  a  high-aliitude  nuclear  explosion  in  the  South 
Pacific,  several  strings  of  street  lights  in  Hawaii,  some  1600  km  from  the  burst,  were  affected.  See 
Appendix  H  (J.  Malik  from  Hoerlin,  1976). 
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surface  of  an  electronic  or  other  system  when  this  is  exposed  to  incident  photons  from  a 
nuclear  burst.  SGEMP  differs  from  "ordinary"  EMP  in  that  the  source  of  electrons  that 
produce  the  field  is  the  system  itself,  rather  than  the  intervening  air.  For  a  review,  see 
Higgins  et  al.,  1978. 


6.3  HIGH-ALTITUDE  NUCLEAR  EXPLOSIONS 


6.3.1  Phenomenology 


6.3. 1.1  Introduction 

The  low-altimde  analysis  of  Section  6.2  applies  generally  for  heights  of  burst  below 
30  km.  At  intermediate  altitudes  (30-60  km)  the  fireball  is  larger  and  more  diffuse,  but  the 
effects  are  still  quite  confined.  The  discussion  of  Fig.  6.1  and  Table  6.3  in  Section  6.1 
indicates  that  from  the  standpoint  of  initial  absorption  of  energy  (which  is  not  necessarily 
the  only  factor)  there  will  be  a  confined  fireball  at  altitudes  below  approximately  60-80  km. 

At  higher  altitudes  things  are  quite  different  Quoting  from  ENW,  1977,  p.74: 

For  burst  heights  up  to  ~  80  km  the  early  fireball  is  approximately  spherical 
(it  later  becomes  toroidal),  but  at  the  higher  altitudes  it  elongates  venically. 

Above  80  km  the  debris  tends  to  be  separate  from  the  X-ray  pancake.  The 
debris  can  rise  to  ^eat  heights,  unlike  the  X-ray  pancake  which  is  always  at 
the  same  height  (since  it  forms  where  downward-traveling  X-rays  run  into 
the  absorbing  atmosphere). 

For  burst  altitudes  above  80  km  there  is  virtually  no  absorption  of  X-rays 
emitted  in  upward  directions.  The  downward  directed  X-rays  are  mostly 
absorbed  in  a  layer  of  air,  the  "X-ray  pancake"  which  becomes  incandescent 
as  a  result  of  energy  deposition.  Tlus  pancake  is  more  like  the  frustum  of  a 
cone  pointing  upward,  10  km  thick  at  a  mean  altitude  of  80  km;  the  radius  at 
this  height  is  approximately  the  HOB  minus  80  km. 

6.3. 1.2  Description  of  W.W.  White,  1989 

The  high-altitude  regime  (above,  say,  100  km)  is  characterized  by  low  air 
mass  density.'*^  in  this  rarefied  environment,  the  mean  free  paths  of  X- 
rays,  UV  photons  and  energized  air  particles  (excluding  ion-magnetic 
coupling)  are  large.  Thus  high-altitude  nuclear  explosion  (HANE) 
disturbances  cover  large  geographic  areas.  For  example,  the  effects  from  a 
megaton-class  explosion  at  200  km  altitude  extend  roughly  700  km 


Air  mass  density  at  100  km  is  more  than  a  million  times  smaller  than  at  the  earth’s  surface. 
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horizontally  out  from  the  burst  point—see,  e.g.,  Fig.6.7/^  In  addition,  Ac 
low  electron-ion  collisional  radiative  recombination  rate  in  this  low-density 
regime  results  in  long-lasting  (tens  of  hours),  widespread  regions  of 
persistent  ionization.^ 


Figure  6.7.  Phenomenology  of  High-Yield  Burst  at  200  km  Altitude  and 
Midlatitude.  A  detailed  Discussion  of  this  Figure  is  given  In  Section  6.3.1. 3. 


The  individual  elements  of  Fig.  6.7  are  discussed  in  the  following  Section  6.3. 1.3. 

Dissociation  due  to  burst  energy  deposition  makes  molecular  deionization  in  these  regions 
unimportant 


Evolution  of  the  disturbed  high-altitude  atmosphere  is  usefully  divided  into 
early-time  and  late-time  regimes,  as  described  below. 

Earlv-time  Regime.  This  refers  to  the  first  minute  or  so  following  a  burst. 
During  this  period  ,  burst  energy  deposition  and  initial  hydrodynamic 
expansion  of  the  disturbed  atmosphere  occur. 

Energy  Deposition.  Roughly  25  percent  of  the  energy  released  by  a  modem 
nuclear  weapon  appears  as  rapidly  expanding  weapon  debris  (kinetic  yield), 
and  75  percent  appears  as  X-rays.  A  small  fraction  of  the  burst  energy  is 
carried  by  neutrons,  betas,  and  gammas— see  Table  6.1.  The  kinetic  yield  is 
partition^  into  a  mix  of  UV  taxation  and  energized  air  and  debris  ions;  the 
relative  proportion  is  mainly  a  function  of  burst  point  air  density.  The 
energized  ions  move  up  and  down  the  (disturbed)  geomagnetic  field  lines— 
see  Fig.  6.8-and  interact  with  the  atmosphere  to  produce  "patch"  regions  of 
enhanced  ionization  and  temperature.  For  mid  to  low  latitude  bursts 
(outside  the  auroral  oval)  deposition  occurs  in  both  magnetic  conjugates. 
For  a  megaton-range  burst,  most  of  the  kinetic  yield  appears  as  UV 

radiation  for  burst  point  air  densities  above  about  2  x  10-13  gm/cm^ 
(corresponding  to  approximately  200  km  altitude  in  ambient  air).  The 
fraction  of  kinetic  yield  appearing  as  energized  ions  increases  rapidly  with 
decreasing  burst  point  air  density. 

We  distinguish  three  basic  high-altitude  nuclear  burst  deposition 
phenomenology  regimes,  according  to  how  the  weapon  debris  interacts  with 
the  atmosphere. 

For  burst  point  air  densities  above  about  2  x  10-13  gin/cm3  (below  about 
200  km  altitude  in  the  ambient  atmosphere),  collisions  between  air  and  fast 
weapon  debris  ions  form  a  debris-air  "piston".  The  piston  sweeps  up  air  as 
it  propagates  outward,  and  through  compression  converts  debris  kinetic 
energy  into  thermal  energy  at  the  shock  front.  Excited  ions  emit  UV 
photons,  a  significant  fraction  of  which  ionize  the  surrounding  air,  creating 
the  UV  fireb^l.  Such  bursts  are  in  the  UV  REGIME. 

At  lower  burst  point  air  densities  (e.g.,  10-15  gm/cm3,  corresponding  to 
about  500  km  altitude  in  ambient  atmosphere),  collision  frequencies  are  too 
small  to  form  a  UV-radiating  blast  wave.  Instead,  debris  kinetic  energy  is 
transferred  to  air  ions  by  a  collisionless  process  mediated  by  the  geomagnetic 
field.  As  ionized  weapon  debris  moves  outward  it  carries  along  tiie  frozen-in 
geomagnetic  field.  This  compresses  the  field  ahead  of  the  debris  and 
energizes  the  air.  The  debris  kinetic  yield  is  thus  partitioned  between  energy 
of  the  air  ions  and  the  deformed  geomagnetic  field  Energized  air  ions  move 
downward  (and  upward)  along  the  geomagnetic  field  until  they  are  stopped 
by  collisions  with  air.^5  Bursts  characterized  by  such  energy  deposition  are 
in  the  KINETIC  ENERGY  PATCH  (KEP)  REGIME.  For  bursts  at  mid  and 
low  magnetic  latitudes,  patch  formation  is  expected  at  both  northern  and 
southern  magnetic  conjugate  regions. 


^5  In  the  100-200  km  altitude  range  for  a  burst  in  ambient  air. 
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At  still  lower  burst  point  air  densities  (below  10“1®  gin/cm^;  the 
corresponding  altitude  in  the  ambient  atmosphere  varies  firom  500  to  800  km 
depending  on  solar  activity),  almost  all  the  debris  kinetic  yield  is  converted  to 
magnetic  energy  of  the  deformed  geomagnetic  field  (the  air  density  is  too 
small  to  significantly  drain  energy  from  the  collisionless  blast  wave). 
Current  understanding  in  this  area  is  somewhat  provisional,  however,  and 
the  physical  processes  likely  to  be  operative  are  still  subjects  of  active 
research.  Furthermore,  no  nuclear  tests  were  ever  conducted  in  this  density 
regime,  so  direct  observational  guidance  is  not  available.  Neveretheless,  it  is 
thought  that  perhaps  half  the  debris  kinetic  energy  will  be  contained  and 
guided  to  the  magnetic  conjugates  by  the  geomagnetic  field.  Formation  of 
ionized  regions  similar  to  IfflPs  is  anticipated.  Depositions  of  this  type  are 
in  the  MAGNETIC  CONTAINMENT  REGIME  (MCR). 

Of  the  three  energy  deposition  regimes,  UV,  KEP,  and  MCR,  UV  is  best 
understood,  and  MCR  least  understood.  In  general,  divisions  between  these 
regimes  are  not  sharp;  an  explosion  will  generally  be  characterized  by  a 
mixture  of  phenomenologies. 

Burst  energy  deposition  takes  place  in  a  short  but  finite  time  (~  10“^- 10^ 
sec).  PROMPT  QffiMISTRY  also  occurs  concurrently  with  the  deposition. 
Processes  include  photodissociation,  photoionization,  charge  exchange, 
coUisional  excitation,  and  so  on. 

Hydrodynamic  Expansion.  In  the  high-altitude  regime,  the  vertical  extent  of 
the  initially  heated  region  is  much  larger  than  the  atmospheric  density  scale 
height.  tIius  the  hydrodynamic  response  to  the  burst-induced  pressure 
gradients  results  in  predominantly  upward  expansion  of  the  nuclear- 
disturbed  air  ("ATMOSPHERIC  HEAVE"-see  Section  6.3.3  ). 

The  relatively  high  temperatures  and  densities  in  the  main  deposition  region 
lead  to  strong  ion-neutnd  coUisional  coupling  in  the  initial  expansion,  ’^us 
the  early-time  expansion  can  be  approximately  described  by  single-fluid 
hydrodynamics.  In  addition,  high  gas  kinetic  pressures  produce  a  high-beta 
plasma**^  in  which  the  magnetic  field  is  carried  along  with  the  hydrodynamic 
flow. 

As  the  disturbed  atmosphere  expands  and  cools,  ion-neutral  coUisional 
coupling  weakens,  plasma  beta  decreases,  and  the  geomagnetic  field  relaxes 
towards  its  ambient  configuration.  After  about  one  minute  a  TWO  FLUID 
MHD  description  of  the  evolving  atmosphere  is  required. 

Also  during  this  period  DEIONIZATION  CHEMISTRY  acts  to  reduce  the 
high  level  of  initial  ionization  pnxluced  by  the  explosion.  Where  molecules 
are  abundant,  recombination  proceeds  rapidly  by  molecular  dissociative 


^  Beta  =  (gas  iHessure)/(magnetic  pressure3^/87t)  >  1 . 


recombinancm/^  Deionization  also  occurs  by  direct  electron-ion  coilisional- 
radiative  recombination. 

Late-Time  Regime.  In  the  LATE-TIME  PHASE,  ions  and  electrons  tend  to 
move  along  the  geomagnetic  field  (now  relaxed  to  near-ambient)  while  the 
neutrals  are  not  directly  affected  by  the  field.  Plasma  and  ioniz;^  weapon 
debris  stream  along  field  lines  to  very  high  altitudes.  A  weak  ion-neutral 
coupling  permits  the  two  fluid  components  to  slip  through  one  another,  and 
a  two  fluid  MHD  approach  is  required  to  adequately  treat  the  late-time 
physics. 

Superposed  on  the  primary  parallel-to-B  motion  is  an  ExB  drift  of  the 
plasma  perpendicular  to  the  field.  The  electrical  field  is  generated  by 
collisional  coupling  of  the  plasma  with  the  neutral  wind,,  the  inertial  force 
due  to  plasma  motion  along  the  curved  geomagnetic  field,  gravitation,  and 
other  less  important  drivers.  Because  the  parallel-to-B  electrical 
conductivity  is  large,  individual  magnetic  flux  tubes  are  nearly  equipotential. 

This  produces  an  ExB  drift  which  preserves  field-aligned  plasma  structures 
and  results  in  a  nearly  incompressible  (acmally,  magnetic-flux-conserving) 
flow.  For  bursts  in  the  northern  hemisphere,  overall  ExB  plasma  drift  is 
northward. 

Fluid  and  kinetic  instabilities  in  the  nuclear-generated  plasma  cause  density 
irregularities  to  develop.  This  spatial  structure  involves  a  broad  spectrum  of 
scale  sizes  covering  many  orders  of  magnitude.  By  late  times  these 
structures  extend  thousands  of  kilometers  along  magnetic  field  lines  forming 
a  system  of  FIELD- ALIGNED  STRIATIONS.  These  field-aligned  density 
irregularities  can  severely  degrade  radar  and  satellite  communication  link 
performance. 

During  the  late-time  phase  the  plasma  plume  rises  to  very  high  altitudes.  Plasma 
density  decreases  largely  through  expansion;  below  10^  electrons/cm^,  recombination 
proceeds  slowly  and  large  regions  of  substantial  ionization  persist  to  very  late  times. 

The  overall  picture  described  here  presents  the  large-scale  structure  at  various  times. 
Note  that  this  picture  is  not  unifoim  and  static: 

(i)  there  are  substantial  winds  and  the  plume  rises  quite  rapidly 

(ii)  there  is  a  great  deal  of  small-scale  non-uniformity.  Thus  as  the  bomb  comes 
apart  there  is  small-scale  turbulent  structure,  and  then  striations  (on  a 
transverse  scale  of  1-10  km)  develop  in  the  first  1-10  sec  and  last  for  long 
times,  ~  30-60  minutes  or  even  longer. 


Ion-molecule  reaction  to  produce  molecular  ions  NO"^  and  02'*’  which  have  a  lower  ionization  potential 
than  the  atomic  ions  N"*"  or  O"^,  followed  by  the  rapid  dissociative  recombination  reactions  NO"*"  +  c'  => 
N  +  O"  and  02"^  c  =>  O'  +  O",  whe.'e  the  primes  denote  possible  electronic  excitation  of  the 
resulting  atoms. 


6-24 


It  must  be  stressed  that  the  details  of  this  structure  represent  the  result  of  20  years  of 
synthesis  of  rather  limited  data,  and  of  computer  modeling  and  analog  experiments  (Barium 
cloud  releases,  Auroral  studies,  etc.). 

The  debris  clouds  from  high-altitude  nuclear  explosions  rise  by  buoyancy  and 
expand  as  they  do  so.  They  rise  much  more  than  do  the  fireballs  from  low-altitude  air 
bursts  because  of  the  very  much  lower  ambient  density.  High-altitude  explosions  send 
their  ionized  debris  a  very  long  way,  typically  along  the  geomagnetic  field  lines.  The 
debris  produces  a  large  number  of  secondary  electrons  which  can  give  rise  to  dramatic 
auroral  displays,  and  some  of  which  can  also  be  trapped  by  the  geomagnetic  field  to 
enhance  the  (Van  Allen)  radiation  belts-see,  e.g.,  Hoerlin,  1976. 

6.3. 1.3  Discussion  of  the  Different  Features  of  Fig.  6.1^^ 

Figure  6.7  shows  the  overall  phenomenology  regions  according  to  current  nuclear 
models.  This  figure  (which  applies  to  a  Megaton-yield  bomb  detonated  at  around  200  km 
at  mid-latitudes— results  are  typical  for  t  =  1-10  sec  after  burst)  is  very  instructive,  and  the 
individual  features  are  discussed  here. 

(a)  X'foy  patch  at  70  km  altitude  with  horizontal  extent  of  up  to  1000  km.  This  is 
produced  by  the  absorption  of  X-rays  emitted  from  the  bomb  in  a  generally 
downward  direction  where  they  are  absorbed  by  the  atmosphere;  X-rays 
emitted  in  an  upward  direction  are  generally  not  absorbed  in  the  very  diffuse 
atmosphere  at  higher  altitudes. 

(b)  Beta-patch  (somewhat  smaller  than  the  X-ray  patch,  but  highly  structured) 
produces  ionization  down  to  65  km,  of  horizontal  extent  ~  100  km.  This  is 
produced  essentially  fiom  the  fission  P-particles  in  the  delayed  radiation. 

(c)  Much  of  the  weapon  debris  will  be  ionized,  especially  at  these  low  densities, 
and  thus  it  will  tend  to  follow  the  geomagnetic  field  lines.  Thus  at  different 
latitudes  its  behavior  may  be  quite  different.  A  conductive  bubble  is  formed 
which  expels  the  geomagnetic  field  lines.  It  creates  a  structured  blast  wave 
whose  boundaries  may  be  unstable  due  to  a  variety  of  mechanisms  (Rayleigh- 
Taylor  and  Field  Curvature  instabilities  are  indicated  at  different  regions  of 
space).  As  a  result  of  these  instabilities,  eventually  the  bubble  will  break  up, 
giving  rise  to  organized  and  irregular  structures  (striations  and  turbulence) 
which  eventually  lead  to  turbulent  mixing. 


48  Sections  6.3.12  and  6.3. 1.3  provide  somewhat  different  perspectives  to  a  very  complex  situation. 
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(d)  The  general  motion  of  the  debris  goes  along  the  geomagnetic  field  lines: 
downward-moving  debris  will  encounter  the  denser  atmosphere  and  thus  will 
be  stopped  and  deposited  as  the  debris  energy  patch,  while  the  upward  moving 
debris  produces  a  plume  which  follows  the  geomagnetic  field  lines  (see 
Fig.  6.8). 

(e)  As  the  magnetic  bubble  expands,  it  generates  UV  radiation  which  produces 
ionization  in  a  uniform  way,  but  in  addition  softer  (non-ionizing)  UV  is 
deposited  in  the  dense  lower  atmosphere,  leading  to  local  heating. 

(f)  Notice  that  some  of  these  regions— those  produced  just  by  electromagnetic 
radiation— are  relatively  uniform,  while  those  that  are  derived  from  charged 
particles  show  structure  due  to  interactions  with  the  geomagnetic  field, 
including  some  instabilities. 

These  are  clearly  very  large-scale  structures  which  survive  for  rather  long  times. 
The  bubble  of  Fig.  6.7  is  typical  of  conditions  ~  1-10  sec  after  burst,  while  the  global  scale 
processes  of  Fig.  6.8  last  for  minutes  to  hours. 

6.3.2  Multiburst  Phenomenology 

At  high  altitudes  where  the  density  is  low,  the  energy  from  a  nuclear  explosion 
travels  a  long  distance  before  being  absorbed  in  the  atmosphere,  and  thus  there  is  the 
possibility  for  interaction  between  neighboring  bursts.  In  terms  of  current  nuclear  arsenals 
of  both  the  US  and  USSR,  the  possibility  of  massive  multiburst  scenarios,  involving  some 
hundreds  to  thousands  of  nuclear  explosions  sufficiently  close  together  (in  space  and  time) 
so  that  interaction  can  occur,  cannot  be  excluded.  Thus  we  ask  in  what  way  such 
multiburst  scenarios  give  rise  to  different  phenomenology  from  that  for  the  case  of  a  single 
burst  as  discussed  in  Section  6.3.1. 

Figure  6.9  gives  the  results  of  a  "Massive  Multiburst"  computation  which  has  been 
performed  recently  (1986-87)  by  MRC  for  DNA  using  the  SCENARIO  code.  A  total  of 
1000  one-Mt  bursts  are  detonated  at  a  rate  of  10  every  10  sec  along  three  ICBM 
trajectories,  corresponding  to  minimum  energy,  lofted  and  elevated  trajectories,  and  a 
spacing  so  close  that  there  is  significant  interaction,  and  yet  far  enough  away  so  that 
fratricide  can  be  avoided.  Figure  6.9a  shows  the  horizontal  spacing  of  the  bursts. 
Figs.  6.9b  and  6.9c  show  the  toul  mass  density  and  electron  density  1  min  into  the 
laydown,  i.e.,  after  60  bombs  have  been  detonated,  while  Figs.  6.9d  and  6.9e  show  the 


total  mass  density  and  electron  density  after  the  last  of  the  weapons  have  been  detonated. 
The  locations  of  the  bursts  are  shown  in  all  these  figures.^^ 


(continued) 

Figure  6.9.  Massive  Multiburst,  Case  MS-1.1  (Source:  White  et  al., 1987a). 

(a)  Vertical  view  showing  the  location  of  the  bursts  over  the  earth; 

(b)  Total  mass  density  at  t  =  1  minute  (60  bursts); 

•  (c)  Electron  density  at  t  =  1  minute; 

(d)  Total  mass  density  at  t  s  16.7  minutes  (1000  bursts); 

(e)  Electron  density  at  t  s  16.7  minutes. 


9  The  quantitative  results  of  this  simulation  are  now  (mid-1990)  considered  obsolete,  although  the 

general  characteristics  aren't  changed.  You  should  ask  DNA/RAAE  for  up-to-date  results. 
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Figure  6.9.  (continued) 


Note  that  the  ambient  ionosphere  is  modified  drastically:  the  ionization  is  enhanced 
by  a  factor  up  to  10^  and  the  neutral  density  is  enhanced  by  a  comparable  or  even  larger 
factor  at  the  higher  altitudes.  The  enhanced  ionization  can  be  understood  quite  simply  in 
terms  of  the  ionization  of  ambient  atoms  and  molecules  by  incident  X-rays  and  UV 
radiation.  The  enhancement  in  neutral  density  is  induced  by  the  buoyant  rise  of  air  heated 
by  the  absorption  of  X-rays  or  debris  kinetic  energy.  This  phenomenon  is  called  Heave 
and  is  discussed  in  Section  6.3.3.  At  the  moment  let  us  simply  note  that  Fig.  6.10  shows 
the  enhancement  in  total  density— both  mass  density  (kg/m^)  and  number  density 
(no./cm^)— due  to  multiple  nuclear  explosions  [cases  P(60)and  H(IOOO)  of  Fig.  6.9]  and 
Fig.  6. 1 1  shows  the  plasma  temperature.^®  Both  these  figures  are  discussed  in  Section 
6.3.3. 

The  detailed  numbers  should  not  be  taken  too  literally;  note  that  the  enhancement  in 
total  and  ionization  density  is  very  large  but  suggests  saturation  behavior,  i.e.,  the  electron 
density  does  not  exceed  10^0-10^^  cm"^  over  large  regions  of  space  while  the  density 
enhancement  below  ~  200  km  is  very  small.  The  heave  phenomenon  will  be  discussed 
next,  in  Section  6.3.3. 

Applications  for  system  performance  are  discussed  in  Section  6.4;  the  ionization 
(which  affects  electromagnetic  signal  propagation)  is  discussed  in  Section  6.4.2,  and  the 
enhancement  in  radiation  ("Redout")  is  taken  up  in  Section  6.4.3.  Both  these  quantities  are 
functions  of  the  plasma  temperature.  Finally,  the  overall  quality  of  the  predictions  and  its 
implications  on  system  performance  are  taken  up  in  Section  6.5. 

6.3.3  Heave 

McNamara  (1935)  describes  atmospheric  heave  as  "the  comparatively  slow, 
generally  adiabatic  motion  of  large  regions  of  the  upper  atmosphere  following  energy 
deposition  by  a  high-altitude  nuclear  detonation."  This  definition  needs  some  qualification 
and  amplification  because  the  upward  speed  increases  greatly  with  altitude,  and  because  the 
motion  of  a  heated  air  parcel  of  radius  R  becomes  ballistic  rather  than  buoyant  at  altitudes 
above  about  50  km-see,  e.g.,  Sowle,  1977,  pp.  480f  and  505f. 


5®  Above  100  km  the  atmosphere  is  significantly  out  of  LTE  (Local  Thermodynamic  Equilibrium)  so  that 
it  is  not  meaningful  to  speak  of  a  Temperature.  However,  one  can  define  a  Plasma  Temperature  in 
terms  of  the  degree  of  ionization  of  the  ambient  atmosphere-  -cf.,  e.g..  Part  1  of  these  Notes,  Section  2. 
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Figure  6.11.  Enhanced  Plasma  Temperature  due  to  Multiple  Nuclear 
Explosions  [Calculations  for  the  1000-burst  case  MS  at 
(a)  t  =  1  min,  (b)  t  =  2  min,  and  (c)  t  =  16.4  min] 
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Figure  6.11.  (continued) 
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In  an  atmosphere  at  constant  temperature  T,  the  density  p(z)  falls  off  with 
increasing  altitude  z  as  follows: 

p(z)  =  p(zo)  exp  -  (z  -  Zo)/H  ,  (6.6a) 

where  the  scale  height  H  is  given  by  the  expression 

H  =  kT/Mg  .  (6.6b) 

H  ~  7  km  at  altitudes  below  100  km  where  the  temperature  is  ~  200-250  K;  above  200  km 
altitude,  where  T  ~  700-1500  K  and  the  atomic  oxygen  is  dissociated  so  that  the  effective 
molecular  weight  is  typically  18  (rather  than  29),  H  ~  50  km. 

An  air  parcel  of  radius  R  is  said  to  rise  buoyantly  if  R  «  H,  and  ballistically  if 
R  »  H.  Buoyant  rise,  which  occurs  at  relatively  low  altitudes,  corresponds  to  relatively 
slow,  adiabatic  rise  of  the  air  parcel  which  maintains  a  uniform  pressure,  while  ballistic 
rise,  which  takes  place  lower  densities  is  relatively  rapid  with  pressure  varying  in  the 
parcel. 

For  definiteness,  consider  bursts  above  300  km,  and  ask  for  the  motion  as  a 
function  of  altitude.  Most  of  the  energy  of  a  bomb  comes  out  as  X-rays  (see  Table  6.1) 
and  downward-traveling  keV  X-rays  are  absorbed  near  80  km  -  see  Table  6.3-  where  the 
fractional  energy  deposition  is  so  small  that  the  atmospheric  rise  velocity  is  quite  small. 
Some  of  the  X-ray  energy  that  is  propagated  obliquely,  as  well  as  UV  energy  emitted  by 
the  burst,  are  absorbed  in  the  100-200  km  altitude  region  and  can  lead  to  substantial  heave 
(i.e.,  a  relatively  large  mass  moving  at  relatively  high  velocity  uh  ~  5-10  km/sec).  Also, 
that  part  of  the  kinetic  energy  of  the  debris  that  moves  down  (mainly  alcmg  the  geomagnetic 
field  lines),  produces  a  shock  in  the  atmosphere  below  ~  250  km  and  is  stopped  at  altitudes 
near  150-200  km,  giving  rise  to  a  large  energy  deposition  as  compared  to  the  ambient 
thermal  energy  and  thus  producing  substantial  heave  over  a  large  horizontal  area. 

Heave  is  particularly  important  for  high-altitude  multiple  burst  scenarios  where 
there  is  a  large  fractional  energy  deposition  at  high  altitudes  and  the  effects  of  heave  cannot 
necessarily  be  modeled  in  terms  of  single  detonations.  There  now  exist  some  computer 
codes  which  model  this  behavior.  Heave  is  illustrated  particularly  dramatically  by  the 
SCENARIO  calculations  of  Figs.  6.9b  and  6.9d,  which  show  how  density  contours  are 
raised  to  altitudes  very  much  above  ambient.  Consider,  for  instance,  contour  A,  which 
corresponds  to  the  ambient  density  at  477  km.  In  Fig.  6.9b  ( corresponding  to  the  plume 
of  60  Mt  detonations)  and  inunediately  above  the  detonations  (shown  by  black  points)  this 
contour  is  raised  to  3000  km;  at  other  locations  it  is  lifted  to  1500  km.  In  Fig.  6.9d  (which 
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shows  the  heaved  atmosphere  following  1000  Mt  explosions)  the  A-contour  is  totally  out 
of  the  picture  (i.e.,  above  300  km).  Again,  contour  G,  which  corresponds  to  96  km 
ambient  altitude,  is  lifted  above  500  km  in  the  60-Mt  plume  of  Fig.  6.9b  and  as  high  as 
1000  km  in  the  1000  Mt  example  of  Fig.  6.9d. 

Figure  6.10  shows  the  overall  character  of  heave,  by  demonstrating  the  ambient 
atmospheric  density  as  a  function  of  altitude  for  low  and  high  solar  activity.  This  makes 
quite  a  difference  at  altitudes  in  the  range  300-1500  km,  where  the  density  is 
10“^3  to  10~^0  times  that  at  sea  level.  The  lines  P(60)  and  H(IOOO)  show  representative 
values  of  density  as  a  function  of  height  for  the  plume  corresponding  to  the  60-burst  case 
of  Fig.  6.9c  and  to  the  total  heaved  atmosphere  arising  from  the  1000-burst  environment  of 
Fig.  6.9f.  In  general,  heave  occurs  for  air  above  150-200  km,  where  the  density  is  so  low 
that  the  energy  absorption  rate  is  large  enough  to  move  air  up  very  rapidly. 

Figure  6.10  demonstrates  different  regions  of  phenomenology;  at  densities  below 
10"^  kg/m3  (above  200  km  in  the  ambient  atmosphere)  there  is  heave,  while  at  densities 
below  10“13  kg/m^  (above  500-800  km  in  the  ambient  atmosphere)  the  phenomenology  is 
characterized  as  the  "Magnetic  Confmemeat  Regime"  (MCR).  The  circles  indicate  U.S. 
high-altitude  nuclear  explosions  for  which  there  is  a  significant  amount  of  data;  note  that 
only  two  lie  in  the  density  region  between  10”^^  and  10“*  kg/m^  in  which  heave  is 
significant  If  the  initial  bursts  of  a  multi-burst  scenario  be  in  the  MCR,  then  provided 
some  energy  absoiption  and  thus  heave  occurs  here  (which  is  most  probably  the  case)  the 
environment  of  greatest  interest  is  likely  to  be  in  the  density  range  of  10“^^  to  10-8  kg/m^. 
However,  we  have  no  test  data  in  the  MCR. 

Note  that  while  the  physics  of  heave  is  straightforward  ("hot  air  rises"),  there  is  no 
quantitative  experimental  evidence  for  it  There  is  qualitative  information  from  three  U.S. 
high-altitude  nuclear  tests,  and  apparently  also  some  evidence  from  solar  flares.  The 
available  data  for  the  highest-altitude  U.S.  nuclear  test  (SF)  have  not  provided  a  complete 
and  quantitative  accounting  of  where  the  debris  went.  We  know  that  the  atmospheric 
density  above  -  500  km  in  the  ambient  atmosphere  is  greater  at  2  p.m.  than  at  2  a.m.  (local 
time)  by  a  factor  2-5,  evidently  as  a  result  of  buoyancy  due  to  solar  heating.  However,  this 
is  clearly  an  example  of  buoyant  rise,  while  the  essential  features  of  heave  come  from  the 
rapid  absorption  of  energy  and  the  reflection  of  a  downward-moving  debris-induced  shock 
wave. 
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A  typical  rise  velocity  of  heaved  air  is  in  the  range  5-10  ktn/sec,  which  is  large 
compared  with  the  ambient  sound  velocity  but  small  conq)ared  to  the  hydromagnetic  Alfven 
velocity^* 

VA  =  B/(4jc  ni  mi)l^  ,  (6.17) 

(-100  km/s)  which  is  the  appropriate  atmospheric  velocity. 

How  long  does  the  heaved  air  stay  up?  It  moves  up  ballistically  under  gravity,  at 
speeds  uh  ~  5-10  km/s,  and  expands  horizontally  to  mix  with  the  ambient  very-low-density 
medium.  Eventually  it  falls  again  under  gravity.  The  fastest  possible  fall  rate  is  equal  to 
the  rise  rate,  i.e.,  the  fastest  fall  time 

tmin~UH/g~  8-16  min  ,  (6.18) 

where  g  =  9.8  m/s^  is  the  gravitational  acceleration. 

6.3.4  Damage  Mechanisms  due  to  High*Altitude  Nuclear  Bursts 

At  high  altitudes  where  there  are  no  fireballs,  the  damage  is  due  directly  to  the 
weapon's  radiative  output,  i.e..  X-rays,  yrays,  and  neutrons—see  Tables  6.1  and  6.2. 

There  is  one  additional  longer-time  damage  mechanism  arising  from  the  presence  of 
electrons  trapped  in  the  radiation  belts.  A  nuclear  explosion  in  the  upper  atmosphere 
generates  a  very  large  number  of  free  electrons,^^  and  at  sufficiently  high  altimdes  a 
significant  fraction  of  these  electrons  can  be  trapped  by  the  geomagnetic  field,  serving  to 
enhance  the  (Van  Allen)  radiation  belts  for  times  of  a  few  days  to  a  year.  The  effects  of  this 
can  be  particularly  serious  for  spacecraft  in  orbits  of  1000-2000  km  radius.  For  a 
discussion,  see  Cladis  et  al.,  1971,  and  Leong  et  al.,  1986. 

6.4  SYSTEMS  IMPLICATIONS  OF  HIGH  ALTITUDE  NUCLEAR 
EXPLOSIONS 

6.4.1  Introduction 

The  detonation  of  nuclear  weapons  in  the  uppier  atmosphere  can  affect  both  System 
Survivability  and  System  Operability.  System  survivability  is  a  measure  of  the  ability  of 

^  ^  In  a  plasfna  there  are  longitudinal  waves  whose  phase  velocity  is  the  normal  sound  speed.  However, 
there  are  also  purely  transverse  waves  whose  phase  velocity  is  the  Alfven  velocity-cf,  e.g.  Jackson, 
1975,  p.  486. 

1  Mt  roughly  doubles  the  total  number  of  electrons  in  the  global  ionosphere. 
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system  elements  (spacecraft,  air  or  ground  components)  to  survive  physically  even  after 
exposure  to  large  fluxes  of  X-rays,  y-rays  and  neutrons.  Tables  6.1  and  6.2  give  the 
exposure,  which  scales  as  yield/(range)2;  it  is  possible  to  harden  many  elements  to 
withstand  substantial  exposure  to  radiation.  System  Operability  is  more  subtle.  Some 
major  examples  are  the  following: 

•  Electronic  Upset 

When  an  electronic  component  encounters  a  radiation  level  significandy  less  (by 
several  orders  of  magnitude)  than  is  needed  to  damage  it,  there  may  still  be  a  false  signal 
which  affects  the  operation  of  the  component  and  of  a  device  made  up  of  such  components. 
Further,  after  a  system  has  been  exposed  to  a  high  level  of  radiation,  by  proper  engineering 
design  and  hardening  there  may  be  an  outage  but  the  system  can  function  again  (though 
possibly  at  reduced  capacity)  once  the  radiadon  has  reverted  to  a  normal  level. 

•  Electromagnetic  Pulse  (EMP) 

The  voltage  or  current  set  up  in  a  transmission  line  or  in  an  electronic  device  may 
interfere  with  the  oj-eration  of  electronic  devices,  in  particular  with  current-driven  low- 
voltage  devices  which  are  much  more  sensitive  to  EMP  than  are  voltage-driven  vacuum 
tubes.  Section  6.2.5  refers  to  more  detailed  calculations. 

•  Communication  Interference  ("Blackout") 

When  an  electromagnetic  signal  is  transmitted  through  the  nuclear-perturbed 
environment  which  has  not  only  a  significantly  enhanced  mean  degree  of  ionization  but  also 
large  fluctuations  in  electron  density,  its  propagation  may  be  degraded  significantly.  A 
simple  discussion  is  given  in  Section  6.4.2. 

•  Enhanced  IR  Atmospheric  Signature  ("Redout") 

An  infrared  sensor  looks  for  a  target  by  its  thermal  signature.  The  nuclear 
penurbed  environment  frequently  radiates  much  more  than  the  ambient  atmosphere^^,  and 
thus  the  target  may  not  be  detectable  at  the  ranges  that  are  expected.  Section  6.4.3  gives  a 
simple  discussion  of  Redout. 


^  ^  Note  that  the  IR  radiation  level  may  also  fluctuate  significantly,  producin''  "clutter"  which  can  interfere 
seriously  with  target  detection. 


6.4.2  Ionization  Effects  on  Radio  and  Radar  Propagation 


6. 4. 2.1  Electromagnetic  Wave  Propagation  in  an  Ionized  Medium 

At  high  altitudes  there  is  a  large  ionized  region  in  the  atmosphere  which  produces 
long-time  "blackout"  for  signals  of  fiequency  f  less  than  the  plasma  frequency 

fp  =  (ne  e2/7tm)l/2  =  9.0  x  IfP  x  n®*^  (Hz)  ,  (6.9) 

where  ne  =  no.  electrons/cm^;  if  ne  =  10^,  fp  =  90  MHz,  so  that  HF  and  VHP  are  severely 
affected.54  Note  that  even  at  much  higher  frequencies  there  will  be  a  variety  of  different 
propagation  effects.  Table  6.6  describes  many  of  the  more  important  effects  on  propa¬ 
gation  due  both  to  (enhanced)  ionization  and  also  due  to  small-scale  spatial/  temporal 
structure  ("striations"). 

6.4. 2.2  The  Disturbed  Ionosphere 

Figures  6.9c  and  6.9e  show  the  ionization  in  a  highly  nuclear  disturbed  (multi¬ 
burst)  environment  If  one  sketches  the  mean  ionization  as  a  function  of  altitude  under 
different  conditions,  one  obtains  the  schematic  results  of  Fig.  6.12,  which  shows  the 
following  as  a  function  of  altitude: 

•  Ambient  ionization  under  representative  mid-latitude  conditions 

•  1(1)  corresponds  to  the  enhanced  ionization  due  to  a  1-Mt  burst  at  ISO-km 
altitude  (from  Keskinen  and  Fedder,  1988,  from  Hain  et  al.,  1985). 

•  P(60)  goes  up  the  plume  shown  in  Fig.  6.9c  foUowing  the  designated  60-burst 
laydown. 

•  H(IOOO)  represents  the  total  heaved  environment  of  Fig.  6.9e  following  the 
designated  1000-burst  laydown. 

•  We  also  show  the  results  of  a  single  burst  MELT  code  run  (from  White  et  al, 
1987b)- 1  MT  at  200  km-respectively  at  30  minutes  and  3  hours  after  the 
burst,  which  should  be  compared  with  the  single-burst  curve  SI(1),  both  to 
indicate  possible  variabilities  of  different  predictions  and  to  remind  the  reader 
that  the  enhanced  high-altitude  ionization  disappears  very  slowly. 


The  nomtion  for  communication  frequency  bands  is  as  follows:  HF  is  3-30  MHz;  VHF  is 
30-300  MHz;  UHF  is  300-3000  MHz  (3  GHz);  SHF  is  3-30  GHz;  EHF  is  30-300  GHz.  For  radar. 
L-band  is  1-2  GHz;  S-band  is  2-4  GHz;  C-l^d  is  4-8  GHz;  X-band  is  8-I2.S  GHz;  Ku-band  is 
12.5-18  GHz;  K-band  is  18-26.5  GHz;  and  K,-band  is  26.5-40  GHz. 


Table  6.6.  Propagation  Disturbances  In  a  Structured  Ionized  Medium 
(from  Bogusch  and  Glugliano,  1982) 


Propagation 

Disturbance 

Effect  on 

Propaqated  Siqnal 

Effects  on 

System  Operation 

Absorption 

Loss  of  RF  energy  due  to  electron-ton 
and  electron-neutral  collisions  in 
ionized  regions  causing  signal 
attenuation. 

Reduction  in  SNR  below  demodulation 
threshold  causing  data  errors  and 
possible  loss  of  tracking  functions. 

Noise 

Random  RF  energy  emission  from 
hot,  ionized  fireball  regions. 

Increased  background  noise  levels 
causing  reduction  in  SNR,  increased 
data  errors  and  track  errors. 

Phase  and 
Doppler  shift 

Signal  phase  variation  due  to  TEC 
(=  Total  Electron  Content) 
changes  along  propagation  path. 

Possible  loss  of  phase/frequency  lock, 
and  difficulty  in  carrier  acquisition  and 
tracking. 

Time  delay 

Signal  time-of-arrival  variation  due  to 
changes  in  TEC  along  propagation 
path. 

Possible  toss  of  timing  lock  and 
difficulty  in  bit  synchronization  and 
tracking. 

Dispersion 

Signal  energy/time  dispersion  due  to 
spectral  components  propagating  at 
different  velocities  in  ionized  regions. 

Signal  waveform  distortion  causing 
degraded  demodulation  and  possibly 
inter-symbol  interference  at  high  data 
rates. 

Faraday 

rotation 

Rotation  of  signal  plane  of  polarization 
due  to  increased  magneto-ionic 
effects. 

Severe  amplitude  fading  with  linear 
polarization;  no  effect  with  circular 
polarization. 

Amplitude 

scintillation 

Rarxtom  fluctuations  in  signal 
amplitude  due  to  scattering  of  signal 
energy  by  ionization  irregularity 
structure. 

Severe  amplitude  fading  below 
demodulation  threshold  causing 
increased  data  error  rates  and  possible 
loss  of  track. 

Phase 

scintillation 

Random  fluctuations  of  signal  phase 
and  Doppler  due  to  TEC  variations  and 
scattering  of  signal  energy  by  ionization 
irregularity  structure. 

Possible  loss  of  phase  lock  and  loss  of 
signal  coherence  leading  to  possible 
catastrophic  derrKxlulatton  failure. 

Angular 

scattering 

Random  fluctuations  of  signal  angle-of- 
arrival  due  to  scattering  of  signal 
energy  by  ionization  irregularities. 

Reduction  of  SNR  due  to  scattering 
loss  in  large  antennas  causing 
increased  data  and  tracking  errors. 

Time  delay 
jitter 

Random  fluctuations  and  dispersion  of 
signal  time  delay  due  to  TEC  varurtton 
and  frequency  selective  scattering. 

Signal  waveform  distortion  and  timing 
jitter  causing  degraded  demodulation 
and  PN  tracking,  and  possibly  severe 
intersymbol  interference  at  high  data 
rates. 
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Figure  6.12.  Enhanced  Ionization  due  to  Single  and  Multiple  Nuclear 
Explosions  as  Function  of  Altitude 


Note  that: 


(i)  The  normal  ionization  shows  a  significant  amount  of  variability  in  both  mean 
level  and  in  both  spatial  and  temporal  structure. 

(ii)  A  single,  high-altitude  nuclear  burst  produces  several  orders  of  magnitude  in 
enhancement  of  mean  ionization  and  also  enhanced  structure. 

(iii)  In  a  multi-burst  situation  there  is  approximate  saturation,  with  maximum  mean 
ionization  in  the  range  of  lO^^-lO^*  electrons/cm^.  This  probably  occurs  for 
some  tens  of  (appropriately  located)  bursts:  the  physical  reason  is  that  the  first 
ionization  potential  of  an  air  atom  or  molecule  is  12-15  eV,  while  the  second 
ionization  potential  is  ~  35  eV  (cf  Part  I  of  these  Notes,  Chapter  2),  so  that  as 
the  air  is  heated  there  will  generally  be  no  more  than  one  free  electron  per  air 
atom  or  molecule.  Reference  to  Fig.  6. 10b  shows  that  on  account  of  heave  the 
peak  total  number  density  is  -  10^*  /cm^,  which  leads  to  a  comparable 
saturation  value  for  the  electron  density. 

6. 4. 2. 3  Some  Design  Parameters  for  EHF  Propagation  in  a  Structured  < 

Ionized  Medium 

To  provide  a  feeling  for  the  magnitude  of  propagation  effects,  in  Table  6.7  we 
show  some  propagation  parameters  for  a  1000-km  path  in  a  medium  of  mean  ionization 
10*0  electrons/cm^  55  and  also  with  with  large  fluctuations  of  RMS  value  Aue^  =  * 

2.5  X  10*^  (electrons/  cm3)2.  The  results  are  obtained  from  the  DNA/Wittwer  code 
SIMPRP,  version  2.  We  show  the  following  parameters,  as  functions  of  signal  frequency: 

•  Attenuation  A  (dB) 

•  Antenna  noise  temperature  Ta  (K) 

•  Signal  decorrelation  time  x©  which  corresponds  to  the  average  time  fluctuation 
in  signal  strength 

•  Frequency-selective  bandwidth  f©  which  measures  the  frequency  separation  at  { 

which  the  fading  histories  of  two  different  frequency  components  begin  to 

become  decorrelated  from  one  another. 


55  At  this  electron  density  the  plasma  frequency  fp  is  0.9  GHz. 
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Table  6.7A  gives  propagation  results  as  a  function  of  frequency,  while  Table  6.7B 
gives  more  results  from  SIMPRP  for  different  values  of  ionization  density  and  path  length 
at  a  single  frequency  (60  GHz)  to  indicate  the  nature  of  the  parametric  variabilities.^^ 

Table  6.7.  Electromagnetic  Propagation  Parameters  In  a 
Structured  Ionized  Medium 


A.  Frequency  Dependence  for  1 000-km  Path  with  ng  =  1 0^  ^  cm~3 


f(GH2) 

1 

10 

20 

40 

60 

A(dB) 

6.1E3 

49 

11.3 

2.6 

1.1 

Ta(K) 

10,000* 

10,000 

9300 

4500 

2200 

To  (sec) 

1 .9  E-5 

1 .9  E-4 

3.8  E-4 

7.6  E-4 

1.1  E-3 

fo(H2) 

1.1 

1.1  E4 

1.8  E5 

2.8  E6 

1.4  E7 

B.  Propagation  at  60  GHz  as  Path  Length  and  Electron  Density  are  Varied 


Electron  density  (cm“3) 

1010 

1010 

10’^ 

1011 

Layer  thickness  (km) 

500 

1000 

500 

10000 

A(dB) 

0.547 

1.09 

54.7 

109 

Ta(K) 

1200 

2220 

10,000 

10,000 

To  (sec) 

1 .6  E-3 

1.15  E-3 

1 .6  E-4 

1.15  E-4 

fo(Hz) 

2.7  E7 

1.4E7 

2.7  E5 

21.4  E5 

“  10,000  is  default  value  for  antenna  temperature. 

^  ^  The  code  can  vary  a  number  of  parameters  such  as  geomagnetic  field,  effective  antenna  size  and  shape, 
etc.  Here  we  use  the  "default"  values;  typically  there  can  be  a  variability  of  order  2  for  many  of  the 
parameter  choices. 
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6.4.3  Redout:  Backgrounds  for  IR  Sensors 

The  IR  signature  of  high-altitude  nuclear  explosions  has  been  studied  for  many 
years,  beginning  some  eight  years  after  the  last  U.S.  Atmospheric  Nuclear  tests,  as 
cryogenically  cooled  sohd-staie  LWIR  detectors  became  available.^^ 

The  presently  accepted  physics  is  intricate: 

•  In  the  normal  atmosphere 

(1)  Most  IR  radiance  is  due  to  the  molecular  minor  species  NO,  CO,  CO2, 
H2O,  and  O3  because  the  major  molecular  species  N2  and  O2  are  homonuclear 
diatomic  molecules  and  thus  have  no  allowed  IR  dipole  transitions. 

(2)  Peak  volume  mixing  ratios^®  of  these  molecules  in  the  20-60  km  altitude 
range  are : 


C02 

3x  1(H 

1  X  l(^5 

H20 

5  X  10-6 

N0-(-N02 

I  X  10-9 

CO 

"small" 

(3)  Triatomic  molecules  are  generally  photodissociated  by  solar  UV  radiation 
at  altitudes  above  80-100  km,  so  that  the  only  important  molecular  radiating 
species  is  NO,  which  is  formed  (to  a  few  percent)  when  air  is  heated  above 
20(X)-2500  K,  and  which  then  does  not  disappear  but  rather  "freezes-in"  as  the 
air  cools  (see,  e.g..  Part  I  of  these  Notes,  Section  2.5). 

•  In  the  nuclear-disturbed  atmosphere  (above  40-60  km,  where  non-equilibrium 
chemistry  is  important): 

(4)  At  very  early  times  when  the  temperature  is  high  and  the  electron  density 
is  very  large,  plasma  radiation^^  predominates. 


During  the  atmospheric  nuclear  test  series  there  were  hardly  any  measurements  with  sensitive  detectors 
beyond  about  3  pm,  because  the  present  generauon  of  IR  sensors  did  not  exist.  Further,  it  was  not 
clear  to  the  scientific  advisors  what  could  be  learned  from  IR  measurements.  [Jack  Carpenter 
(Visidyne),  private  communication.]  For  a  review  of  IR  data  during  the  tests  see  Reed  et  al.,  1980. 

^  ®  Defined  in  standard  way  as  the  ratio  of  the  volume  (not  mass)  of  CO2,  etc.,  to  the  volume  of  air. 

Electron-ion  and  elecpon-neutral  free-free,  free-bound,  and  bound-bound  radiation  (see,  e.g..  Part  I  of 
these  Notes,  Section  3.6). 
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(5)  This  is  particularly  true  at  high  altitudes  (>  150  km)  where  there  are  few 
radiating  air  molecules.  However,  current  calculations  suggest  that  there  is 
significant  radiation  above  150  km  due  to  CO2  and  possibly  other  molecules 
that  are  heaved  up,  as  well  as  due  to  NO  produced  in  air  that  is  heated  above 
2000  K.  Thus  the  assumption  that  plasma  radiation  predominates  over 
molecular  radiation  may  not  be  generally  correct. 

(6)  There  will  be  some  radiation  from  debris  species,  arising  from  the  nuclear 
weapon  and  its  carrier  vehicle,  but  this  is  not  thought  generally  to  be  very 
important  (except  possibly  in  particular  spectral  regions)  because  there  is  so 
very  little  of  this  material  present  in  comparison  to  the  heated  air. 

Figure  6. 1 3  gives  a  representative  model  for  the  IR  radiance  from  a  1  Mt  explosion 
at  200  km  as  viewed  by  a  limb-viewing  sensor  looking  through  the  atmosphere  at  85-km 
and  120-km  tangent  height  (reference  to  Fig.  6.7  indicates  that  at  85-km  tangent  height  one 
views  the  X-ray  and  beta-patches^®,  while  at  120  km  one  sees  mainly  the  debris  patch  and 
the  ionized  and  excited  air).  At  high  altitudes  and  early  times  (when  the  ionization  is  large) 
there  is  a  strong  continuum  due  to  plasma  radiation,  while  at  late  times  (when  the  ionization 
is  relatively  weak)  the  band  radiation  due  to  NO  (at  2.7  and  5.4  pm),  due  to  CO2  (at  2.8, 
4.3,  and  15  pm)  and  due  to  O3  (at  9.6  pm)  predominates. 

Note  that  at  85  km  there  is  relatively  little  nuclear  enhancement  (factor  10^  versus 
10^  at  140  km)  because  little  of  the  nuclear-induced  excitation  penetrates  down  as  low  as 
85  km. 

At  later  times  and  lower  altitudes  the  atmospheric  absorption  due  to  CO2  becomes 
important.  Note  that  the  time  variation  of  the  radiance  is  much  greater  at  120  km  than  at 
85  km.  I  show  these  figures  here  to  urge  the  reader  to  refer  to  current  Defense  Nuclear 
Agency  (DNA)  literature  for  the  most  current  model  predictions,  bearing  in  mind  that,  in 
the  absence  of  nuclear  test  data,  the  predictions  are  based  almost  entirely  on  computer 
modeling  and  analogy  with  auroral  and  other  simulation  experiments. 


Which  are  shown  in  Fig.  6.7  and  discussed  in  Section  6.3.I.3. 
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LOG  SPECTRAL  RADIANCE  (W/cm2/sr/micron) 


LiMg^VIEWlNG  RADIANCE 

TAN  ALT/NORSE  851217 


Figure  6.13.  Llmb<Viewlng  IR  Radiance  Following  a  1-Mt  Detonation  at  200-km 
HOB.  (a)  85-km  Tangent  Height,  (b)  I40*km  Tangent  Height. 

[Source:  NORSE  Code-Stephens  (1986)]. 
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LOG  spectral  radiance  (W/cm2/sr/micron) 


LIMB-VIEWING  RADIANCE 


Fig.  6.13  (Continued) 
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6.5  QUALITY  OF  THE  PREDICTIONS 

This  discussion  emphasizes  high-altitude  backgrounds  (ionization  and  infrared)  for 
multiburst  scenarios  because  this  is  the  area  of  greatest  uncertainty  and  importance  for 
Strategic  Defense  applications.  By  contrast,  for  low-altitude  nuclear  explosions  there  is  a 
relatively  large  amount  of  data  and  our  understanding  of  the  phenomenology  and  effects  is 
relatively  good  (see  ENW  and  EM-1  for  results).  Further,  since  fireballs  are  spatially 
confined  at  low  altitudes  because  of  the  relatively  high  atmospheric  density,  for  multiple 
bursts  one  simply  superposes  the  single  burst  effects  linearly  at  low  altitudes. 

At  intermediate  altitudes  (30-90  km),  non-equilibrium  chemistry  becomes  very 
important.  In  the  higher  regions,  where  there  is  not  a  well-confined  fireball  because  of  the 
low  I'mbient  atmospheric  density  there  is  the  possibility  of  multiburst  IR  signatures  larger 
than  at  higher  altitudes  (Carpenter,  1989).  Note  that  the  data  base  is  quite  poor;  there  have 
been  just  three  U.S.  atmospheric  nuclear  explosions  (indicated  on  Fig.  6.10)  between  30- 
and  90-km  burst  height. 

At  high  altitudes  (>  90  km)  the  chemistry  becomes  simpler  because  there  are 
relatively  few  chemically  or  radiatively  interesting  molecules,®^  and  again  there  is  little 
data-Just  three  well-instrumented  and  analyzed  U.S.  bursts  (which  arc  indicated  in  Fig. 
6.10),  At  these  high  altitudes  and  correspondingly  very  low  densities  the  bomb  energy  is 
distributed  over  very  large  regions  of  space,  multiburst  interactions  will  be  important-and 
there  are  simply  no  data. 

There  is  a  particular  problem  in  the  infrared;  while  there  are  extensive  radio/radar 
and  visible  optical  data,  and  also  data  in  the  short-wave  IR,  there  are  essentially  no 
observations  beyond  3-5  pm^^  Thus,  the  predicted  IR  signatures  come  from  computer 
simulations  which  are  based  on  visible  diagnostics  during  the  1962  and  1958  atmospheric 
test  series. 

Table  6.8  provides  a  tentative  listing  of  the  level  of  understanding  of  nuclear 
environment  parameters,  going  from  those  instances  in  which  the  understanding  is  good 


^  ^  All  molecules  other  than  N2  are  essentially  dissociated  above  100  km;  some  1-3%  of  nitric  oxide  (NO) 
is  made  when  air  is  heated  above  2S00  K,  and  this  tends  to  remain  ("frozen-in”)  as  the  air  cools.  In 
addition,  some  air  can  be  heaved  up  from  below  100  km,  so  that  there  is  the  possibility  of  small 
amounts  of  the  DR  radiating  species  CO2  and  O3,  as  well  as  of  CO  produced  by  the  dissociation  of 
CO2  (also  smaller  quantities  of  H2O  and  OH  from  the  1-5  parts  per  million  of  H2O). 

See  footnote  51. 
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(e.g.,  survivability,  in  particular  X-ray  hardness)  to  those  cases  where  the  understanding  is 
poor  (e.g.,  UVA^IS  sensor  operability 

Table  6.8.  Levels  of  Understanding  of  Nuclear  Environment 


State  of  Knowledge 


Examples 


Effects  of  Uncertainty 


Good 

Survivability: 

X-ray  hardness 

P&Y  Hardness 

Fair/Moderate 

Communications: 

ionization 

Moderate 

Structure 

Can  be  large 

Moderate/Poor 

Sensor  Operation: 

Radar 

See  Comm 

IR 

Mod/Large 

UVA/IS 

Large? 

Structure 

Can  be  large 

Poor/Zero 

(We  could  be  missing  some  physics  because  there  are  no  multiburst 

high-attitude  data) 

In  conclusion,  the  system  architect  and  designer  dealing  with  the  nuclear-perturbed 
environment  should  recognize  that  he  is  faced  with  three  independent  sets  of  uncertainties, 
all  of  which  may  be  quite  large: 

(1)  Uncertainties  in  weapon  output  for  specified  yield. 

(2)  Uncertainties  in  threat  scenarios. 

(3)  Uncertainties  in  the  nuclear  effects  codes. 


^3  For  a  more  extended  discussion  of  these  uncertainties,  see  Bauer.  1990b. 
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7.0  SOME  TOPICS  IN  THE  PROPAGATION  OF 
ELECTROMAGNETIC  RADIATION 
THROUGH  THE  ATMOSPHERE 

7.1  INTRODUCTION 

The  propagation  of  electromagnetic  radiation  through  the  atmosphere  is  important 
for  a  variety  of  applications.  Here  we  present  the  physics  for  a  number  of  different 
situations.  We  begin  by  discussing  the  effects  of  molecular  absorption,  which  is  very 
strongly  wavelength-dependent  with  transmission  windows  alternating  with  absorption 
bands.  Figure  7.1  demonstrates  the  atmospheric  transmittance  through  a  1  km  (horizontal) 
path  at  heights  HI  =  0,  10  km,  20  km,  and  30  km,  at  wavelengths  out  to  5  ^lm.  Note  the 
strong  wavelength  dependence  and  also  how  much  better  the  transmittance  becomes  at  the 
higher  altitudes  where  the  atmospheric  density  is  lower,  roughly  by  a  factor  4  for  each 
10  km  variation  in  altitude.  The  results  shown  here  come  from  the  USAF  LOWTRAN  and 
FAS  CODE  computer  models,  which  are  appropriate,  respectively,  for  low  resolution 
(20  cm~i)  applications,  and  high  resolution  (<  1  cm”D,  which  is  equivalent  to  a  single 
spectral  line  and  is  thus  appropriate  for  laser  propagation  and  other  applications  calling  for 
very  high  spectral  resolution.  Section  7.2  sketches  the  physics  contained  in  the 
LOWTRAN  and  FASCODE  atmospheric  transmission  codes—not  how  to  run  them. 
Section  7.3  discusses  some  characteristics  of  atmospheric  aerosols,  which  are  particles  of 
dust  and  chemical  pollutants  as  well  as  water  or  ice  droplets.  This  section  includes  a 
discussion  of  the  basics  of  atmospheric  visibility  and  of  Mie  scattering  calculations. 

Clouds  occur  very  frequently;  optically  thick  clouds  which  block  the  transmission 
of  visible  or  infrared  radiation  occur  some  20-80  percent  of  the  time,  depending  on 
location,  altitude,  season,  and  time  of  day.  Table  7. 1  gives  mean  values  for  high  and  total 
cloudiness  for  January  and  July  at  selected  locations.  Section  7.4  gives  a  very  simple 
discussion  of  the  physics  of  water/ice  particle  clouds.  In  Section  7.5  we  discuss  a  very 
specific  application,  the  analysis  of  Cloud-Free  Lines  of  Sight  (CFLOS). 

For  the  propagation  of  High-Energy  Laser  Beams  (or  of  other  high-energy  and/or 
high-resolution  light  beams)  there  is  an  additional  series  of  problems.  Section  7.6  treats 
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key  characteristics  of  atmospheric  turbulence  which  serves  to  broaden  narrow  beams,  while 
Section  7.7  discusses  some  nonlinear  processes  related  to  the  propagation  of  high-energy 
laser  radiation,  in  particular  thermal  blooming  and  Stimulated  Raman  Scattering  (SRS). 
The  object  of  the  present  discussion  is  to  present  a  very  simple  introduction  to  all  these 
topics  to  enable  the  reader  to  go  funher  on  his  own,  starting  with  the  reference  material 
cited  here. 


Figure  7.1.  Atmospheric  Transmittance  Models;  LOWTRAN  and  HITRAN 
(now  FASCODE)  Compared  (Source:  R.W.  Fenn,  AFGL, 
private  communication,  1983) 
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Table  7.1.  High  and  Total  Cloudiness  at  Representative  Locations  In  the 
Northern  Hemisphere  (Source:  Malick  and  Allen,  1978,  1979) 


Location 

Coordinates 

High/Total 

Cloudiness 

Longitude 

Latitude 

January 

July 

China  Lake,  CA 

36“N 

117“W 

.17/. 38 

.12/. 18 

Grand  Forks,  NO 

48“N 

95“W 

.39/. 63 

.31/. 56 

Maui ,  HI 

21“N 

156“W 

.14/. 40 

.12/. 50 

Hudson  Bay 

60“N 

88“W 

.06/. 36 

.08/. 29 

N.  Atlantic  S 

52“N 

35“W 

.24/. 81 

.13/. 70 

N.  Atlantic  N 

62“N 

30“W 

.18/. 76 

.16/. 72 

Jan  Mayen  Is 

71“N 

10°W 

.20/. 81 

.16/. 85 

Thule 

76'’N 

68“W 

.10/. 35 

.11/. 73 

Barrow,  AK 

7rN 

ise^w 

.08/. 34 

.11/. 64 

Arabian  Sea 

8°N 

65“E 

.02/. 23 

.16/. 55 

Teheran 

36“N 

52°E 

.14/. 38 

.02/. 22 

Ionian  Sea 

39"N 

18»E 

.07/. 54 

.01/. 06 

Moscow 

56“N 

39'E 

.22/. 61 

.24/. 46 

Tyuratam 

46®N 

64 '’E 

.16/. 49 

.13/. 30 

Lop  Nor 

40''N 

9PE 

.22/. 48 

.24/. 57 

Vladivostok 

43“N 

132‘’E 

.11/. 43 

.22/. 66 

Japanese  Trough 

SS^N 

ISO^E 

.16/. 67 

.12/. 37 

Anadyr 

64'’N 

177°E 

.28/. 59 

.21/. 75 

Murmansk 

69“N 

34°E 

.22/. 70 

.16/. 66 

Notes: 

a.  High/total  cloudiness  means,  for  example,  that  at  China  Lake  in 
January  high  clouds  occur  0.17  of  the  time  and  total  cloudiness 
occurs  0.38  of  the  time. 

b.  These  data  come  largely  from  downward  viewing  satellites  such 
as  NOAA-6  and  OMSP,  which  tend  to  underreport  optically  thin 
ciouos. 

c.  Clouds  are  reported  as  present  when  at  least  1/10  of  the 
appropriate  field  of  view  is  covered  by  clouds. 

d.  High  clouds  are  those  above  7  km,  with  the  altitude  determined 
by  the  effective  radiative  temperature  in  the  10  to  12  urn 
infrared  band  as  compared  with  the  atmospheric  temperature/ 
altitude  profile. 

e.  High  clouds  are  thus  mainly  moderately  thick  cirrus  or  cirro- 
stratus,  plus  some  cumulonimbus  (thunderclouds)  at  the  lower 
latitudes  (<30*). 
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7.2  TRANSMISSION  THROUGH  THE  ATMOSPHERE:  MOLECULAR 

ABSORPTION 

7.2.1  Introduction 

Transmission  through  the  atmosphere  is  good  in  those  spectral  regions  in  which  the 
molecular  and  atomic  species  present  (which  are  listed,  for  instance.  Table  2-1  of  Part  I)  do 
not  absorb  (or  emit,  by  Kirchoff s  Law),  the  so-called  "window  regions",  and  bad  where 
atmospheric  species  absorb.  Atoms  have  electronic  transitions,  generally  with  excitation 
energies  somewhat  below  the  first  ionization  potential,  which  is  generally  in  the  range  of 
10-15  eV— see,  e.g..  Table  2-3  of  Part  Molecules  have  electronic  structure  somewhat 
different  from  that  of  the  atoms  that  make  them  up,  but  in  addition  they  also  have  rotational 
and  vibrational  structure,^^  and  they  have  dissociation  energies  in  the  range  1-10  eV— see, 
e.g..  Table  2-3  of  Part  I. 

Table  7-2  characterizes  the  different  spectral  ranges  in  terms  of  the  energy  of  the 
transitions  and  in  terms  of  the  dominant  transitions  that  are  responsible  for  absorption  in 
them.  We  specify  wavelengths  in  nm  (1  nm  =  10~^  m),  and  fim  (1  iim  =  10“^  m),^^ 
and  energies  in  eV;  for  quantitative  work  it  is  customary  to  use  energies  in  cm"* 
(1  eV  =  8067  cm"*)  instead  of  these  two  sets  of  units.  The  labeling  of  the  spectral  ranges 
is  generic  and  not  precise: 

•  The  eye  is  sensitive  to  wavelengths  in  the  Visible  spectral  range,  400-700  nm, 
and  (not  coincidentally)  much  of  the  solar  irradiance  is  in  this  spectral  range. 

•  Shorter  wavelengths  are  called  Ultraviolet  (UV),  which  is  sometimes  divided 
into  several  regions; 

•  EUV,  X  <  100  nm,  is  actually  not  of  great  interest  at  present. 


But  note  that  in  general  there  are  also  much  higher  excited  states,  since  an  atom  such  as  O  with  eight 
elecuons  has  an  ionization  potential  of  0.53  keV  for  the  most  suongly  bound  electron,  as  compared 
with  the  first  ionization  potential  of  13.6  eV,  and  furthermore  multiplets  such  as  the  0(^P)  ground 
state  has  transitions  of  much  lower  energy  ( -  .02  eV)  between  the  three  levels  of  the  triplet. 

For  a  rotational  quantum  number  J  (defined  in  Part  I,  especially  Appendices  A  and  B),  the  energy  is 
hcBeJ(J+l);  for  O2,  J  =  1,  this  is  2.89  cm"*  or  .00036  eV.  For  vibrational  quantum  number  v,  the 
vibrational  energy  is  hc(0e(''+l/2);  for  O2,  v=l  this  is  2370  cm"*  or  0.29  eV. 

Note  that  wavelengths  are  often  measured  in  A  =  10"*  cm  =  0.1  nm-  see,  e.g..  Fig.  3-6,  p.  3-31. 


•  Far  UV,  100  nm  <  X  <  200  nm  includes  some  wavelengths  where  the 
atmosphere  absorbs  very  strongly;  the  O2  molecule  has  the  ver>’  strong 
"Schumann-Runge"  absorption  in  the  130-180  nm  spectral  range. 

•  Mid  UV  (MUV)  is  sometimes  defined  as  200  nm  <  A,  <  300  nm;  ozone  O3  is  a 
very  strong  absorber  in  the  230  -  270  nm  spectral  range. 

•  Near  UV  (NUV)  may  be  defined  as  300  nm  <  A  <  400  nm. 

•  The  division  of  the  infrared  (IR)  is  sometimes  made  into  near  IR  or 
photographic  IR;  some  special  photographic  films  are  sensitive  out  to  ~  1000 
nm  or  1  }xm.  The  SWIR  ("short  wavelength  IR")  corresponds  to  wavelengths 
less  than  about  3  pm,  the  MWIR  ("mid-wave  IR")  is  traditionally  taken  as  3 
pm  <  1  <  8  pm,  while  the  LWIR  ("long  wavelength  IR")  corresponds  to 
wavelengths  above  8-10  pm. 

•  Beyond  A  =  100  pm  =  0.1  mm  we  enter  the  Millimeter  Wave  Region  which  is 
not  discussed  here-see,  e.g.,  Kulpa  and  Brown,  1979. 


Table  7-2.  Listing  of  Different  Spectral  Regions 


Name 

Wavelength 

range{nm) 

Energy  range 
hc/l(eV) 

Absorption 
mainly  due  to: 

EUV 

<  100 

>  12.4 

FUV 

100-200 

12.4-6.2 

electronic  transitions  (atomic) 

MUV 

200-300 

6.2-4.13 

electronic  transitions  (molecular) 

NUV 

300-400 

4.13-3.1 

electronic  transitions  (molecular) 

VIS 

400-700 

3. 1-1.8 

a  few  molecular  electronic  transitions 

NIR 

700-1000 

1.8-1.24 

a  few  molecular  electronic  transitions 

SWIR 

1-3  pm 

1.24-0.41 

molecular  vibrational  transitions 

MWIR 

3-8  pm 

0.41-0.155 

molecular  vibrational  transitions 

LWIR 

8-25+  pm 

<0.155  eV 

vibrational  transitions:  H2O  rotation 

7.2.2  The  Nature  of  the  Problem 

Some  excellent  photographs  of  molecular  electronic  band  spectra  are  given  in 
Herzberg,  Spectra  of  Diatomic  Molecules,  1950,  Ch.  II.  They  show  a  great  richness  of 
phenomena,  with  continua  as  well  as  both  coarse  and  fine  spectral  structure.  The  coarse. 
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banded  structure  (-  10^  cm“'  to  10^  cm~^  in  total  energy  width)  is  due  to  molecular 
vibrations,  while  the  fine  structure  due  to  molecular  rotational  lines  ha«  lines  spaced 
10~^-102  cm“l  apart  In  general,  the  vibrational  structure  has  sharp  edges  on  one  side  with 
a  gradual  falloff  on  the  back  side  (either  toward  higher  or  lower  frequency).  A  typical 
electronic  band  system  covers  perhaps  10,()(X)  cm~^  to  20,000  cm"'  (~  100  nm  in  visible), 
so  that  it  may  have  ~  10^  individual  lines,  of  varying  spacing  and  strength.  From  the 
discussion  of  Section  3.4  of  Part  I  of  these  Notes,  we  see  that  in  the  lower  atmosphere 
(below  10  to  40  km)  most  spectral  lines  are  collision  broadened  with  a  typical  line  width 
-0.1  cm"',  although  at  higher  altitudes  Doppler  broadening  predominates  (see  Goody  and 
Yung,  1989,  p.  97f). 

In  the  IR  there  are  molecular  vibration-rotation  bands  due  to  molecules  such  as  NO, 
CO,  CO2,  h20,  O3,  etc.  The  radiating  species  are  typically  optically  active  molecules 
present  in  much  smaller  concentrations  than  the  major  species  N2  and  O2  which  have  such 
high  molecular  symmetry  that  they  are  not  important  radiators.  Which  species  are 
important  radiators  depends  strongly  on  the  spectral  region  and  temperature  range. 

The  applications  considered  determine  how  the  band  structure  should  be  treated.  If 
one  asks  for  the  emission  from  high-temperature  gases  at  low  spectral  resolution,  very  little 
detail  is  necessary.  In  particular,  for  an  optically  thin  situation  in  which  self-absorption  and 
scattering  are  unimportant,  the  total  integrated  intensity  is  all  that  matters,  and  thus  a 
Planck-type  mean  [see  Eq.  (3.23)  of  Section  3.3.3  of  Part  I]  is  adequate.  Thus  the 
description  in  terms  of  electronic  system  f-numbers  with  vibrational  structure  analyzed  in 
terms  of  a  Franck-Condon  factor  is  satisfactory  for  determining  the  total  radiation  from  air 
or  combustion  products  at  high  temperatures. 

However,  if  one  asks  for  the  atmospheric  transminance  through  1-  to  10-km  path 
lengths  (or  through  the  whole  atmosphere  in  vertical  viewing),  especially  at  high  spectral 
resolution  in  the  IR,  where  there  are  numerous  absorption  regions  alternating  with 
"windows,"  much  more  detail  is  needed.  In  particular,  when  one  views  not  a  gray^^  body 
through  a  long  atmospheric  path  but  rather  a  rocket  plume  or  other  array  of  hot  CO,  CO2, 
and  H2O  molecules  viewed  through  the  atmosphere  containing  cold  H2O  and  CO2 
molecules  (see,  for  instance.  Section  4.5. 1  and  Fig.  4.6  in  Ch.  FV  of  Part  I),  or  for  laser 
transmission  at  very  high  spectral  resolution,  very  much  detail  is  needed. 


A  grey  body  is  defined  as  one  having  an  emissivity  e  different  from  one,  but  not  varying  with 
wavelength  Thus  a  black  body  is  the  limit  of  a  grey  body  having  e  =  1. 
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During  the  1950s  and  1960s  a  variety  of  very  sensitive  solid-state  detectors  of  IR 
radiation  as  developed,  and  thus  many  specific  questions  regarding  atmospheric 
transmission  arose;  and  with  the  increase  in  computer  power  during  this  time  frame  a 
number  of  Band  Models  were  devised  and  tested  against  experiment.  These  models  tend 
to  assume  spectral  lines  of  either  equal  or  random  spacing,  of  uniform  or  varying  strength, 
sometimes  with  a  configuration  like  the  "Q-Branch"  of  a  band  (see,  e.g.,  Herzberg,  1950, 
p.  48).  Two  limits  are  the  weak-  and  strong-line  approximations.  In  the  weak-line  limit 
the  line  is  optically  thin  even  at  the  line  center,  so  that  the  emission  (or  absorption)  is 
proportional  to  the  number  of  molecules  in  the  path  [Aav  ~  (nfd)^*].  By  contrast,  in  the 
strong-line  limit  the  line  is  optically  thick  at  the  center,  so  that  increasing  the  number  of 
molecules  in  the  path  simply  increases  the  emission  (or  absorption)  in  the  line  wings,  and 
of  course  this  is  different,  depending  on  the  line  shape  (Doppler,  Lorentz,  or  Voigt-see, 
e.g..  Fig.  3.5);  the  strong-line  approximation  traditionally  corresponds  to  the  Lorentz-line 
limit  Aav  ~  (nfd)*^^.^^  For  a  review  of  band  models,  see,  e.g..  Goody  &  Yung,  1989, 
Ludwig  et  al.,  1973. 

During  the  1970s  the  LOWTRAN  transmittance  computer  code  was  developed  at 
the  U.S.  Air  Force  Geophysical  Laboratory  (AFGL)  for  automating  low-resolution  IR 
transmittance  calculations.^®  It  was  originally  meant  to  apply  to  the  lower  atmosphere,  as 
an  empirical  model  of  transmittance  through  perhaps  a  1-  to  10-km  horizontal  atmospheric 
path  at  altitudes  up  to  10-15  km,  with  Lorentz-broadened  atomic  spectral  lines.^^  As  more 
sophisticated  and  detailed  questions  were  asked  simultaneously  with  a  continued  rapid 
growth  in  available  computing  power,  line-by-line  r  aculations  became  feasible  and  so  the 
FASCODE  computer  program  was  developed  at  AFGL,  based  on  an  Atlas  of  Spectral 
Lines  (the  latest  version  of  which  is  referenced  by  Rothman  et  al.,  1987).  By  now  this 
Atlas  contains  some  300,000  individual  lines,  and  in  the  IR  line-by-line  calculations  can  be 


Aav  is  the  absorption  coefficient  averaged  over  the  relevant  spccUal  interval,  n  is  the  number  of 
absorbing  molecules  per  unit  volume,  f  is  the  oscillator  strength  of  the  relevant  transition,  and  d  is  the 
geometrical  path  length  under  consideration.  See  Eq,  (3.38a)  or  branch  QA  of  Fig.  3.5. 

cf.  Eq.  (3.38c),  or  branch  QC  in  Fig.  3.5. 

This  code,  which  uses  a  very  simple  empirical  band-type  model,  was  first  developed  as  a  graphical 
procedure— see  McClatchey  et  al.,  1970,  or  IR  Handbook,  1985,  p.  5-46ff,  which  was  then  automated 
for  digital  computation. 

^ '  Thus,  at  high  altitudes,  where  Voigt  or  Doppler  broadening  applies  for  most  specual  lines  rather  than 
Lorenu  broadening,  one  has  to  check  to  be  sure  that  the  LOWTRAN  results  are  correct  for  the 
particular  problem  under  consideration-see  Section  7.2.5. 

I 
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carried  through  with  the  FAS  CODE  computer  program.  This  could  not  have  been 
considered  20  years  ago. 

The  point  is  the  following.  In  an  optically  thin  medium,  if  there  are  n  absorbers  per 
unit  volume,  each  having  absorption  cross  section  Oabs.  then  for  a  path  length  d  one  can 
define  a  dimensionless  optical  thickness 

X=naabsd  (7.1) 

and  the  transmittance  x  would  be 

t=e-X  =e-noabsd  .  (72) 

In  the  context  of  atmospheric  transmittance,  Eq.(7.2)  is  known  as  "Beer's  Law". 

However,  for  a  high-temperature  molecular  gas  such  as  air,  the  absorption 
coefficient  as  averaged  over  a  "reasonable"  spectral  interval  does  not  follow  a  simple 
exponential  law,  so  that  Beer's  Law  of  Eq.(7.2)  is  generally  not  obeyed.  This  arises 
because  in  a  designated  spectral  range^^  the  spectral  absorption  varies  drastically  with 
frequency,  so  that  the  signal  saturates  in  some  portions  of  the  spectral  range  and  not  in 
others.  It  is  here  that  the  LOWTRAN  and  FASCODE  computer  codes  (or  something 
equivalent)  must  be  used. 

It  must  be  stressed  that  while  the  LOWTRAN  and  FASCODE  programs  enable  a 
wide  variety  of  atmospheric  transmittance  and  radiance  calculations  to  be  carried  out,  they 
are  not  appropriate  for  "all  conceivable"  calculations.  Thus  for  a  sophisticated  application  it 
is  important  to  check  to  be  sure  that  the  computer  program  used  is  appropriate.^^ 

7.2.3  The  LOWTRAN  Code 

This  section  should  be  read  after  Sections  3.4  and  3.5  of  Part  I  of  these  Notes, 
where  the  concepts  and  terms  used  here  are  explained.  This  material  is  presented  here 
(from  McClatchey  et  al.,  1970,  and  Kneizys  et  al.,  1980)  to  explain  to  the  user  what  is 
actually  in  the  code.  This  code  works  with  Av  =  20  cm~*  spectral  resolution  (although  a 


^  ^  Such  as  may  be  characterized  by  a  particular  filter. 

One  current  example  is  the  following.  For  the  earth's  surface-atmosphere  radiation  budget,  molecular 
absorption  in  the  LWIR  contributes  significantly.  Kralz  and  Cess,  1988,  have  compared  line-by-line, 
narrow-band  and  broad-band  models  for  the  9.6  and  14  pm  bands  of  ozone,  and  find  that  the  Malkmus 
narrow-band  model  can  provide  a  good  fit  to  a  line-by-line  description,  which  is  not  the  case  for 
methane. 


5  cm“^  version  exists).  It  is  basically  a  transmittance  code,  although  the  radiance  can  also 
be  provided. 

(1)  For  a  path  length  AS  one  defines  an  average  transmittance 

x(av,  Av;v)  =  (1/Av)  i  x(v)  Sv  (7.3) 

and  in  general—because  Beer's  Law  does  not  hold-  there  exists  no  function 
K(v)  such  that 

x(av,Av;v)  =  exp[-K(v)  AS]  (7.4) 

because  of  the  various  saturation  effects. 

(2)  The  average  transmittance  factor  x(av,  Av;v)  for  a  given  path  length  AS  (which 
will  now  be  written  just  as  x)  is  written  as  the  product  of  a  number  of 
components: 

X  =  x(Molecular).x(Continuum).x(Rayleigh).x(Aerosol)  .  (7.5) 

x(Aerosol)  is  discussed  in  Section  7.3.2;  x(R.*yieigh)  (scattering  by 
atmospheric  molecules  at  short  wavelengths)  is  treated  in  Section  7.3.4; 
x(Continuum)  (due  principally  to  N2  and  H2O  molecules)  is  discussed  in  the 
LOWTRAN  documents  (e.g.,  Kneizys  et  al.,  1980),  while  x(Molecular)  is 
broken  up  into  contributions  from  H2O,  O3  and  UMG  =  Uniformly  Mixed 
Gases  (CO2,  N2O,  CH4,  CO,  02,etc). 

x(Molecular)  =  x(H20).x(03).t(UMG)  .  (7.6) 

(3)  To  evaluate  x(Molecular)  one  defines  an  Equivalent  Optical  Depth 

EOD  =  Cv.(0*.AS  ,  (7.7) 

where  Cy  is  a  spectral  absorption  coefficient  (an  empirical  expression  is  used 
in  the  code,  see,  e.g.,  //?  Handbook,  1985,  Fig.  5-14,  p.  5-56),  AS  is  the  path 
length,  and  to*  is  an  equivalent  density  along  the  path,  characterized  at  altitude 
z  by  pressure  P(z)  and  temperature  T(z). 

CO*  =  [  (0  (P(z)/Po)  (T(/r(z))l/2]n  .  (7.8) 

where  co  =  concentration  of  the  absorber  and  Pq,  T©  refer  to  STP.  The 
exponent  n  has  the  value  n  =  0  for  the  weak-line  limit  and  n  =  1  for  the  strong¬ 
line  limit  (which  have  been  explained  in  the  discussion  of  band  mcxlels  above 
(Section  7.2.1)--see  also  Goody,  1964,  p.  1270-  The  LOWTRAN  Ccxle  uses 
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n  =  0.9  for  H2O 


0.75  for  UMG  • 

0.4  for  O3  , 

and  the  function  x(EOD)  as  used  in  LOWTRAN  is  shown  in  Fig.  7.2.  Note 

that  it  is  very  different  from  the  simple  Beer's  Law  form  exp  (-EOD)  because 

of  the  effects  of  spectral  saturation.  In  particular,  at  large  optical  thickness  the  • 

holes  in  absorption  coefficient  between  the  individual  spectral  lines  are  still 

frequently  optically  thin  so  that  the  transmittance  decreases  much  less  rapidly 

with  EOD  than  does  the  simple  form  exp  (-EOD).  The  latest  version  of 

LOWTRAN  (7,  Kneizys  et  al.,  1988)  extends  out  to  the  UV.  The  electronic 

band  systems  that  arc  important  in  the  visible  and  UV  range  arc  treated  in  terms  ® 

of  band  models,  like  the  vibrational  band  systems. 


Figure  7.2.  LOWTRAN  Empirical  Transmittance  Function  versus  EOD 
[Cv  (0*  AS  of  Eq.  (7.S),  from  Kneizys  et  al.,  1980], 
as  compared  with  the  simple  form  exp  (-EOD) 
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7.2.4  FASCODE,  a  Line-by-Line  Model 

This  code  uses  the  same  methodology  for  species  concentrations  and  path  geometry 
as  does  LOWTRAN,  but  the  spectral  information  is  treated  differently,  using  the  AFGL 
HITRAN  Atlas  of  Spectral  Lines  (Rothman  et  al.,  1987)  which  as  of  December  1987 
contains  some  300,000  lines  of  28  different  molecular  species  for  wave  numbers 
V  <17,900  cm“*  (X  >  0.57  iim).^"*  Where  lines  are  available,  the  band  systems  are  treated 
on  a  "Line-By-Line"  basis  using  the  full  Voigt  line  shape.  This  applies  mainly  to  the 
vibration-rotation  band  systems  in  the  IR. 

The  current  version  is  FASCODE  2,  to  be  replaced  in  Summer  1990  by 
FASCODE-3;  this  goes  out  to  v  <  50,000  cm“*  (0.2  pm),  so  that  it  contains  also  some 
electronic  band  systems  which  are  treated  by  using  band  models  the  same  as 
LOWTRAN.75 

At  present  LOWTRAN  and  FASCODE  are  coordinated  so  that  species 
concentrations,  geometry,  aerosol  models,  continua,  etc.,  are  all  treated  similarly;  the  only 
difference  is  the  treatment  of  molecular  band  systems  for  which  lines  exist  in  the  AFGL  line 
atlas.  Figure  7. 1  illustrates  some  different  applications  of  these  two  codes.^^  LOWTRAN 
models  apply  to  ambient  atmospheric  temperatures  only,  but  FASCODE  can  in  principle^^ 
treat  such  problems  as  the  propagation  of  spectrally  very  narrow  signals  such  as  laser 
radiation  through  a  designated  region  of  the  atmosphere,  or  the  propagation  of  radiation 
from  a  rocket  plume  [which  is  largely  hot  (2000  K)  H2O  and  CO2]  through  a  long 
atmospheric  path  containing  cold  (250  K)  air,  including  H2O  and  CC)2. 

7.2.5  Atmospheric  Transmittance/Radiance  Above  50  km  Altitude 

The  LOWTR  AN/FAS  CODE  series  of  models  were  originally  designed  for  use 
below  ~  10  km  altitude  where  the  amiosphere  is  in  local  thermodynamic  equilibrium  and 


A  1990  HITRAN  Tape  is  in  course  of  preparation. 

See  Clough  et  al.,  1989;  AFGL  Users  Manual  for  FASCODE-3  in  preparation. 

Both  codes  can  now  be  run  on  a  PC-AT/286  or  equivalent  microcomputer;  there  is  extensive 
documentation  on  how  to  run  the  code.  Indeed,  there  are  commercial  versions  for  both  LOWTRAN  and 
FASCODE  which  can  be  run  on  personal  computers. 

But  note  that  the  HITRAN  atlas  is  intended  for  atmospheric  applications,  i.e.,  room  temperature  or 
below-say  600  K  maximum.  Accordingly,  a  high-temperature  (>  1000  K)  gaseous  source  may  have 
substantial  contributions  from  spectral  lines  which  are  not  contained  in  the  HITRAN  line  atlas.  Thus, 
FASCODE  results  may  look  reasonable  but  be  incorrect  A  new  HITEMP  line  atlas  is  being  prepared 
to  supplement  this  deficiency  for  particular  molecules  in  specific  absorption  bands.  However,  a  very 
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where  the  assumption  of  Lorentz  or  collisional  spectral  line  shape  is  generally  good.  At 
altitudes  above  20-30  km  both  these  assumptions  tend  to  break  down:  as  the  density 
decreases,  the  spectral  lines  tend  to  be  Doppler  rather  than  Lorentz  broadened^*  and  the 
minor  chemical  species  which  are  important  optical  radiators/absorbers  are  frequently  not  in 
Local  Thermodynamic  Equilibrium,  so  that  they  may  be  described  as  having  different 
values  for  translational,  rotational,  vibrational  and  excitational/dissociative  effective 
temperatures. 

Thus  at  altitudes  above  30-50  km  one  should  consider  both  Doppler  line  shapes  and 
the  effects  of  NLTE  (Non-Local  Thermodynamic  Equilibrium),  rather  than  continue  to 
apply  LOWTRAN  /  FASCODE.  One  model  which  is  currently  becoming  operational  is  the 
AFGL  SHARC  (=  Strategic  High  Altitude  Radiation  Code),  which  calculates  atmospheric 
radiation  and  transmittance  for  paths  from  60  to  300  km  altitude  in  the  2-40  |J.m  spectral 
region.  The  initial  version  includes  the  five  strongest  IR  radiators,  NO,  CO,  H2O,  O3,  and 
CO2  -see  Sharma  et  al.,  1989. 

7.3  ATMOSPHERIC  VISIBILITY  AND  AEROSOLS 
7.3.1  Introduction:  Visibility 

In  the  visible  spectral  range,  how  far  one  can  see  is  generally  a  function  of  the 
number  density  of  small  particles  ("atmospheric  aerosols")  rather  than  of  the  gaseous  air 
species,  most  of  which  are  transparent  in  the  visible.^^  If  there  are  n  particles  per  unit 
volume,  each  of  mean  extinction  (=  absorption  +  scattering)  cross  section  o,  one  defines 
the  extinction  coefficient  k  (length'^) 

k  =  n  a  .  (7.9) 

An  effective  visual  range  Vc  is  given  by  the  "Koschmieder  relation" 

Vc  =  ln(100/C)/na  ,  (7.10) 


large  number  of  lines  are  needed  in  every  particular  specual  region  with  HITEMP,  thus  before 
auempting  to  use  this  code  you  should  check  to  be  sure  the  Atlas  has  been  completed  for  this  region. 

See  Pan  I,  Section  3.4:  Doppler  broadened  lines  have  a  gaussian  profile  which  falls  off  rather  rapidly 
away  from  the  line  center,  by  contrast  to  the  slow  fall-off  of  Lorentz  broadened  lines,  whose  line 
wings  can  affect  the  absorption  very  far  from  the  line  centers. 

Note  however  that  03  and  NO2,  both  of  which  are  often  present  in  small  concentration  in  polluted  air, 
do  absorb  in  the  visible. 
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where  C  is  the  (percent)  contrast  ratio 

contrast  ratio=  n(T+B)  -  I(B)]/I100  x  1(B)]  ,  (7.11) 

at  which  the  human  observer  can  detect  a  target.  Here  1(B)  is  the  intensity  of  the 
background  and  I(T+B)  is  the  intensity  of  a  target  seen  against  the  background.  C  =  2  to  5 
(percent)  are  normal  values  for  a  human  to  perceive  the  target;  traditionally,  meteorological 
visibility  is  defined  as  V2  for  a  contrast  ratio  of  2  percent  so  that  ln(50)  =  3.91,  but  in  fact 
an  "average"  observer  will  detect  targets  only  at  a  contrast  ratio  of  5  percent  or  greater,  so 
that  V5  =  ln(20)/na  =  3.00/na,  is  a  more  representative  visual  range.  Typical  values  of 
"good"  visibility  are  in  the  range  10-100  km.  For  a  discussion,  see,  e.g..  List,  1951, 
p.  452f. 

In  Section  7.3.2  we  discuss  very  briefly  the  different  types  of  aerosols  that  occur  in 
the  atmosphere,  while  Section  7.3.3  discusses  "Mie  Scattering,"  i.e.,  expressions  for  the 
absorption  and  scattering  cross  section  as  a  function  of  particle  size,  in  particular  of  the  Mie 
parameter 

q  =  2icaA  ,  (7.12) 

where  a  =  particle  radius  and  X  =  wavelength  of  light,  and  of  the  complex  refractive  index 
of  the  medium,  m  =  Ni  -  iNi.  Section  7.3.4  discusses  the  particular  problem  of  Rayleigh 
scattering  by  air  molecules  which  is  important  in  the  near  UV  because  if  q  «  1,  the 
scattering  cross  section  is  proportional  to  q^. 

For  the  problem  of  water  droplets  or  ice  crystals,  there  is  additional  physics  due  to 
the  condensation  and  evaporation  of  the  droplets.  In  this  context  it  is  customary  to  talk 
about  clouds  rather  than  aerosols;  elements  of  cloud  physics  are  discussed  in  Section  7.4. 

7.3.2  Atmospheric  Aerosols 

Atmospheric  aerosols  are  typically  of  0.03  to  0.5  pm  mean  radius.  They  can  be  of 
quite  different  types,  namely  soot  and  other  pollutants,  mineral  particles,  water-soluble 
material  such  as  sea  salt,  and  biotic  particles  such  as  bacteria  and  fine  plant  debris,  as  well 
as  water  and  ice  particles.  Table  7.3  gives  values  for  the  complex  refractive  index 

m  =  Ni-iN2  (7.13) 

of  some  of  these  materials  at  different  wavelengths;  Ni  is  the  real  pan  of  the  refractive 
index,  while  N2  is  the  absorption  coefficient.  Note  that  water  (and  ice)  absorbs  very  little 
in  the  visible  while  soot  is  a  highly  absorbing  material  at  all  wavelengths.  Dust  is 
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intennediate,  and  in  general  N2  or  the  absorption  coefficient  increases  with  increasing 
wavelength  in  the  IR.  Figure  7.3  gives  some  typical  values  for  the  atmospheric  extinction 
coefficient  in  the  visible  (X,  =  0.55  pm)  as  function  of  altitude.  Note  that  there  is  a  great 
deal  of  extinction  in  the  "Planetary  Boundary  Layer,"  roughly  the  lowest  1  km  of  the 
atmosphere,  and  a  secondary  maximum,  the  Junge  Layer,  in  the  15  to  20  km  altitude  range. 
The  Junge  layer  is  largely  of  volcanic  origin,  is  highly  variable,  and  is  made  up  mainly  of 
H2SO4/H2O  droplets  of  ~  1  pm  mean  size.  For  more  information,  see  USAF,  1985, 
pp.  18-9ff. 


Table  7.3.  Complex  Refractive  Index  m  =  Ni  - 1  N2  for  Different  Materials. 
(Source:  USAF,  1985,  pp.  18*16ff) 


Wavelencjth 

Water 

Dust-Like 

Soot 

(urn) 

Ni 

N2 

N, 

N2 

Ni 

N2 

0.25 

3.35E-8 

d 

3.00E-2 

d 

0.450 

0.55 

mm 

1 .96E-9 

8.00E-3 

119 

0.440 

1.06 

4.18E-6 

mm 

8.00E-3 

1.750 

0.440 

2.70 

1.188 

1 .90E-2 

RH 

1 .30E-2 

1.830 

0.520 

4.50 

1.332 

1 .34E-2 

mi 

1.40E-2 

1.940 

0.590 

10.0 

1.218 

5.08E-2 

0.162 

0.720 

20.0 

1.480 

0.393 

■ni 

0.220 

0.850 

7.3.3  Mie  Scattering  by  Small  (Spherical)  Dielectric  Particles*® 

We  examine  the  interaction  of  radiation  of  wavelength  X  with  spherical  particles  of 
radius  a,  and  introduce  the  Mie  parameter  q  =  27ta/X  of  Eq.  (7. 12).  The  particles  are  made 
up  of  a  medium  of  complex  refractive  index  m  =  Ni  -  i  N2.  The  scattering  cross  section 
o(sca)  and  the  absorption  cross  section  o(abs)  are  expressible  as  functions  of  q  and  m;  it  is 
possible  to  get  useful  results  in  a  variety  of  limits  that  are  presented  here  first.  We  also 
present  the  differential  scattering  cross  section  or  angular  distribution  P(0).  We  discuss 
three  limits,  first  the  "Rayleigh  Limit"  of  small  particles,  q  =  27ta/X  «1,  then  the 
"Geometrical  Optics"  limit  of  large  particles,  q  »  1,  and  then  the  limit  of  particles  of 
various  values  of  the  size  parameter  q  with  refractive  index  close  to  one. 


*  ®  The  present  results  come  mostly  from  Van  de  Hulst,  1957;  they  arc  presented  here  so  that  the  reader  can 
get  a  simple  semi-quantitative  model  for  scattering  from  small,  spherical  particles.  The  discussion 
emphasizes  dielecuic,  i.e.,  scattering  particles  rather  than  particles  (such  as  soot  or  metals)  that  are  • 

primarily  absorptive. 
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VERTICAL  AEROSOL  DISTRIBUTIONS 


Figure  7.3.  Representative  Extinction  Coefficients  for  Atmospheric  Aerosols 
as  Function  of  Altitude  (USAF,  1985,  p.  18*14) 

7.3.3. 1  Very  Small  Particles:  q«l  (Rayleigh  Limit) 

Here  the  scattering  cross  section  is 

o(sca)/;ta2  =  (8/3)q4  ((m2  -  l)/(m2  +  2)|2  ,  (7.14) 

and  if  m  close  to  1  and  also  N2  small,  so  that  there  is  very  little  absorption,  Eq.  (7.14) 
becomes 

o(sca)Ma2  =  (32/27)  q^  (Ni-1)2  .  (7. 15) 

In  this  limit,  provided  m  is  close  to  1,  and  N2/(Ni-l)  «1-  we  have 

o(abs)/7ia2  =  (8/3)q  N2  .  (7.16) 

Eq.  (7.16)  gives  o(abs)  ~  a^,  so  that  in  this  limit  the  absorption  cross  section  is 
proportional  to  the  volume  or  mass  of  material.  Thus,  every  individual  atom  or  molecule 
contributes  to  absorption/emission  if  the  particle  is  sufficiently  small  so  that  the  outer 
molecules  do  not  shield  the  inner  ones.  For  large  particles  where  this  shielding  is 
imponant,  a(abs)  ~  a2.  Note  that  in  the  Rayleigh  limit,  o(sca)/7ra2  ~  q^,  so  that  there  is 
very  little  scattering  if  q  «1. 
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7. 3. 3. 2  Very  Large  Particles:  q  »  1  (Geometrical  Optics  Limit) 


Here  ^ 

<y(sca)/ra2  =>  2  as  q  becomes  very  large  ,  (7.17) 

and 

a(abs)Ma7  =>  1  as  N2  q  becomes  very  large  .  (7.18)  ^ 

The  physical  meaning  of  these  two  independent  conditions  is  the  following.  The 
scattering  cross  section  for  large  particles  is  twice  the  geometrical  particle  size  because  of 
diffraction  or  shadow  scattering,  while  the  absorption  cross  section  can  be  as  large  as  the 
geometrical  particle  size  for  the  case  in  which  all  radiation  falling  on  the  particle  is  ^ 

absorbed.*' 


For  the  case  of  m  close  to  1 ,  and  N2  small,  it  is  convenient  to  define 

u  =  2qlm-ll~2q(Ni-l)  ,  (7.19) 

and  now 

a(sca)/;ca2  =  2  -  (4/u)  sin  u  +  (4/u2)  (1  -  cos  u )  (7.20) 

a(abs)/jca2  =  (8/3)qN2  .  (7.21) 


7. 3. 3. 3  Angular  Distribution,  P(0)*2 

If  q«l ,  so  that  the  particles  are  small  compared  to  the  wavelength  of  light,  there  is 
very  little  scattering  because  from  Eq.  (7.13)  a(sca)  ~  q^,  and  what  scattering  there  is  has 
the  Rayleigh  angular  distribution 

P(0)  =  (1  +  cos2  0)  .  (7.20) 

Here  P(0  deg)  =  2  =  P(180  deg)  and  P(90  deg)  =  1,  so  that  there  is  relatively  little  angular 
variation. 


By  contrast,  if  q»l  (geometrical  optics  limit)  and  m  is  close  to  1,  there  is  little 
refraction  and  absorption,  and  most  of  the  scattering  is  in  a  forward  direction.*^  If  q  >  1, 


*'  For  an  example  in  which  both  limits  (7.15)  and  (7.16)  apply  simultaneously,  see,  e.g..  Van  de  Hulst, 

1957,  p.  181,  Fig.  33. 

*  2  We  define  6  =  0  as  forward  scattering  (i.e.,  direct  transmission),  and  0  =  180  deg  as  back  scattering. 

As  an  example,  note  that  a  large  glass  bead  (or  cloud  water  droplet)  does  not  scatter  much  incident 

radiation  out  of  the  direct  beam,  so  that  most  of  the  radiation  incident  on  the  bead  is  transmitted,  with  • 

just  slight  refraction. 
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different  light  rays  going  through  the  particle  can  have  significant  differences  in  optical  path 
within  the  particle.  These  lead  to  interference  between  the  different  beams  and  thus  the 
angular  distribution  has  a  series  of  maxima  and  minima— see,  e.g.,  Fig.  7.4,  which  shows 
a  polar  plot  of  the  scattering  of  linearly  polarized  light  by  a  dielectric  sphere  of  refractive 
index  Ni  =  1.25.  The  quantities  ii  and  12  correspond  to  orthogonal  polarizations  of 
scattered  radiation.  However,  if  instead  of  having  a  medium  consisting  of  particles  all  of 
the  same  size  as  in  Fig.  7.4,  there  is  a  smooth  particle  size  distribution,  then  the  violent 
fluctuations  in  angular  distribution  are  smoothed  out  due  to  the  distribution  in  particle  sizes, 
and  one  gets  the  much  smoother  distribution  shown  in  Fig.  7.5,  which  has  actually  been 
calculated  for  water  droplets  (cloud  particles)  at  different  wavelengths  between  visible  and 
IR.  Note  that  the  shorter  the  wavelength,  i.e.,  the  larger  the  value  of  qa  =  27caa/7.,  the  more 
forward  scattering,  and  the  larger  is  the  scattering  cross  section.  (Here  aa  and  thus  qa  refer 
to  the  average  particle  radius  and  thus  average  value  of  q).  Thus  for  X  =  20  pm,  we  have 
qa  =  1.26,  and  the  cross  section  approaches  the  Rayleigh  form,  although  the  ratio  of 
forward  to  backward  scattering  cross  section  is  still  ~  30  (versus  3000  for  X  =  2.7  pm).  At 
the  shortest  wavelengths  there  is  some  back  (near-0  =180  deg)  scattering. 

Normally  one  cannot  find  a  simple  analytic  fit  to  the  various  curves  of  Fig.  7.5,  but 
the  Henyey-Greenstein  function 

Phg(0)  =  a  (1  -  g2)/(l  +  g2-  2g  cos0)3/2  (7.23) 

is  frequently  used  for  applications  such  as  Monte  Carlo  calculations  where  a  simple  analytic 
form  is  needed.  Here  g  =  (cos0)av.  and  it  and  the  scale  factor  A  are  chosen  to  fit  the  data. 
The  Henyey-Greenstein  function  does  not  show  the  (small)  enhanced  back  scattering 
shown  in  the  visible  in  Fig.  7.5— that  would  require  at  least  one  more  parameter. 

7.3.4  Rayleigh  Scattering  by  Air  Molecules  in  the  Visible  and  Near-UV 

At  wavelengths  below  about  0.4  pm,  atmospheric  transmission  is  reduced  due  to 
scattering  by  air  molecules  which  here  behave  like  very  small  particles  [with  q  «  1,  so  that 
a(sca)  ~  l/X^  and  is  very  small,  from  Eq.  (7.15)].  This  effect  is  responsible  for  the  blue  of 
the  sky.*^  The  effective  extinction  coefficient  is 

ks(^)  =  (29  jc2/33)  (Ni  -  1)2  n  (7ta2)3A4  ,  (7.24) 


When  we  look  away  from  the  sun  at  the  sky,  we  see  scattered  solar  radiation.  Because  o(sca)  -  lA.'*, 
there  is  much  more  short  wavelength  (blue)  sunlight  scattered  than  long  wavelength  (red)  sunlight  (by  a 
factor  ~  2^  =  16)  and  thus  the  sky  looks  blue. 
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where  n  =  number  of  air  molecules  per  unit  volume,  na^  =  mean  cross-sectional  area  of  a 
particle,  and  Ni  is  the  real  part  of  the  mean  refractive  index  of  air,  some  typical  values  are 

Ni  -  1  =  3.0  X  1(H  (0.25  urn) 

2.8  X  1(H  (0.5  pm) 

2.7  X  10-4  (2  pm)  . 


See  USAF,  1985,  p.  18-7  for  more  detail. 


7.4  PHYSICS  OF  CLOUDS 


The  saturated  vapor  pressure  (SVP)  of  water  increases  very  strongly  with 
increasing  temperature,  and  at  temperatures  below  freezing  the  SVP  over  ice  is  smaller  than 
that  over  water—see  Fig.  7.6.  Thus,  if  we  take  an  "air  parcel"  and  cool  it  adiabatically  (and 
without  changing  the  water  vapor  mixing  ratio),  the  relative  humidity  increases.  If  we  cool 
the  parcel  further,  eventually  it  saturates,  and  once  the  moisture  content  exceeds  saturation 
at  the  parcel  temperature,  the  excess  of  water  condenses.  Reference  to  Fig.  7.6  shows  that 
at  temperatures  below  0°C  the  SVP  over  water  is  greater  than  that  over  ice,  so  that 
condensation  over  ice  occurs  at  lower  water  vapor  mixing  ratios  than  over  water.  At  air 
temperatures  of  0°C  to  -20°C,  condensation  occurs  either  as  water  dioplets  or  as  ice 
crystals,  but  below  -20°C  ice  generally  predominates. 


Figure  7.6.  Saturated  Vapor  Pressure  of  Water  Vapor  Over  Water  and 
Over  Ice  as  Function  of  Temperature  (Source:  List,  1951,  p.  351  f) 


The  latent  heat  of  evaporation/condensation  of  water  is  of  the  order 
2.8  X  lO^  Joule/kg-see  Fig.  7.7--which  is  sufficiently  large  for  the  enthalpy  of  moist  air 
(with  a  water  vapor  mixing  ratio  of  order  several  time  \0-^)  to  be  significantly  larger  than 
that  of  dry  air.85  Furthermore,  the  settling  speed  of  small  water  droplets  (~  10  pm  or  less) 
shown  in  Fig.  7.8  is  sufficiently  small  compared  with  vertical  atmospheric  convection 


This  has  considerable  implications  for  the  climatic  difference  between  dry  and  moist  regions,  including 
such  things  as  the  day-night  temperature  difference  which  is  much  smaller  over  water  than  in  a  desert 
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speeds,  which  are  usually  in  the  range  1-10  cm/sec,  so  that  assemblies  of  air  containing 
water  droplets  or  ice  crystals-i.e.,  clouds— are  effectively  (if  not  thermodynamically) 
stable. 


I 


I 


■ 


ft 


ft 


ft 


Because  Cp(ice)  =  0.5  caVgm  K 

Cp(water)  =  l.Ocal/gmK 

Cp(watcr%-apor)  =  0.45caVgmK 

(cf:  AIP  Handbook Jnd  Ed.,p.4-189;  Hdbook  of  Chem.  &  Phys.,44th  Ed„  p.2436. 
Slater,  Intro,  to  Chem.  Physics,p.l47)  we  have  the  following  cycle: 


H2O  (liquid,  273  K)  H2O  (gas, 250  K) 


Note  that,  because  Cp(dry  air)  =  0.24  cal/gm  K,  the  enthalpy  change  corresponding  to  the 
evaporation  of  1  gm  ice  is  of  the  same  order  as  the  enthalpy  required  :c  raise  the 
temperature  of  1  kg  dry  air  by  3  K. 

Figure  7.7.  Latent  Heat  of  Evaporation/Condensation  of  Water 
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Clouds  are  always  produced  by  the  adiabatic  cooling  of  a  (moist,  warm)  air  mass. 
This  process  can  take  place  under  a  variety  of  different  conditions,  and  so  it  is  appropriate 
to  begin  by  identifying  the  concept  of  atmospheric  stability.  Figure  7.9  shows  a  generic 
atmospheric  temperature  profile  which  indicates  three  distinct  altitude  regions. 

Below  approximately  1  km  there  is  the  Planetary  Boundary  Layer  (PBL),*^  a 
region  in  which  the  temperature-altitude  profile  changes  dramatically  between  normal 
daytime  and  night  time  conditions.  Above  the  PBL,  and  up  to  a  temperature  minimum  near 
1 1  km,  is  the  Free  Troposphere,  a  region  in  which  the  temp)erature  normally  decreases  with 
height  (although  there  may  be  additional  temperature  layers  within  this  region).  Above  the 
temperature  minimum  near  1 1  km,  the  Tropopause,  comes  the  Stratosphere,  a  region  in 
which  the  temperature  increases  with  height. 

A  region  in  which  the  temperature  decreases  strongly  with  increasing  height  is 
dynamically  unstable  in  that  if  a  parcel  of  air  is  displaced  slightly  from  its  initial  position  it 
moves  further  in  the  direction  of  its  initial  displacement,  either  up  or  down.  The  free 
troposphere  is  such  an  unstable  region,  and  this  is  why  convective  or  cumulus  clouds  are 
often  found  in  this  region.  The  night  time  PBL  and  the  stratosphere  are  dynamically  stable; 
thus  generally  there  are  no  clouds  in  the  stratosphere,  and  typically  stable,  layered  clouds 
("stratus")  just  at  the  inversion  topping  the  PBL. 

It  can  easily  be  shown  that  a  parcel  of  air  is  stable  if  the  vertical  temperature  gradient 
for  an  adiabatic  displacement,  9T/dzladiab>  is  sufficiently  small: 

dT/dzIartiah  <  ~  E/Cp  —  Yad  »  .25) 

where  g  =  acceleration  due  to  gravity  and  Cp  =  specific  heat  of  air  at  constant  pressure. 
Conversely,  an  air  mass  is  unstable  if 

3T/9zladiab  > Yad  •  (7.26) 


A  good  definition  of  the  PBL  (provided  by  R.  PenndorO  is  "a  layer  in  which  a  noticeable  influence  of 
the  earth's  surface  on  wind  occurs  by  friction.  The  upper  boundary  is  determined  by  the  altitude  at 
which  the  wind  vector  is  approximately  equal  to  the  geostrophic  wind  (i.e.,  the  wind  determined  by  the 
overall  motion  of  the  atmosphere).  Within  the  PBL,  the  wind  direction  and  speed  are  affected  by 
friction  near  the  ground  in  establishing  a  universal  function,  the  "Ekman  spiral.”  The  upper  boundary 
of  the  PBL  is  frequently,  but  not  always,  characterized  by  a  temperature  inversion.  Clouds  are 
generally  missing  in  the  PBL,  but  haze  or  mist  may  form.  Only  just  below  the  upper  boundary 
turbulent  clouds  may  form.  There  is  a  lot  of  turbulence  in  the  PBL  due  to  heating  of  the  ground,  and 
adiabatic  or  even  stronger  temperature  gradients  may  develop  as  a  result  of  radiative  processes  involving 
the  ground.  If  the  temperature  inversion  at  the  top  of  the  PBL  is  strong,  it  may  block  mixing  of  the 
air  from  above.” 
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Temperature  Ocean 

Figure  7.9.  Generic  Temperature  Profile  In  the  Atmosphere 


Representative  values  for  the  Adiabatic  Lapse  Rate  Yad  =  “  g/Cp  arc 

Yad  =  “  9.8  °C/kni  for  dry  air  (dry  adiabatic  lapse  rate)  (7.27 a) 

-  6.5  °C/km  for  moist  air  (moist  adiabatic  lapse  rate)  .  (7.27b) 

In  general  there  are  three  major  types  of  clouds: 

•  Cumulus  (Cu):  Convective  clouds  of  extensive  vertical  development  under 
unstable  conditions.  They  occur  generally  in  the  free  troposphere,  although 
very  strong  thunderheads  or  cumulonimbus  can  rise  to  the  tropopause  or  even 
into  the  lowest  stratosphere. 

•  Stratus  (St):  Layered  clouds  of  uniform  horizontal  development,  which  are 
frequently  found  near  the  top  of  the  PBL. 

•  Cirrus  (Ci):  Thin,  high-altitude  clouds  made  up  of  ice  crystals,  which  typically 
occur  in  the  upper  troposphere  and  may  move  into  the  low  stratosphere. 
Cirrus  form  in  a  variety  of  ways,  such  as; 

(a)  by  blowoff  from  cumulonimbus,  or 

(b)  when  a  warm  and  therefore  moist  air  mass  approaches  a  cold  air  mass 
near  a  frontal  zone. 
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There  is  evidendy  no  unique  characterization  of  these  different  morphological  types 
of  clouds.  Figure  7.10  sketches  different  types  of  clouds,  as  viewed  from  above. 
Figures  7.10a,  b,  c  show  "representative"  cumulus,  stratus,  and  cirrus  clouds,  while 
Fig.  7.10d  shows  an  actual  cloud  configuration  observed  one  day  over  the  U.S.  Midwest 
from  a  NASA  high-altitude  research  aircraft;  notice  that  actual  clouds  tend  to  be  a 
combination  of  the  various  idealized  types.  Lehr  et  al.,  1965,  present  an  excellent 
elementary  discussion  of  meteorology,  including  good  sketches  of  clouds  and  a  discussion 
of  different  cloud  types.*^ 

Table  7.4  gives  some  representative  quantitative  characteristics  of  different  types  of 
clouds.  There  are,  of  course,  many  other  types  of  clouds,  and  generally  low-altitude 
clouds  occurring  in  the  relatively  warm  and  moist  atmosphere  have  a  higher  condensed 
water  content  (>  0.3  give?  of  condensed  water)  than  do  higher-altitude  clouds  that  form  in  a 
colder  and  drier  ambient  atmosphere.  Note  that  if  a  cloud  is  to  form,  the  ambient  air  must 
be  saturated  with  water  vapor  in  addition  to  the  condensed  water  in  the  droplets,  and 
typically  the  water  vapor  accounts  for  most  of  the  atmospheric  moisture  in  the  cloud,  while 
the  particles  may  make  up  5  to  20  percent  of  total  H2O  content. 

The  cloud  particle  size  distribution  is  quite  variable  but  is  generally  continuous,  so 
that  the  mean  particle  scattering  cross  section  has  the  character  of  Fig.  7.5,  i.e.,  the  violent 
oscillations  as  a  function  of  q  =  2jra/X  shown  in  Fig.  7.4  are  damped  out  because  of  the 
variation  in  particle  size.  One  representative  particle  size  distribution  is  that  introduced  by 
Deirmendjian,  1964: 

dN(a)/da  =  Ndabe-«a  ,  (7.28a) 

where  a  =  particle  radius,  Nd  is  proportional  to  the  number  of  particles  per  unit  volume,  and 
b  and  c  are  parameters.  Thus,  the  mean  radius  aav  is 

aav  =  b/c  ,  (7.28b) 

and  the  mean  number  density  is 

nav  =  Nd  b! /c»>+l  .  (7.28c) 


An  excellent  description  of  clouds  is  given  by  Scorer,  1972.  Note  that  there  is  also  "subvisual  cirrus" 
which  occurs  not  infrequently  at  high  altitudes.  This  is  10-100  times  less  dense  than  "regular"  cirrus 
and  cannot  be  seen  with  the  unaided  eye;  it  shows  up  on  some  lidars  and  can  interfere  with  long-range 
near-horizontal  viewing. 
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(a)  Typical  cumulus  clouds 


(b)  Typical  stratus  clouds.  Relatively  smooth  surface 
producing  benign  background 

Figure  7.10.  Some  "Representative"  Cloud  Types 


(c)  Typical  cirrus  clouds.  "Fail  streaks"  of  the  relatively  large 
Ice  crystals  making  up  these  clouds  are  typical. 


(d)  Sketch  of  an  actual  cloud  configuration  photographed  by  a  WB-57  aircraft 
flying  at  60,000  ft  over  the  Midwest.  Note  cumulus,  stratus,  and 
cumuilform  edges,  haze  and  cirriform  clouds. 


Figure  7.10.  (continued) 
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For  water  clouds,  aav  =  4  ,  b  =  6,  c  =  1.5  are  plausible  values,  while  for  ice 

clouds,  aav  =  32  pm  ,  b  =  6,  c  =  0.1875  |im~^  may  be  used. 

Note  that  clouds  are  highly  variable  but  not  random.  Cumulus  clouds  will  change 
as  observed  over  a  10-  to  30-min  period,  which  is  the  typical  lifetime  of  an  individual 
cumulus  convection  cell.  Cirrus  change  noticeably  over  times  of  1  to  3  hours,  while  stratus 
clouds  are  relatively  long-lived,  with  lifetimes  are  in  the  3-10  hr  range.  The  bottom  surface 
of  a  cloud  is  normally  fairly  sharp  as  compared  with  the  top,  because  its  position  is 
determined  by  the  local  condensation  temperature  in  the  atmosphere  of  the  air  parcel  which 
gives  rise  to  the  particular  cloud  system;  by  contrast,  the  top  of  a  convective  cloud  system 
shows  the  characteristic  shape  due  to  turbulent  convection. 

Table  7.5  gives  the  time-space  scales  of  different  types  of  meteorological 
phenomena  which  reflect  on  the  underlying  atmospheric  motions.  Atmospheric  motions 
are  often  described  as  a  combination  of  advection  (=  convection)  and  turbulent  diffusion, 
but  there  is  no  very  clear  distinction  between  them,  except  perhaps  time-space  scales. 
Representative  wind  speeds  are  10-100  m/sec  horizontally,  and  0.01-0.1  m/sec  vertically; 
Fig.  7.11  shows  "representative"  horizontal  winds,  which  vary  strongly  with  altitude. 
Atmospheric  diffusion  is  highly  scale-dependent;  Fig.  7.12  shows  horizontal  mean 
diffusive  (i.e.,  non-convective)  motions.  Note: 

(a)  The  wide  range  of  variability  of  data  from  the  mean  curve  HI  to  the  lower  and 
upper  bounds  I  and  IV  respectively. 

(b)  The  scales  of  atmospheric  motions  from  Table  7.4  correspond  quite  well  with 
the  bounds  I  and  FV  of  the  figure. 

(c)  These  data  are  mainly  tropospheric  (1-10  km  altitude)  but  apply  up  to 
20-25  km. 

(d)  An  older  range  of  variability  (Bauer,  1974)  with  mean  curve  11  and  bounds  I 
and  III  should  not  be  used,  since  it  is  based  on  older  observations  which  in 
general  are  biased  towards  quiet  atmospheric  conditions.  (The  older 
observational  techniques  had  limited  sensitivity,  so  that  under  non-quiescent 
conditions  the  tracer  cloud  would  be  dispersed  so  rapidly  that  no  measurements 
could  be  made.) 
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Table  7.5.  Representative  Scales  of  Atmospheric  Motions 
(Source:  Bauer.  1983) 
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Figure  7.11.  Steady-State  Horizontal  Cumulative  Wind  Speed  Profiles  (m/sec) 
for  Cape  Canaveral  Space  Center  (Source:  Turner  and  Hill,  1982)  (At  i  km 
altitude,  99%  of  the  time  the  wind  speed  Is  <  24  m/sec). 
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Figure  7.12.  Horizontal  Dispersion  as  Function  of  Travel  Time.  The  Blocks  Hi, 
H2,...,R1,  R2,...  come  from  Table  7.5.  (Source:  Bauer,  1983) 


7.5  CLOUD-FREE  LINE-OF-SIGHT  (CFLOS)  ANALYSES 


7.5.1  Introduction 

We  recognize  that  clouds  occur  very  frequently— see,  e.g..  Table  7.1— and  are 
highly  variable-see,  e.g..  Fig.  7.10.  Visible  and  infrared  radiation  cannot  penetrate 
optically  thick  clouds  directly,  although  the  brightness  of  visual  illumination  under 
complete  cloud  cover  is  as  much  as  1/8  to  1/4  of  the  full  solar  illumination.  For  instance,  to 
take  a  photograph  under  cloudy  skies  one  increases  the  aperture  by  just  2-3  stops,  but 
reference  to  Table  7.3  shows  that  water  clouds  are  very  much  more  absorptive  in  the 
infrared  than  in  the  visible. 

A  number  of  specific  questions  arise  on  how  often  there  are  clouds  in  the  field  of 
view  of  the  human  eye,  or  of  a  visible  or  infrared  sensor,  which  interfere  with  its  seeing  a 
given  object  ("target").  This  is  normally  treated  as  a  probabilistic  issue,  but  we  shall  begin 
by  recalling  what  clouds  look  like  in  a  litde  more  detail  than  the  generic  model  presented  in 
the  previous  section,  and  how  the  clouds  interfere  with  optical  sensing  of  a  target. 

Table  7.6  is  a  "do-it-yourself  example  of  what  clouds  "actually"  look  like.  I 
decided  to  look  out  of  my  office  window  or  equivalent  for  19  successive  days  (actually 
between  13  and  31  May  1988)  at  10  a.m.  and  3  p.m.  and  ask: 

(a)  What  fraction  of  the  sky  is  covered? 

(b)  What  type  of  cloudiness  is  observed? 

While  the  results  are  not  quantitatively  significant,  yet  they  give  the  interested  reader  a  good 
feeling  for  how  much  cloud  cover  there  actually  is  and  for  the  complexity  and  variability  of 
real  clouds,  before  he  plunges  into  a  probabilistic  analysis  of  Cloud-Free  Line-of-Sight 
(CFLOS). 

One  further  remark  should  be  made  here.  For  surveillance  from  high  altitude  one  is 
interested  in  local  holes  in  cloud  cover  which  enable  one  to  observe  a  designated  target 
from  above.  A  "hole"  in  cloud  cover  is  an  inverse  of  the  presence  of  a  cloud  only  in  a 
mathematical  rather  than  a  physical/meteorological  sense.  Cumulus  clouds  are  ~  1-10  km 
in  linear  dimension,  spaced  a  comparable  distance  apart,  but  stratus  cloud  decks  rarely  have 
cloud-free  holes.  Rather,  a  hole  in  cloud  cover  is  often  to  be  interpreted  in  terms  of 
Fig.  7.10,  as  a  gap  between  overlapping  layers  of  low,  middle,  and  high  clouds  of 
different  types.  Statistics  of  clear  versus  cloudy  intervals  under  different  levels  of  total 


cloud  cover  are  listed  in  Table  7.7.  A  somewhat  different  representation  of  Cloud-Free 
Intervals  (CFIs)  at  different  locations  is  given  by  Malick  and  Allen,  1978, 1979,  and  this  is 
discussed  in  Section  7.5.3  below. 

Having  given  this  phenomenological  introduction  to  what  clouds  actually  look  like, 
let  us  now  address  the  probabilistic  issue  which  has  been  addressed  by  a  number  of 
different  data  bases  and  methodologies.  Here  are  listed  some  specific  questions  of  interest, 
representative  data  bases  used  in  answering  them,  and  relevant  analyses.  The  examples 
listed  here  are  intended  to  be  representative  rather  than  complete. 

Specific  questions  that  are  often  addressed  by  the  DOD  community  deal  with  the 
probability  of  the  following: 

•  Clear  and  cloud-free  conditions,  including  effects  of  haze  or  thin  cirrus 

•  Cloud-Free  Field  of  View  (CFFOV) 

•  Cloud-Free  Interval  (CFI)  on  a  given  path 

•  Dynamic  Cloud-Free  Line-of-Sight  (DCFLOS)— which  considers  the  temporal 
motion  of  cloud,  target,  or  sensor.** 

For  current  work  in  this  field,  the  reader  is  referred  to  the  annual  reports  of  the 
CIDOS  (=  Cloud  Impacts  on  DOD  Operations  and  Systems)  Conference.*^ 

7.5.2  Data  Bases 

•  De  Bary  and  Moeller,  1963,  German  Weather  Service  aircraft  flights  over 
Germany  in  the  1930s. 

•  Lund  and  Shanklin,  1972,  AFGL,  all-sky  camera  pictures  at  Columbia, MO. 
(This  data  base  is  being  updated,  beginning  in  1988,  by  using  six  all-sky 
vidicon  cameras  at  various  locations  in  the  USA.) 

•  Bertoni,  1977,  AFGL,  extensive  pilot  reports  of  clear  and  cloud-free  lines  of 
sight  in  the  northern  hemisphere  (data  should  be  used  with  great  caution). 

•  Avaste  et  al.,  1974— see  also  Ellingson,  1982— Soviet/East  European  data  base. 

•  National  Weather  Service-a  wide  variety  of  NWS  and  USAF  data  are  available 
from  National  Climatic  Center,  Asheville,  NC. 


**  This  problem  is  so  complicated  that  analyses  such  as  that  of  Nelson  and  Wetherbe,  1976  are  not  used 
extensively. 

Published  by  AFGL;  contact  Dr.  Wm,  Snow,  617-377-5952. 
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Table  7.6.  A  Do-lt'Yourself  Actual  Example  of  Cloudiness: 
Cloudiness  During  13-31  May  1988  in  the  Washington,  DC,  Area 


Date  (May.  1988) 

10  a.m. 

3  p.m. 

■Hi 

Cover*  (tenths) 

mm 

Cover  (tenths) 

F 

mm 

Cs 

8 

Cs 

6 

S 

mm 

Sc 

8 

Sc 

8 

Su 

15 

Sc 

6 

Sc 

8 

M 

16 

St 

10 

As 

5 

Tu 

17 

Ac 

5 

Ac 

7 

W 

18 

Sc 

10 

Sc 

9 

Th 

19 

Cu 

5 

Cu 

5 

F 

20 

St 

10 

St 

10 

S 

?1 

Sc 

6 

Cu 

6 

Su 

22 

Cu 

2 

Cu 

3 

M 

23 

Ci 

1 

Cu 

3 

Tu 

24 

Ci 

1 

Cu 

6 

W 

25 

Ns  (rain) 

10 

St 

4 

Th 

26 

-- 

0 

- 

0 

F 

27 

-- 

0 

- 

0 

S 

28 

Cu 

2 

Cu 

6 

Su 

29 

-- 

0 

Cu 

2 

M 

30 

Ci 

2 

Cu 

3 

Tu 

31 

- 

0 

Cu 

2 

Average  Cloudiness;  10  a. m.:  0.45;  3  p.m.:  0.49 


Notation  used  is  standard:  Ci  (cirrus),  Cs  (cirrostratus),  Cc(  cirrocumulus)  are  thin  high  clouds;  As 
(attostratus),  Ac  (altocumulus)  are  thicker  "middle*  clouds  (1-5  km  altitude);  Ns  (nimbostratus)  are  rain 
clouds;  Cu  (cumulus)  are  convective  clouds  of  relatively  large  vertical  extent;  Cb  (cumulonimbus)  are 
thunderheads;  St  (stratus)  and  Sc  (stratocumulus)  are  low  clouds,  largely  of  horizontal  extent. 

In  general,  stratiform  clouds  (St,Sc)  are  evaporated  by  solar  radiation  as  the  day  passes,  i.e.,  there  is 
more  stratus  in  the  morning  than  in  the  afternoon.  However,  cumuliform  clouds  (Cu,  Cb)  are  produced 
by  sofar-induced  convection,  so  that  they  are  more  likely  to  occur  during  the  afternoon/evening  than  in 
the  morning. 


Table  7.7  Cloud  Cover:  Mean  Horizontal  Sizes  of  Clouds  and  Spaces  Between 
Them.  [Source:  Spiridonova  (l976)-see  Matveev,  1984,  p.  253  ff] 


•  Analysis  of  TV  cloud  pictures  from  Cosmos  and  Meteo.'  satellites  between  1967 
and  1971,  over  latitudes  from  the  equator  to  60  deg  N.and  longitudes  between 
160  deg  E  and  70  deg  W  (North  America  and  Adjacent  Pacific  Ocean, 
approximately  one-third  of  the  area  of  the  northern  hemisphere).  Daylight  only, 
minimum  resolution  4  km. 

•  Data  are  broken  down  by  season  (January  or  July),  by  land  or  sea,  and  by 
Scattered,  Broken,  Overcast,  and  Total  Cloud  cover,  which  correspond, 
respectively,  to  2/10-5/10,  5/10-8/10,  8/10-10V10.  and  10/10  cloud  cover, 
where  10V10  is  "total  cloud  cover  with  gaps."  The  figures  are,  respectively,  cloud 
size  and  cloud  spacing,  all  in  km.  (Thus,  for  scattered  cloud  cover  over  ocean  in 
January,  the  mean  cloud  size  is  13  km  with  10  km  standard  deviation,  and  a  mean 
cloud  spacing  of  40  km.  Evidently,  for  10/10  cloud  cover  over  the  ocean  in 
January,  the  mean  size  of  a  solid  cover  is  190  km.) 

Mean  Cloud  Size/Spacing 


Cloud  cover 


Scattered(2-5/10) 
Broken  (5-8/10) 
Overcast  (8-10*/10) 
Total(  10/10) 


Ocean 

msm 

HBfflHIil 

13/40 

8/78 

29/21 

17/60 

46/19 

49/13 

190 

76 

Janua 


16/47 


29/31 


140/60 


10/64 


21/21 


Standard  deviation  of  cloud  sizes: 


Scattered  (2-5/10) 

10 

9 

15 

9 

Broken  (5-8/10) 

27 

29 

26 

19 

•  USAF  3DNEPH  and  RTNEPH-see  Fye,  1978,  and  Kiess  and  Cox,  1988. 
These  use  mainly  satellite  imagery  from  DMSP  and  NOAA  polar-orbiting 
•  satellites,  supplemented  by  ground  and  pilot  reports  to  provide  global  coverage 

(from  1975  on)  on  a  spatial  resolution  of  26  nm  (48  km)  updated  every  six 
hours. 


7.5.3  Some  Representative  Analyses 

•  Lund,  1965,  1966,  AFGL,  clear  LOS  from  sunshine  and  cloud-cover 
observations. 

•  Lund  and  Shanklin,  1973,  AFGL,  use  Lund  and  Shanklin,  1972  data  base. 

•  Lund,  Grantham,  and  Elam,  1975-h,  AFGL/ETAC,  global  atlas  using  Lund- 
Shanklin  and  other  ground-based  data. 

•  Rapp,  Schutz,  and  RotLiguez,  1973  conducted  an  error  analysis  of  various 
"early"  CFLOS  analyses,  and  concluded  that  the  Lund-Shanklin  model  was  the 
best  available  one. 

•  Malick  and  Allen,  1978,  1979,  SRI/DARPA,  combine  3DNEPH  global  data 
and  the  Lund-Shanklin  analysis. 

•  De  Violini,  Shlanta,  and  Elam,  1980,  USN,  global  marine  atlas  using 
3DNEPH  data. 

•  Katz,  DeBold,  and  Perez-Esandi,  1980,  NSWC,  surface-to-air  data  using 
surface  marine  observations  from  NWS  and  Lund-Shanklin  analysis. 

•  USAF-AWS  provides  a  variety  of  outputs  from  RTNEPH. 

Figure  7.13  shows  the  probability  of  a  cloud-free  line  of  sight  (PCFLOS)  as  a 
function  of  elevation  angle  0  and  sky  cover  S.  This  compares  the  modification  by  Malick 
and  Allen,  1978,  1979,  of  the  Lund  and  Shanklin,  1972,  work  which  may  be  considered 
as  the  "U.S.  Standard"  model  with  the  Soviet/Eastem  European  model  of  Avaste  et  al., 
1974  (see  also  Ellingson,  1982).  Note  that  for  near- zenith  viewing  (9  large). 


PCFLOS  =  1-S  , 


(7.29) 


but  for  oblique  viewing  (0  <  30-40  deg)  the  PCFLOS  becomes  very  much  smaller  (all 
curves  converge  to  zero  at  0  =  0  ).  The  difference  between  the  U.S.  and  Soviet  models 
tends  to  reflect  on  the  fact  that  the  Soviet  Union  is  a  very  cloudy  and  high-latitude  country. 

The  analytic  fit  of  Malick  and  Allen,  1978,  postulates  an  average  cloud  thickness- 
to-width  ratio  b  which  depends  on  the  fractional  sky  cover  S  through  the  relation 


b  =  0.55  -  S/2 


(7.30) 
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The  fractional  sky  cover  S  applies  to  visual  observation  over  the  hemisphere  looking  up 
from  the  ground,  as  contrasted  to  the  fractional  cloud  cover  Q  as  seen  looking  down  from 
space.  According  to  Malick  and  Allen,  1978,  these  two  quantities  are  related  by  the 
empirical  relation 

Q  =  (S/4)  (1  +  3  S)  .  (7.31) 

Now  Malick  and  Allen  parameterize  the  Lund-Shanklin  results  as  a  function  of  elevation 
angle  0  as  follows: 

PCFLOS  =  (1  -  Q)  (1  +  b/tan0)  .  (7  32) 

The  Malick  and  Allen  model  is  attractive  because  it  provides  simple  and  generally  usable 
results,  but  it  clearly  oversimplifies  the  real  situation.  However,  where  this  kind  of  simple 
model  is  not  adequate,  it  may  be  appropriate  to  use  numerical  results  from  the  RTNEPH 
data  base. 

7.5.4  Horizontal  Viewing 

Most  CFLOS  analyses  involve  upward  looking  observers  or  sensors  who  measure 
sky  cover,  or  else  downward  (nadir)  viewing  satellites  that  measure  earth  cover.  This  is 
satisfactory  for  viewing  most  of  the  sky  dome,  but  it  fails  for  (near)-horizontal  viewing,^ 
such  as  for  an  IRST  (Infrared  Search  and  Track)  application.  Here  it  is  more  appropriate  to 
use  results  from  the  NASA  SAGE  (Stratospheric  Aerosol  and  Gas  Experiment)  and  SAM 
satellites  which  use  near-visible  sensors  viewing  the  sun  horizontally  through  the  "eanh 
limb",  and  thus  provide  direct  experimental  data  on  the  probability  of  a  cloud-free  line  of 
sight  along  a  horizontal  earth  limb  path  as  a  function  of  tangent  height’* 

The  SAGE  and  SAM  limb-viewing  solar  photometers  are  flown  on  separate 
satellites  to  provide  global  coverage  and  are  intended  primarily  to  study  stratospheric 
aerosols  and  gases  (see  McCormick  et  al.,  1979).  However,  observations  go  down  to 
tangent  heights  as  low  as  5  km  and  they  show  an  enhanced  extinction  at  the  lower  altitudes 
which  is  due  to  high-altitude  clouds-see  Figure  7.14.  Livingston  and  Malick,  1983,  have 
analyzed  the  initial  1.5  years  of  SAGE  data  to  determine  the  PCFLOS  as  a  function  of 
tangent  height.  Figure  7.15  shows  their  results  as  a  function  of  latitude  and  season. 


Thus  note  that  Fig.  7.13  shows  no  results  for  elevation  angles  6  <  10  deg  above  the  horizon. 
’  *  The  effective  path  length  in  the  aunosphere  is  of  the  order  of  200  km. 
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Figure  7.14.  Limb-viewing  observations  of  aerosoi-cloud  extinction  from  the 
NASA  SAM-II  sensor  at  X  =  1  jam  (SAGE  is  similar).  The  maximum  In  extinction 
In  (a)  at  7  km  is  presumably  due  to  high-altitude  clouds. 

(Source;  McCormick  et  ai.,  1982) 


TANGENT  ALriTUOE  —  km  ASL 


Figure  7.15.  Seasonally  Derived  1.0  um  Transmission  Curves  From  SAGE 
Plotted  as  a  Function  of  Tangent  Altitude  and  Probability  for 
Three  Northern  Hemisphere  Latitude  Bands. 

(Source:  Livingston  and  MalIcK,  1983) 
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It  is  possible  to  extract  information  on  the  frequency  of  occurrence  of  a  cloud-free 
line  of  sight  in  horizontal  viewing  as  a  function  of  tangent  height  on  the  assumption  that  the 
low -altitude  extinction  shown  in  Fig.  7.14  is  due  to  high-altitude  clouds  located  near  the 
tangent  region,  i.e.,  that  the  frequency  of  cloudiness  decreases  monotonically  with 
increasing  altitude.  Some  results  for  summer  and  winter  at  latitudes  of  30-50  deg  N  are 
given  in  Fig.  7.16.  As  an  example  of  how  to  use  these  data,  consider  winter  conditions, 
and  a  transmittance  of  0.5.  Then  the  minimum  tangent  height  ztan  which  can  be  achieved 
10  percent,  50  percent,  and  90  percent  of  the  time  is,  respectively,  4.4  km,  8.0  km,  and 
11.0  km,  and  for  a  sensor  at  height  Zq  viewing  a  target  at  the  same  altitude,  the  range 
achievable  s,  is  related  to  zq,  Ztan.  and  the  radius  of  the  earth,  R©  =  6400  km  by  the  simple 
geometrical  relation 

s  =  2y[2Ro(Zo-z^)]  .  (7.33) 

Thus  we  see  that  under  the  listed  conditions,  the  maximum  range  achievable  10  percent, 
50  percent,  and  90  percent  of  the  time  is,  respectively,  620  km,  450  km,  and  230  km  for 
both  target  and  sensor  at  height  z©  =  12  km  (39,000  kft.). 


7.6  ATMOSPHERIC  TURBULENCE 

The  atmosphere  is  not  a  homogeneous  optical  medium.  The  temperature,  pressure, 
and  humidity  are  not  spatially  or  temporally  uniform,  and  this  leads  to  random  variations  in 
the  index  of  refraction  n  (r,t)  seen  by  a  propagating  electromagnetic  wave.  In  the 
atmosphere,  the  refractive  index  n  differs  only  very  slightly  from  unity;  a  typical  value  at 
sea  level  is  1.0003,  but  fluctuations  on  the  order  of  one  percent  in  the  3  x  10~4  residual  are 
sufficient  to  produce  beam  wander,  beam  spread  and  intensity  fluctuations  in  a  laser  beam, 
and  also  limit  astronomical  and  other  "seeing"  through  the  atmosphere.  These  phenomena 
arise  because  of  fluctuations  in  the  residual,  (n-1),  due  to  atmospheric  turbulence,  which 
accumulate  along  an  optical  path  through  the  atmosphere. 

The  refractive  residual  (n-1)  is  proportional  to  the  density,  and  thus 

(n-l)  =  AP/T  ,  (7.34) 

where  A  ~  7.9  x  l(>-5  (K/millibar)  for  red  light  (cf.,  HELPH-I,  p.  5-3).  It  is  appropriate  to 
express  changes  8n  in  terms  of  changes  60  in  potential  temperature 


Tangent  Height  (km) 


0.0  0.2  0.4  0.6  0.8  1.0 

Transmittance  in  Horizontal  Viewing 

(a)  Winter 


Transmittance  in  Horizontal  Viewing 


■e -  90% 

-  75% 

■a —  50% 

-• -  25% 

■■ -  10% 


(b)  Summer 


Figure  7.16.  Probability  of  Attaining  a  Given  Tangent  Height  In  Horizontal 
Viewing  as  a  Function  of  Horizontal  Atmospheric  Transmittance. 
Results  for  Winter  and  Summer,  30  deg  to  50  deg  N  latitude 
Bands  Computed  from  the  Results  of  Fig.  7.15. 
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where  z  =  altitude  and  Yad  =  adiabatic  lapse  rate-  see  Eq.(7.27)-because  (unlike  the 
temperature  T),  the  potential  temperature  0  is  conserved  in  the  atmosphere.  Thus 

5n  =  -  A(P/T2)  80  ,  (7.36) 

while  at  constant  altitude  z: 

80  =  8T  .  (7.37) 

In  a  turbulent  atmosphere,  the  local  temperature  is  a  random  variable,  and  it  is 
customary  to  introduce  a  structure  function 

r>r(r2-ri)  =  <[T(r2)-T(ri)]2>  (7.38) 

where  the  brackets  <  >  imply  an  ensemble  average  in  which  all  possible  point  pairs  ei,  12 
are  averaged.  By  imposing  the  condition  of  isotropy  of  turbulence,  we  remove  dependence 
on  the  vector  coordinates  rj  and  n  and  only  need  to  consider  the  magnitude  r  =  1 12  -  Ii  1  of 
the  difference  between  the  two  coordinates.  Kolmogorov,  1972,  demonstrated  by  the  use 
of  dimensional  analysis  that  as  long  as  r  is  restricted  to  the  "inertial  subrange,"  Eq.  (7.38) 
has  a  r2/3  power  law  dependence 

DT(r)  =  Ci^r2/3  .  (7.39) 

Likewise,  from  Eqs.  (7.36)  and  (7.37)  the  index  of  refraction  structure  function 

Dn(r)  is 

Dn(r)  =  Cn2  =  Of- r2/3  (AP/T2)2  .  (7.40) 

The  parameter  Cn^  is  frequently  used  to  describe  atmospheric  turbulence. 

Figure  7.17  is  a  sample  profile  of  the  variation  of  Cn^  from  the  ground  to  30  km 
altitude.  It  also  shows  the  simultaneously  measured  wind  speed  profile.  This  profile  was 
taken  during  a  jet  stream  passage  and  is  thus  an  instructive  example  of  some  of  the 
characteristics  of  atmospheric  turbulence.  A  typical  feature  shown  is  the  occurrence  of  the 
largest  values  of  Cf  at  the  surface  with  a  rapid  falloff.  In  addition  to  the  jet  stream 
enhancements,  the  other  significant  characteristic  of  atmospheric  optical  turbulence 
demonstrated  here  is  the  stratification.  Turbulence  always  occurs  in  layers  in  which  Cf 
can  exceed  the  "background"  value  by  more  than  100.  In  some  applications  of  vertical 
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Figure  7.17.  Optical  Turbulence  In  the  Atmosphere:  a  Sample  of  the  Parameter 
Cn‘  as  Function  of  Altitude.  (Data  Provided  by  R.R.  Beland,  AFGL) 


propagation  this  layering  is  not  important,  but  in  others,  particularly  for  horizontal 
propagation,  the  layering  may  be  very  significant.^^ 

The  discrete  layers  of  atmospheric  turbulence  shown  in  Fig.  7.17  are  characterized 
by  the  Richardson's  number  Ri— see  Section  5.6.4— rather  than  by  the  Reynolds'  number 
Re  of  Section  5.6.3— which  characterizes  aerodynamic  rather  than  atmospheric  turbulence. 
For  a  simple  discussion  of  atmospheric  turbulence,  see  USAF,  1985,  p.  19-6ff. 

7.7  NON-LINEAR  EFFECTS  FOR  HIGH-INTENSITY  LIGHT  BEAMS 

7.7.1  Introduction 

When  a  high-intensity  beam  of  radiation  passes  through  a  medium,  even  a  small 
fractional  absorption  can  lead  to  significant  effects  on  propagation.  For  instance,  consider 
a  10-MW  laser  beam  propagated  from  a  3-m  diameter  mirror  in  the  sea  level  atmosphere. 
Even  if  we  assume  an  absorption  coefficient  as  low  as  0.01  m“^  which  corresponds  to  a 
mean  molecular  absorption  cross  section  of  4  x  10“24  cm^  or  10“^  times  gas  kinetic,  the 
heating  rate  is  170  K/sec  because  Cp  =  1.0  J/g  K.  It  is  thus  clear  that  the  propagation  of  a 
high-intensity  beam  of  radiation,  such  as  from  a  high-energy  laser,  deposits  a  large  amount 
of  energy  in  the  atmosphere  that  can  lead  to  significant  changes  in  the  ambient  medium 
which  in  turn  affects  beam  propagation. 

To  some  extent  it  is  possible  to  correct  adaptively  for  these  changes  in  the 
propagation  medium.  These  corrections  will  not  be  discussed  here  (see,  e.g.,  APS,  1987) 
but  I  introduce  the  relevant  atmospheric  concepts  for  adaptive  correction  for  turbulence  and 
other  disturbances  of  the  propagating  atmosphere.  These  are  the  (phase)  coherence 
distance  Xq  and  the  isoplanatic  angle  %.  In  the  absence  of  turbulence  or  other  disturbances 
of  the  medium,  the  phase  variation  4)(r)  across  the  apenure  of  a  beam  is  determined 
uniquely  by  its  initial  distribution  across  the  launching  aperture;  explicitly,  if  a  coherent 
laser  beam  of  wavelength  X  is  launched  from  an  aperture  of  diameter  D,  after  traveling  a 
path  length  L  its  area  at  the  target  will  be 

'^ideal  “  »  (7.41) 

but  in  a  turbulent  medium  of  coherence  distance  Tq,  the  effective  beam  area  will  be 

Note  that  frequently  used  models  of  Cn^  are  smoothed  fits  to  averaged  data  and  thus  do  not  exhibit  this 
layering. 


A~Aidea][  l+(D/ro)2]  •  (7.42) 

^  A  typical  value  for  Tq  in  the  visible  is  5  to  10  cm,  and  it  scales  as  The 

isoplanatic  angle  defines  a  cone  in  the  atmosphere  where  the  wavefront  is  coherent  in  an 
angular  sense.^3  To  correct  for  perturbations  in  the  propagating  medium  by  adaptive  optics 
techniques,  it  is  necessary  to  define  a  reference  wavefront  which  cannot  be  displaced  from 
^  the  outgoing  wavefront  by  more  than  Oq.  Typical  values  of  6o  are  3-10  prad  in  the  visible, 

and  it  scales  as  It  is  to  be  compared  with  the  point-ahead  angle  0pa  required  to  track 
and  point  at  fast-moving  objects.  At  geostationary  satellite  altitudes,  6pa  ~  15-20  jirad. 

Here  we  give  a  brief  discussion  of  thennal  blooming;  another  process  which 
#  become  significant  at  flux  levels  in  excess  of  IMW/cm^  is  stimulated  Raman  scattering,  in 

particular  rotational  Raman  scattering  in  N2.  For  a  discussion,  see,  e.g.  APS,  1987, 
p.SlOSf. 


7.7.2  Thermal  Blooming 

The  numerical  example  at  the  beginning  of  the  previous  subsection  points  out  that 
intense  laser  beams  propagating  through  the  atmosphere  produce  significant  heating  of  the 
medium  due  to  absorption  of  the  radiation  by  molecules  and  aerosols.  This  heating  will 
cause  the  density  and  thus  the  refractive  index  to  decrease  in  the  region  of  the  beam, 
according  to  the  relation 

(n-l)  =  kGDP  .  (7.43) 

where  the  Gladstone-Dale  constant  koD  ~  0-22  cmVg  throughout  the  IR  or  roughly  twice 
that  large  in  the  visible  (cf.,  IR  Handbook,  1985,  p.  24-16).  The  effect  of  this  change  in 
refractive  index  is  to  produce  an  aberrated  negative  lens  whose  profile  is  related  to  the  beam 
profile.  This  produces  a  deflection  and  defocusing  of  the  beam  in  a  complex  fashion  which 
is  called  "thermal  blooming."  Laser  wavelengths  for  atmospheric  propagation  are  chosen 
where  absorption  is  weak,  but  the  numerical  example  cited  above  demonstrates  that  residual 
absorption,  together  with  the  intense  radiation,  can  lead  to  blooming.  The  absorbed  energy 
per  unit  volume  is  the  product  of  the  beam  intensity,  the  absorption  coefficient,  and  a  time 
which  can  depend  on  beam  motion  through  the  air  and  heat  transport  due  to  wind  and/or 
slewing  of  the  beam. 


93 


For  turbulence-degraded  phase  fronts  measured  in  two  directions  0  and  (0  +  0o)  the  correlation 
coefficient  between  the  two  will  be  e~k 
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The  physics  of  blooming  is  reasonably  well  understood,  if  complex.  There  are  two 
limits  of  time  dependence,  namely  a  steady-state  behavior,  which  is  generally  importunr  for 
continuous  wave  laser  propagation,  and  a  transient  or  t^  behavior,  iiequently 

dominates  in  the  case  of  pulsed  laser  propagation.  For  a  discussion  of  CW  blooming,  see 
APS,  1987,  p.  S105f;  a  complementary  discussion  is  given  in  HELPH-I,  Ch.  6. 
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FORMAL  ASPECTS  OF  THE  EQUATIONS  OF 
GAS  DYNAMICS 

This  Appendix  sketches  how  one  can  derive  the  equations  of  gas  dynamics-which 
express  conservation  of  mass,  momentum,  and  energy-  from  the  Boltzmann  Equation.^  ’ 
Solutions  in  successive  powers  of  the  ratio  of  the  gas-kinetic  mean  free  path,  Igk,  to  a 
characteristic  dimension  of  the  flow,  L,  give  first  the  Euler  equations  {O(lgk/L)0}  which 
have  no  explicit  collisional  terms,  then  the  Navier-Stokes  Equations  {0(lgk/L)^}  which 
introduce  the  transport  coefficients  of  viscosity,  diffusivity,  and  thermal  conductivity, 
which  are  proportional  to  the  first  power  in  Igk-^^ 

The  distribution  function  f  =  is  a  solution  of  the  Boltzmann  equation: 

df/dt  +  (y.grad)  f  +  (X.grady)  f  =  Dc  f  ,  (E.  1) 


where  X  =  force  per  unit  mass  (which  is  here  put  equal  to  zero)  and  Dc  is  the  usual 
collision  term.  If  there  are  n  particles  (of  mass  m)  per  unit  volume,  then  the  density  p,  the 
flow  velocity  u,  and  the  temperature  T  are  defined  by  the  standard  expressions: 


p  =  m  /  f  d^vsm.n 

(E.2) 

a  =  (l/n)/fyd3y 

(E.3) 

T  =  (l/nR)/v2fd3y  =  p/Rp  , 

(E.4) 

where  the  V  is  the  velocity  relative  to  the  mean  flow: 

V=y-ii  . 

(E.5) 

The  pressure  tensor  Pij  (i,j  =  1,  2,  3)  and  the  heat  flux  vector  a  are  given  by  the 
expressions: 


See,  e.g.  Chapman  &  Cowling,  1952,  or  Hirschfelder,  Curtiss,  and  Bird,  1954,  1964,  for  a  discussion 
of  this  fundamental  equation. 

9  E-2  The  next  iteration,  with  terms  0(tgitA-)2  gives  the  Burnett  and  Grad-Maxwell  13  moment 

approximations,  but  this  is  rarely  attempted. 
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(E.6) 

(E.7) 


Pij  =  m  1  Vi  Vj  f  d^v  =  p  8ij  -  Xij 

a  =  (ni/2)/Y  V2fd3v  . 

where  ty  is  the  (symmetric)  stress  tensor,  and  conventionally 

5jj  =  lifi==j  ;0ifi=/=j  .  (E.8) 

To  obtain  the  macroscopic  flow  equations,  one  takes  the  Boltzmann  equation  (E.  1), 
multiplies  by 

=  m  ;  =  mY ;  =  (m/2)V2  +  internal  energy  ,  (E.9) 

and  integrates  over  3-dimensional  velocity  space,  Jd^v.  On  account  of  the  conservation 
laws  for  mass  (c),  momentum  (m),  and  energy  (e)  before  and  after  collisions,  the  collision 
term  vanishes: 

/  Dc  f  d^y  =  0  ;  a  =  c,m,e  ,  (E.  1 0) 

and  thus  one  gets  the  result 

0  =  3(n\|/3))/  0t  +  div  n<\j/(2)^>av  -  n  <dv/(^V3t>  av  -  n<(y.grad)\i/^)>av  - 

n  <X.gradv  >  av  ;  a  =  c,m,e  ,  (E.  1 1 ) 

where 

<K>av  =  (l/n)jKfd3y  .  (E.12) 

If  one  now  substitutes  Eq.  (E.7)  in  (E.9),  then  one  gets  the  nnacroscopic  flow  equations: 


Continuity  (c): 

d  p/d  t  E  i=i3  9(r  ui  )/9xi  =  0 

(E.l  3) 

Momentum  (m): 

d  ui/9 1  +  Ej=i^(uj  9  ui  /  9  xj  (1/  r )  9  Pij/9  xj]  =  0 

(E.14) 

Energy  (e): 

9  p/9 1  Ei=i3  [9  (p  ui)/9  xi  -h  (2/3)  9  qi/9  xi  + 
(2/3)Zj=i3pij9ui/9xj]  . 

(E.l  5) 

With  the  present  customary  form  for  the  collision  term  Dcf  in  the  Boltzmann 
equation  (E.l)  which  leads  to  the  pressure  tensor  Pij  and  the  heat  flux  vector  q,  the  system 
of  equations  (E.l 3)  to  (E.l 5)  is  incomplete  in  that  there  arc  just  five  scalar  equations  but 
many  more  than  five  unknowns  [p  ,  p,  uj  make  five,  plus  Xij  gives  five  more  (a  symmetric 
tensor  has  six  independent  components,  but  we  have  taken  out  p),  plus  qi  which  makes 
13  in  all— cf.  Hirschfelder,  Curtiss,  and  Bird  (1954),  p.  492f].  Thus  one  must  make 
appropriate  assumptions  for  Xjj  and  qj  in  the  equations  of  gas  dynamics  (E.l 3- 15),  and  this 
is  done  by  expanding  the  collision  term  Dcf  in  powers  of  the  ratio  of  gas-kinetic  mean  free 
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path,  Igk  to  a  characteristic  flow  length  L.  The  lowest  order  (Euler)  approximation  is 
obtained  by  putting  all  the  collision  terms  equal  to  zero; 

Tij(0)  =  0;qi(0)  =  0(i.j  =  1,2,3)  .  (E.16) 

The  Euler  equations  are  useful  for  explaining  fluid  flows  when  the  only  applied  forces  are 
normal,  i.e.pressure.  However,  there  is  no  friction  as  long  as  Eq.  (E.16)  applies,  and  thus 
the  fluid  cannot  cannot  sustain  shearing  forces,  such  as  the  drag  force  on  a  surface  in  the 
fluid. 

The  next  order  term  in  the  expansion,  O(lgkA-),  gives  the  Navier-Stokes 
approximation: 

Tij(i)  =  ^iaui/axj;  qi(i)=-^aT/axi  ,  (e.i?) 

which  introduces  the  viscosity  coefficient  |i  and  the  thermal  conductivity  k.  These 
quantities  depend  on  collisions  of  gas  molecules  with  one  another  and  with  surfaces,  and 
can  be  expressed  in  terms  of  the  elementary  kinetic  theory  of  gases  as  follows^-^: 

lt  =  (l/3)pvalgk  (E.18) 

k  =  fliCy  (E.19) 

D  =  (l/3)valgk  ,  (E.20) 

where  the  thermal  mean  speed  Va  and  the  sound  speed  a  are  given  as^*^ 

Vg  =  y  8kT/jtM  (E.21a) 

a  =  V7kT/M  .  (E.21b) 

Here  Cv  is  the  specific  heat  at  constant  volume  per  unit  mass,  D  is  the  diffusion  coefficient, 
Igk  is  the  gas-kinetic  mean  free  path 

.  (E.22) 

where  Ogk  is  the  gas-kinetic  collision  cross  section  [~(3-6)  x  10~^5  cm2— see,  g  g  ^ 
Section  2.5,  Eq.  (3.28)  for  another  application]  and  n  =  p/M  =  no.  of  molecules  per  unit 

The  numerical  coefficients  in  Eqs.  (E.18)  and  (E.20)  are  respectively  5ji/32  =  0.489  and  37t/16  =  0.587 
for  a  rigid  sphere,  while  f  in  Eq.  (E.19)  is  the  "Eucken  factor”  (9y-  5)/4  -  2-2.5  with  y  =  Cp/Cy  (see, 
e.g.,  Hirschfelder,  Curtiss,  and  Bird,  1954,  p.  13f0. 

Note  that  in  Eq.  (E.21)  k  is  the  Boltzmann  constant,  not  the  thermal  conductivity  which  is  here  written 
as  k. 


volume.  For  the  other  numerical  coefficients  in  Eq.  (E.21)  and  (E.22),  see,  e.g.,  Jeans 
(1925,  1982). 

The  equations  of  fluid  mechanics  are  non-linear  and  very  difficult  to  solve  in 
general.  Traditionally  aerodynamicists  use  similarity  solutions  to  a  great  extent,  with  a 
number  of  dimensionless  parameters  which  are  discussed  briefly  in  Section  5.5  and  Table 
5.3  in  the  text. 
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Here  we  present  some  details  of  the  calculation  of  reentry  vehicle  optical  signature 
from  Section  5.9.  The  geometrical  profile  of  the  RV  is  shown  in  Fig.  F.l  which  shows  the 
eight  segments  into  which  the  body  is  broken  up  for  purposes  of  computation.  The 
computation  of  the  thermal  emission  signal  in  terms  of  the  contributions  of  these  eight 
segments  is  outlined  in  Table  F.l. 
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Figure  F.1.  Simplified  Geometrical  Profile  of  an  RV 
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Table  F.l.  Thermal  Emission  Signal  Computation 


A.  To  Compute  Thermal  Radiance  (W/m2-sr>M.m)  at  the  Higher  Altitudes 


Segment  j 

1 

2 

3 

4 

5 

6 

7 

8 

Fractional  area  oj 

.011 

.029 

.080 

.12 

.16 

.20 

.24 

.17 

h  =  3QQkft.  (91  kmi 

Tj(K) 

810 

560 

477 

422 

394 

377 - 

— 

— 

B(0.5.Tj) 

1.4E-6 

2E-13 

2E-17 

9E-21 

7E-23 

3E-24 

B(3.5.Tj) 

1400 

150 

41 

13 

7 

4 

B(12,Tj) 

140 

64 

42 

30 

24 

20 

aj  B(>..Tj) 

aj  B(0.5,Tj) 

1.6E-8 

5E-15 

2E-18 

IE-21 

IE-23 

2E-24 

1.6E-8 

ttj  B(3.5Jj) 

16 

4 

3 

2 

1 

3 

29 

ajB(12,T|) 

2 

2 

3 

4 

4 

13 

27 

h  =  15Qkn.  (46  km) 

Tj(K) 

2590 

2200 

1420 

1090 

920 

870 

810 

780 

B(0.5,Tj) 

5.7  E4 

8  E3 

6 

1  E-2 

1  E-4 

2E-5 

1  E-6 

4E-7 

B(3,5.Tj) 

5.8  E4 

4.1  E4 

1.3  E4 

5.3  E3 

2.6  E3 

2.0  E3 

1.4  E3 

1.2  E3 

B(12.Tj) 

810 

660 

360 

240 

180 

160 

140 

130 

aj  B(0.5.Tj) 

660 

230 

0.5 

- 

*  900 

Oj  B(3.5,Tj) 

670 

1210 

1070 

640 

420 

400 

340 

190  *  5000 

ajB(12,Tj) 

9 

19 

29 

29 

28 

32 

33 

22  •  200 
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Table  F.1  (continued) 


B.  Thermal  Emission  Signal  for  Different  Conditions 


m 

h  (kft) 

300 

150 

75 

25 

(km) 

91 

46 

23 

7.6 

T  (K) 

Variable- 

see  A,  above 

2590 

700 

m 

•  =  0.5  pm:  Zj  aj  B(0.5,Tj) 

1.6  E-8 

900 

57000 

5.3  E-9 

*e 

LT.e 

1.5  E-8 

840 

53000 

4.9E-9 

m 

•  X  =  3.5  pm:  Zj  aj  B(3.5,Tj)* 

29 

5000 

58000 

640 

LT.B 

38 

6500 

76000 

840 

• 

•?t=12pm:  ZjOjB(12.Tj)* 

27 

200 

810 

110 

LT.e 

45 

340 

1400 

180 

* 


W/m^-sr-^im 

W/sr-nm,  i.e.,  we  have  multiplied  Z  a  B(X,T)  by  e  Ap.  £  from  Table  5.4  and  Ap 
Eq.  (5.24). 
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THE  ILLUMINATION  OF  A  FLAT  PLATE 
BY  EARTHSHINE 


G-1 


APPENDIX  G 

THE  ILLUMINATION  OF  A  FLAT  PLATE 
BY  EARTHSHINE 


The  problem  for  a  flat  plate  is  defined  in  Fig.  5.7"See  also  Eqs.  (5.29)  and  (5.30). 
Because  of  its  complexity,  I  display  how  one  evaluates  the  integral  over  energy  radiated 
from  various  elements  of  the  earth’s  surface: 

E*T£  =  1  dAE  £Ea)  B(X,Te)  cos  (Pi  [RE(0a)]-2  ■  (G.  1 ) 


After  diffuse  reflection  from  the  target,  the  effective  radiance  seen  by  the  sensor  is 

Lt£  =  E*t£  cos()>2  (p/7t)  .  (G.2) 

To  evaluate  the  integral  in  Eq.(G.l),  change  the  coordinate  01  to  (0a  +  P)-'See 
Fig.  G.l: 

01  =  0a  +  P 


s  =  Re  tan  0a ;  ds  =  Re  tan20a  d0a  (G.3) 

dAE  =  2  7t  Re^  tan  0a  sec20a  d0a 
Re=  Re  sec0a ;  teC^a)  =  '^E  sec0a 

and  thus 

E*t£  =  2  7t  e(X)  B(X,Te)  It(Xe.P)  >  (G.4) 


where 

It(xe.P)  =  {IqP  +  (1/2)  tan0a  cos(0a  +  p)  e'b^^a  d0a  .  (G.5) 

Here, 

b  =  ln(l/XE)  ,  (G.6) 

and  the  region  of  integration  is  split  because  a  portion  of  the  (27i)  steradian  seen  from  the 
earth  is  screened  at  an  angle  of  attack  P  ^  0.  It(xe>P)  is  well  behaved  and  falls  off  rapidly 
with  increasing  0(x  unless  b  =  0  ,  i.e.,  unless  te  =  1  • 
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the  Evaluation  of  Eq.  (G.1) 
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Figure  G.1.  Sketch  for 


I  have  evaluated  iT('tE.P)  for  the  representative  case  j3  =  7i/4,  when  Eq.  (G.5) 
becomes 

lT('rE.45°)  =  tan  <()  (cos  <))  -  sin  <}> )  d0/2t(2  .  (G.7) 

Numerical  values  for  IjCh)  are  shown  in  Fig.  G.2  for  P  =  45°. 


Figure  G.2. 


The  Integral  1(b) 
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ELECTROMAGNETIC-RADIATION  EFFECTS  ON 
ELECTRICAL  AND  ELECTRONIC  SYSTEMS 
(J.  Malik,  LANL,  from  Hoerlin,  1976) 


APPENDIX  H 

ELECTROMAGNETIC-RADIATION  EFFECTS  ON 
ELECTRICAL  AND  ELECTRONIC  SYSTEMS 
(J.  Malik,  LANL,  from  Hoerlin,  1976) 

The  electromagnetic  radiation  in  the  radio-frequency  portion  of  the  spectrum 
(EMP)  can  cause  problems  in  electronic  systems.  The  pulse  from 
detonations  above  about  30  km  is  caused  by  the  deflection  of  Compton 
electrons  produced  by  gamma-ray  interaction  with  the  earth's  magnetic  field 
in  the  deposition  region  (20-40  km).  The  resulting  transverse  current  in  the 
large  area  of  gamma-ray  deposition  produces  a  large  coherent  radiating 
element.  With  appropriate  yield,  detonation  altitude  and  magnetic  azimuth, 
the  electric  fields  over  large  areas  at  the  earth's  surface  can  exceed  10^  V/m. 
Such  fields  can  cause  detrimental  effects  on  some  types  of  electrical 
systems.  The  pulse  width  is  less  than  a  microsecond. 

Starfish  produced  the  largest  fields  of  the  high-altitude  detonations;  they 
caused  outages  of  the  series-connected  street-lighting  systems  of  Oahu 
(Hawaii),  probable  failure  of  a  microwave  repoeating  station  on  Kauai, 
failure  of  the  input  stages  of  ionospheric  sounders  and  damage  to  rectifiers 
in  communication  receivers.  Other  than  the  failure  of  the  microwave  link, 
no  problem  was  noted  in  the  telephone  system.  No  failure  was  noted  in  the 
telemetry  system  used  for  data  transmission  on  board  the  many 
instrumentation  rockets. 

There  was  no  apparent  increase  in  radio  or  television  vepairs  subsequent  to 
any  of  the  JI  detonations.  The  failures  were  generally  in  the  unprotected 
input  stages  of  receivers  or  in  rectifiers  of  electronic  equipment;  transients 
on  the  power  line  probably  caused  the  rectifier  failures.  There  was  one 
failure  in  the  unprotected  part  of  an  electronic  system  of  the  LASL  Optical 
Station  on  top  of  Mount  Haleakala  on  Maui  Island.  With  the  increase  of 
solid-state  circuitry  over  the  vacuum-tube  technology  of  1962,  the 
susceptibility  of  electronic  equipment  will  be  higher,  and  the  probability  of 
more  problems  for  future  detonations  will  be  greater.  However,  if 
detonations  are  below  line-of-sight,  the  fields  and  therefore  system 
problems  will  be  much  smaller. 


APPENDIX  I 

SCATTERING  OF  SUNLIGHT 
FROM  A  CLOUD 


I.l  INTRODUCTION 

For  the  development  of  a  space-based  surveillance  system  in  the  atmospheric 
2.7  p.m  absorption  band  the  issue  of  false  alarms  due  to  the  scattering  of  sunlight  from 
clouds  was  addressed.  An  extensive  experimental  program  using  instruments  mounted  on 
high-altitude  balloons  and  on  U-2  high-altitude  aircraft  was  conducted  by  Baird- Atomic- 
see  Elliott,  1961,  and  Bauer,  1963, 1964,  also  Espinola,  1970.^-' 

The  physics  is  quite  different,  depending  on  the  wavelength  between  visible  and 
LWIR  and  on  the  effective  scattering  angle.  At  wavelengths  <  4  p.m  a  target  can  be  seen 
principally  because  it  scatters  sunlight,  while  at  the  longer  wavelengths  (and  at  night)  the 
principal  contribution  to  the  target  signature  comes  from  thermal  emission  from  the  body, 
with  a  relatively  minor  contribution  from  the  scattering  of  earthshine  (thermal  emission 
from  the  surface  of  the  earth  and  lower  atmosphere).  For  large  scattering  angles,^'^  which 
is  the  case  considered  here,  single  scattering  will  predominate,^'^  while  for  small  scattering 
angles  (near-forward  scattering)  there  are  multiple  scattering  effects.  Also,  the  size  of  the 
cloud  can  make  a  significant  difference  to  its  signature.  Here  we  discuss  some  of  these 
factors  in  a  preliminary  way,  namely: 

a.  Small  versus  large  clouds 

b.  Geometrical  factors 

c.  Effective  diffuse  reflectivity. 


For  more  recent  data  and  analyses,  see  Kisliuk,  1979,  and  Humphrey,  1984. 

^■2  The  scattering  angle  0  is  so  defined  that  0  =  0  corresponds  to  forward  scattering,  0  =  180  deg  to  back- 
scauering.  See  Fig.  I- 1. 

For  an  empirical  analysis  of  multiple  scattering  effects  in  this  situation,  see  Young,  1978. 


for  the  application  of  reflection  rather  than  transmittance  of  sunlight.  The  scenario 
considered  here  is  discussed  in  Section  J.3,  but  first  we  distinguish  between  small  and 
large  clouds. 

1.2  SMALL  VERSUS  LARGE  CLOUDS 

Given  the  scattering  and  extinction  properties  of  cloud  particles,  one  must 
distinguish  between  a  small  cloud  in  which  all  particles  scatter  the  incident  radiation,  and  a 
large  cloud  in  which  the  incident  radiation  is  extincted  within  the  cloud  so  that  not  all 
particles  serve  as  scatterers  of  the  incident  radiation. 

For  a  finite  ("cubical")  cloud  of  linear  dimension  L  made  up  of  n  particles  per  unit 
volume,  each  of  extinction/scattering  cross  sections  Oext.  <7sca  we  define  an  extinction 
distance  where 

p  =  n  Oext  •  (I-l) 

Consider  a  small  cloud  of  dimension  L,  which  is  made  up  of  L^n  (independent)  scatterers, 
so  that  its  total  scattering  cross  section  is 

Q(small;0)  =  L^n  (dasca/^^)6  .  (1-2) 

where  the  mean  differential  cross  section  of  a  cloud  droplet  is  given  from  Fig.  7.5.  By 
contrast,  for  a  Jorge  cloud,  not  all  L^n  particles  are  available  to  provide  effective  scattering 
because  the  interior  panicles  of  the  optically  thick  cloud,  i.e.,  those  which  are  more  than  a 
distance  P“^  from  the  surface,  are  shielded  from  the  incident  solar  radiation.  Thus,  the 
effective  total  number  of  scatterers  is  L^P"!  n  ,  and  so  the  total  scattering  cross  section  for 
a  large  cloud  is 

QGarge;e)  =  L2p-1  n  (dasca/dQ)0  =  (L^/Oext)  (dOgca/^^Je  ■  (L3) 

Equation  (1.3)  is  only  schematic  in  that  there  are  some  geometrical  (angular)  factors  that 
have  not  been  included  here. 

1.3  GEOMETRICAL  FACTORS 

Figure  1. 1  shows  the  geometry  of  the  scattering  process.  We  consider  the  coplanar 
problem  with  all  of  SA,  DA,  NA  in  the  same  plane  so  that 

i  +  e  +  0  =  7t  (I.4a) 

and  also  the  condition  of  specular  reflectance 
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i  =  e  (1. 4b) 

which  give  the  maximum  amount  of  scattering  for  the  case  of  specular  reflection  (which  is  ( 

not  necessarily  the  case  here). 

It  is  clear  that  the  actual  scattering  cross  section  Q  (large;  9)  of  Eq.  (1.3)  actually 
depends  on  the  macroscopic  angles  i,  e,  which  determine  the  relation  between  the  incident 
and  reflected  light  beams  relative  to  the  macroscopic  cloud  as  well  as  on  the  scattering  angle  * 

6  for  a  particular  particle,  as  represented  by  Fig.  7.4,  7.5,  or  1-2.  It  is  useful  to  present  a 
simple  parameterization  in  terms  of  a  Lambert's  Law  diffuse  reflection  coefficient  ri(6), 
even  though,  as  is  pointed  out  in  Bauer  1964,  the  description  of  a  cloud  as  a  Lambertian 
reflector  is  very  poor,  except  possibly  in  the  visible  where  there  is  a  great  deal  of  multiple  ^ 

scattering  at  small  values  of  6  so  that  the  single-scattering  .  alysis  is  inappropriate. 

We  can  define  by  the  following  expression  (from  Bauer,  1 964'-5) 


flmax(0) 

=  7C  [u  (i.e.)  OexJ"^ 

(1.5) 

u  (i.e.) 

=  2  cos  i  cos  e/sin(0/2) 

(L6) 

i 

=  e  =  (7i-e)/2 

a4c) 

Results  at  several  different  wavelengths  are  given  in  Fig.  1.2;  note  that  for  0  <  30  deg, 
T)  >  1  which  clearly  indicates  that  the  formula  for  Tiinax(0)  here  is  inappropriate.  At 
larger  angles  the  reflectivity  is  fairly  small.*'^  The  physics  of  this  is  fairly  easy  to 
understand  [cf.,  e.g.,  Bauer,  1964,  below  Eq.  (3)].  Because  the  complex  refractive  index 
N 1  -  iN2  of  the  water  droplets  is  not  very  different  from  one,  the  reflectivity  at  normal 
incidence  from  a  plane  slab,  which  is  given  by  the  expression 

R  =  [(Ni-l)2-»-N22]/[(Ni-H)2-f  N22]  (1.7) 

is  quite  small,  comparable  with  the  relative  energy  scattered  through  angles  greater  than 
90  deg  (see  Table  I.l). 


Young,  1978,  presented  an  empirical  correction  for  multiple  scattering  effects. 

It  must  be  stressed  that  the  expression  for  u  is  not  a  good  one,  although  clearly  something  of  this  type 
is  needed.  There  are  problems  both  with  the  geometry  and  also  with  effects  of  single  versus  mutliple 
scattering.  Readers  are  solicited  for  a  better  expression! 

I-®  The  increased  reflectivity  in  the  visible  for  9  >  120  deg  is  due  to  rainbow  backscattering-see  Fig.  7.4 
7.54,  or  1-2. 
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Figure  1.2.  Scattering  From  a  Large  Cloud  at  Different  Wavelengths: 
Effective  Diffuse  Reflectivity 


Table  1.1.  Scattering  Parameters  at  Different  Wavelengths 


Wavelength  (urn) 

0.7 

2.7 

10 

Qav  ~  27taav^7.  (Sav  ~  4p.m) 

36 

9.3 

2.5 

Osca*  (10"®  cm2) 

1.45 

1.38 

Oabs*  (10-®cm2) 

- 

1.41 

CText**  (10-®  cm2) 

1.45 

1.45 

1.45 

m  «  Ni  -  iN2*" 

1.331-3.05E-8 

1.188-1 .90E-2 

1. 218-5. 08E-2 

R  [from  Eq.(A.7)] 

0.020 

0.0074 

0.010 

*  Bauer  et  al.,  1984,  p.  3-4--ffom  E.  Shettle,  AFGL. 

* '  "oext  generally  agrees  with  2  Ogeom  "  1  AS  to  within  a  factor  of  2  (Bauer  et  al.,  1 984,  p.  3-4). 
USAF,  1985,  p.  18-16  (from  E,  Shehle). 
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The  effective  reflectivity  R  due  to  a  slab  of  water  (as  against  an  assembly  of 
particles)  can  be  compared  roughly  with  Tj  (0  =  180  deg).  One  would  expect  the  actual 
reflectivity  due  to  the  particles  to  be  smaller  than  R,  especially  in  the  LWIR  where  the  Mie 
scattering  parameter  qav  for  mean  particle  radius  aav 

Qav  ~  2  71  aavA.  (1-8) 

is  close  to  one,  since  if  q  «1,  there  is  very  little  scattering. 

It  must  be  stressed  that  the  significance  of  the  effective  diffusivity  is  (]uite  different 
at  varying  wavelengths.  In  the  visible  a  target  is  seen  by  solar  scattering,  and  at 
X  =  2.7  )im  the  solar  scattering  is  still  greater  than  thermal  emission  in  the  daytime;  but  at 
X,  =  10  |im  a  cloud  (or  other  target)  is  seen  by  its  thermal  emission  supplemented  by  the 
scattering  of  earthshine,  rather  than  by  the  scattering  of  sunlight 
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